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AUTOMATED  calculation  of  protection  factors 

FOR  THE  SODIUM  CHLORIDE  RESPIRATOR 
QUANTITATIVE  FIT  TEST  INSTRUMENT 


INTRODUCTION 

The  purpose  of  this  report  is  to  present  an  automated  procedure  for  cal¬ 
culating  a  respirator's  protection  factor  (PF)  afforded  to  the  respiratory 
tract  and  eyes  against  chemical  warfare  (CW)  agents  in  particulate,  aerosol, 
or  vapor  form.  The  participating  nations  of  the  Air  Standardization  Coordi¬ 
nating  Committee  (ASCC)  have  drafted  a  North  Atlantic  Treaty  Organization 
(NATO)  Standardization  Agreement  (STANAG)  recommending  a  sodium  chloride  quan 
titative  fit  test  scheme  for  this  procedure  [1].  The  Chemical  Defense  Estab¬ 
lishment  (CDE)  in  the  Ministry  of  Defense  (MOD)  of  the  United  Kingdom  is 
credited  with  developing  the  basic  sodium  chloride  (NaCl)  quantitative  fit 
test  technology  [2-11].  The  United  States  Air  Force  School  of  Aerospace  Medi 
cine  (USAFSAM),  which  has  had  several  years  of  laboratory  experience  with  the 
sodium  chloride  quantitative  fit  test  instrumentation  [1],  has  reviewed  and 
endorsed  this  draft  agreement. 

Our  experience  with  the  sodium  chloride  quantitative  fit  test  instrumen¬ 
tation  has  shown  one  area  to  be  of  particular  concern;  namely,  the  method  of 
reducing  the  collected  data  (respirator  sodium  chloride  penetration  concentra 
tion  for  a  particular  exercise  protocol)  and  calculating  a  protection  factor. 
Tht.  STAN AG  agreement  draft  proposes  that: 

'"Lon  the  output  of  the  flame  photometer  fluctuates  during 
o  measurement  of  penetration,  the  maximum  output  is  to  be 
used  to  calculate  the  Protection  Factor.  For  this  pur¬ 
pose,  occasional  transient  increases  in  output  the  dura¬ 
tion  of  each  of  which  does  not  exceed  2  sec  may  be 
ignored.”  [1] 

In  addition,  various  interested  organizations  have  developed  and  reported  on 
the  following  data  reduction  schemes  for  respirator  quantitative  fit  testing 
[2.13-43]: 

a.  selection  of  the  overall  maximum  output  peak 

b.  arithmetic  average  of  the  maximum  output  peaks 

c.  arithmetic  average  of  the  maximum  output  peaks  and  minimum 
valleys  (midpoint) 

d.  visual  estimation  of  the  midpoint  between  the  maximum  output 
peaks  and  minimum  valleys 

e.  time-averaged  or  integrated  value. 


This  report  develops,  with  some  rigor,  an  automated  procedure  to  reduce 
the  NaCl  leak  test  data  and  calculate  a  protection  factor.  An  overview  of  the 
USAFSAM  sodium  chloride  respirator  quantitative  fit  test  instrument  is  fol¬ 
lowed  by:  a  discussion  on  conventional  protection  factor  calculations;  the 
use  of  a  voltage- to- frequency  (V/F)  converter  circuit  to  do  time-averaged 
integration;  and,  finally,  a  least-squares  curve  fit  computer  program  to  cal¬ 
culate  a  protection  factor. 


USAFSAM  SODIUM  CHLORIDE  RESPIRATOR 
QUANTITATIVE  FIT  TEST  INSTRUMENT 

The  CDE  sodium  chloride  respirator  quantitative  fit  test  method,  de¬ 
scribed  in  British  Standards  4400  and  2091,  has  been  adapted  and  modified  by 
USAFSAM  to  measure  the  protection  factor  in  the  respiratory  and  eye  compart¬ 
ments  of  aircrew  chemical  defense  respirators  [3.11J.  This  instrument  gener¬ 
ates  a  solid  aerosol  of  sodium  chloride  crystals  as  the  challenge  atmosphere. 
The  concentration  of  the  challenge  atmosphere  in  an  aircrew  respirator  is 
measured  using  a  hydrogen  flame  photometer,  and  the  result  is  displayed  on  a 
strip-chart  recorder.  This  technique  allows  protection  factors  as  high  as 
106  to  be  calculated. 


Instrument  Description 

Illustrated  in  Figure  1  are  the  primary  components  used  in  the  USAFSAM 
sodium  chloride  respirator  quantitative  fit  test  instrument  [1,3,11,33].  The 
sodium  chloride  solid  aerosol  challenge  atmosphere  is  generated  hy  atomizing  a 
s oo in"  chloride  solution,  drying  the  liquid  aerosol  in  a  drying  tube,  and 
delivering  the  dry  cloud  to  the  top  of  a  transparent  plastic  hood.  A  subject, 
having  donned  a  respirator,  enters  the  hood  and  performs  a  series  of  breathing 
and  head  movement  exercises.  A  sampling  pump  is  used  to  draw  a  portion  of  the 
atmosphere  from  the  interior  compartment  of  the  respirator.  This  sample  is 
vaporized  in  a  hydrogen  flame  photometer,  and  the  output  signal  is  displayed 
on  a  strip-chart  recorder. 

Production  of  the  Sodium  Chloride  Challenge  Atmosphere 

The  challenge  atmosphere  is  generated  by  atomizing  a  5t  aqueous  sodium 
chloride  solution,  prepared  by  dissolving  SO. 00  grams  of  sodium  chloride  in 
1000  mt  of  distilled  water.  Operating  a  Oautrebande  atomizer  at  70  psi 
(482.3  kPa)  (clean  dry  compressed  air)  yields  a  flow  through  the  atomizer  of 
35  liters/ttdn  (STP).  The  liquid  aerosol  from  the  atomizer  is  then  injected 
perpendicularly  into  the  air  stream  flowing  through  a  mixing  and  drying  tube. 
The  50  liter/min  (STP)  source  of  clean  dry  air  evaporates  the  water  from  the 
liquid  sodium  chloride  aerosol  and  produces  a  solid  sodium  chloride  aerosol 
challenge  atmosphere.  The  mass  median  aerodynamic  diameter  (NMAD)  of  the  dry 
sodium  chloride  crystals  range  from  0.4  -  0.6  urn.  The  sodium  chloride  chal¬ 
lenge  atmosphere  is  delivered  through  a  short  length  of  Tygon  tubing  to  the 
top  center  of  a  transparent  plastic  hood  that  covers  the  subject  from  head  to 
waist.  The  hood's  vertical  displacement  is  controlled  by  releasing  a  locking 
pin  and  turning  a  cable-connected  crank.  The  diameter  of  the  hood  permits  a 
subject  to  move  freely  his  arms,  shoulders,  and  head.  An  adjustable  cloth  col¬ 
lar  attached  to  the  bottom  edge  of  the  hood  is  drawn  snugly  around  the  sub¬ 
ject's  waist  and  serves  to  contain  the  challenge  atmosphere  [1,3,11,33]. 
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Figure  1*  USAFSAM  sodium  chloride  respirator  quantitative  fit  test  instrument. 


Measurement  of  the  Sodium  Chloride  Respirator  Leakage 

When  evaluating  a  full-face  aircrew  chemical  defense  respirator,  a  pri¬ 
mary  concern  is  the  penetration  of  the  sodium  chloride  challenge  atmosphere 
into  the  visual  compartment .  In  order  to  make  this  measurement,  an  aluminum 
tube,  approximately  1  in.  (2.54  cm)  long  and  0.25  in.  (0.635  cm)  i.d.,  is  so 
fitted  through  and  sealed  to  the  protective  device's  visor  that  the  distance 
from  the  cornea  to  the  open  end  of  the  aluminum  sampling  tube  (interior  to  the 
respirator's  visor)  is  not  greater  than  0.8  in.  (2  cm)  [1]. 

The  concentration  of  the  sodium  chloride  challenge  atmosphere  that  has 
leaked  into  the  visual  compartment  is  determined  by  continuously  sampling  gas 
from  this  site  and  vaporizing  it  in  the  flame  photometer.  Sampling  is  accom¬ 
plished  by  attaching  one  end  of  a  short  length  of  Tygon  tubing  to  the  open  end 
of  the  aluminum  sampling  tube  (exterior  to  the  respirator's  visor),  and  then 
passing  the  opposite  end  of  the  plastic  tubing  through  a  sealed  port  in  the 
top  of  the  hood.  A  pump--of  the  metal  bellows  positive  displacement,  continu¬ 
ously  sampling  type--is  connected  to  the  open  end  of  the  hood's  port,  and  is 
used  to  draw  a  gas  sample  from  the  respirator's  visual  compartment  (constant 
7-1  iter /min  flow).  Before  the  7-1  iter /min  sample  is  injected  into  the  hydro¬ 
gen  flame  photometer,  a  6-1  iter /min  amount  is  bled  to  the  ambient  atmosphere 
through  a  calibrated  orifice.  The  resulting  1-liter/min  sample  is  diluted 
with  a  14-liter /min  flow  of  clean  dry  compressed  air,  and  this  mixture  is  then 
injected  into  the  photometer  for  analysis.  The  determination  of  the  7-liter/ 
min  mask  sampling  rate,  6-1  iter /min  bleed-off,  and  subsequent  14-1  iter /min 
dilution  of  the  1-liter/min  portion  of  the  penetration  sample  was  based  on  two 
experimental  observations:  First,  the  7-liter /min  sampling  rate  was  selected 
to  minimize  the  negative  pressure  within  the  respirator's  visual  compartment; 
i.e.,  a  greater  sampling  rate  was  observed  to  distort  the  penetration  measure¬ 
ment.  Second,  the  sodium  chloride-air  mixture  composition  was  metered  to 
yield  optimum  performance  of  the  hydrogen  flame  photometer  for  PF's  ranging 
from  l(r  to  10°  [33]. 

The  flame  photometer  used  to  analyze  the  sampled  sodium  chloride  aerosol 
from  the  interior  of  a  respirator  consisted  of  two  primary  components:  a 
burner  to  vaporize  the  sodium  chloride  crystals,  and  an  electronic  package  to 
quantify  the  actual  leakage  concentration. 

Hydrogen  is  used  as  the  combustion  gas  for  the  burner  because  of  its 
pure,  almost  colorless  (pale  blue)  flame.  With  the  particular  gas  jet  being 
used,  the  hydrogen  source  is  regulated  to  flow  at,  a  nominal  0.05  -  0.07  psi 
(0.34  -  0.48  kPa);  these  conditions  produce  a  1-in.  (2.54  cm)  vertical  flame 
[33]. 


The  electro-optic  scheme  used  to  quantify  a  respirator's  sodium  chloride 
leakage  is  shown  in  Figure  2  [33].  In  operation,  the  sampled  sodium  chloride 
crystals  are  vaporized  in  the  hydrogen  flame  photometer.  The  589  nanometer 
(nm)  wavelength  optical  bandpass  filter  is  used  to  detect  the  yellow  sodium 
emission  lines  and  to  reject  all  other  undesirable  light  energy  emissions. 
Detection  and  quantification  of  the  intensity  of  the  filtered  yellow  light  is 
accomplished  with  a  photomultiplier  tube  (PMT).  In  principle,  the  PMT  output 
current  Is  directly  proportional  to  the  Intensity  of  the  yellow  light 
impinging  upon  Us  cathode.  Thus,  since  the  Intensity  of  yellow  light  pro¬ 
duced  by  the  hydrogen  flame  Is  proportional  to  the  concentration  of  sodium 
chloride  In  the  medium  surrounding  the  flame,  the  PMT  output  current  Is 
also  proportional  to  the  sodium  chloride  concentration  around  the  flame. 
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processing  scheme  used  in  the  USAFSAM  sodium  chloride  respirator 


Electronic  processing  of  the  PMT  output  current  is  necessary  for  three 
reasons: 

a.  The  PMT  output  current  must  be  converted  to  a  voltage  of  suffi¬ 
cient  magnitude  so  that  the  leakage  can  be  displayed  on  a  strip- 
chart  recorder. 

b.  The  PMT  output  current  must  be  filtered  to  reduce  the  random 
high-frequency  noise  components. 

c.  The  PMT  output  current  must  he  scaled  logarithmically  to  accom¬ 
modate  a  wide  dynamic  range  (four  orders  of  current  magnitude) 
on  a  single  span  strip-chart  recorder. 


The  electronic  processing  scheme  used  in  the  USAFSAM  sodium  chloride  respira¬ 
tor  quantitative  fit  test  instrument  is  shown  in  Figure  3  [33],  A  Harris 
Semiconductor  Corporation  (HA2-2905-5)  integrated  circuit  is  used  as  the 
active  element  in  the  current-to-voltage  converter;  an  Analog  Deviaes  Corpora¬ 
tion  (Model  755N)  integrated  circuit  is  used  as  the  logarithmic  amplifier;  and 
a  Signetics  Corporation  general  purpose  741  operational  amplifier  integrated 
circuit  is  used  in  the  interface  circuit  to  drive  the  strip-chart  recorder. 

The  low- pass  filter  and  logarithmic  amplifier  are  shown  in  Figure  4  [33]. 
Illustrated  in  Figure  5  are  the  electronics  associated  with  the  strip-chart 
recorder  interface  circuit  [33].  The  PMT  detector  is  an  International  Tele¬ 
phone  and  Telegraph  (ITT)  special  purpose  16-stage,  electrostatically  focused, 
model  FW130  tube.  With  an  operating  anode-to-cathode  potential  of  1300-2200 
volts,  the  PMT's  dark  current  is  100  times  less  than  the  PMT  background  cur¬ 
rent  (PMT  background  current  is  defined  to  be  that  current  produced  by  the 
tube  when  the  hydrogen  flame  is  surrounded  by  a  medium  free  of  sodium  chlo¬ 
ride).  With  the  associated  electronics  package  (Figs.  4  and  5),  the  following 
transfer  function  relates  the  PMT  output  current  to  the  strip-chart  recorder 
voltage: 


v  =  2  t  log10 


I 


_i _ 

1.1  yA 


} 


(1) 


whpre 


v  =  output  voltage  to  strip-chart  recorder  (in  volts),  and 

i  =  PMT  output  current  (in  microamperes) 

In  addition,  this  transfer  function  Is  Illustrated  graphically  In  Figure  6 
[33].  Of  particular  attention  is  the  fact  that  the  PMT  output  current  spans 
11  nanoamperes  (nA)  to  110  microamperes  (yA),  and  the  associated  strip-chart 
recorder  voltage  spans  0  to  4  volts;  l.e.,  a  1-volt  output  change  occurs  per 
decade  of  input  current  change. 

In  order  to  calculate  the  protection  factor  for  a  subject’s  chemical 
defense  respirator,  the  foregoing  Information  concerning  the  dynamics  of  the 
instrument  must  be  coupled  with  one  more  factor— the  calibration  procedure. 
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Figure  6.  Transfer  function  plot  relating  PMT  output  current  to  strip-chart  recorder  voltage. 
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Calibration  Procedure 

The  calibration  procedure  being  used  for  the  USAFSAM  sodium  chloride 
respirator  quantitative  fit  test  instrument  utilizes  .the  serial  dilution  tech¬ 
nique  [1].  This  technique  is  based  on  the  fact  that  the  concentration  of 
sodiun  chloride  in  the  challenge  atmosphere  is  directly  proportional  to  the 
concentration  of  the  sodium  chloride  in  the  atomizer  solution.  Thus,  since  a 
5%  sodiun  chloride  solution  standard  is  used  to  produce  the  challenge  atmo¬ 
sphere,  a  series  of  dilutions  of  this  standard  can  be  used  to  produce  known 
PMT  voltage  responses.  If  the  serial  dilutions  are  carefully  selected,  a 
series  of  PMT  voltages,  directly  proportional  to  these  concentrations,  can  be 
measured.  In  the  end,  it  will  be  possible  to  interpolate  between  these  known 
PMT-voltage  response  vs.  serial-dilution  concentrations,  and  to  calculate  the 
sodium  chloride  leakage  (concentration)  (and  thus  PF)  based  on  the  associated 
strip-chart  recorder  voltage  fluctuations.  The  judicious  selection  of  serial 
dilutions  was  made  to  produce  calibration  standards  ranging  from  10"5  to  10*1 
of  the  basic  5%  sodium  chloride  challenge  atmosphere  solution.  Table  1 
depicts  the  relative  concentrations  and  component  amounts  [33]  for  the  chal¬ 
lenge/calibration  sodium  chloride  solutions  (see  “Author's  Note,"  below). 


TABLE  1.  SODIUM  CHLORIDE  CHALLENGE/CALIBRATION  SOLUTIONS  [Ref.  33] 

Relative  Mass  of  sodium  Volume  of  distilled 

concentration  chloride  (grams)  water  (liters) 


10°  (5S  challenge) 

10-1 

50.00 

1.0 

5.000 

1.0 

10-2 

0.5000 

1.0 

10-3 

0,05000 

1.0 

10-4 

0.005000 

1.0 

10-5 

0.0005000 

1.0 

10-5 

0.00005000 

1.0 

When  calibrating  the  instrument,  to  avoid  contaminating  a  weaker  solution 
by  a  stronger  one,  the  operator  begins  with  an  atomizer  containing  the  10'6 
solution  and  advances  to  the  10°  solution.  The  strip-chart  voltage  amplitude 
for  each  calibration  sample  is  annotated  and  retained  for  subsequent  analy¬ 
sis.  A  typical  plot  (semilogarithmfc)  of  the  PMT  output  voltage  response  for 
each  of  the  sodium  chloride  concentration  standards  is  illustrated  in  Fig¬ 
ure  7;  and  Table  2  depicts  this  relationship  with  actual  data. 


AUTHOR'S  NOTE:  The  dual  function  of  the  10°  NaCl  solution  is  to  serve  as 
a  calibration  solution  and  the  challenge  solution. 


17 


Sodium  chloride  calibration/challenge  concentration  vs.  PMT  output  voltage  (semilogarithmic 


TABLE  2.  SODIUM  CHLORIDE  CHALLENGE /CAL I BRAT I ON  SOLUTION  CONCENTRATION  VS. 
PMT  OUTPUT  VOLTAGE 


Sodium  chloride  challenge/ 
calibration  solution  concentration 


PMT  output  voltage 
(volts) 


10°  (5%  challenge) 

3.355 

io-i 

2.915 

10-2 

2.310 

10“3 

1.500 

10‘4 

0.545 

10“5 

0.  j.55 

io-6 

0.105 

CONVENTIONAL  PROTECTION  FACTOR  CALCULATIONS 

A  general  discussion  of  a  protection  factor  is  presented  here,  along  with 
a  description  of  the  conventional  method  (hand  calculation)  used  to  calculate 
the  PF's  associated  with  aircrew  chemical  defense  respirator  quantitative  fit 
testing.  Through  this  information,  the  reader  can  evaluate  the  respective 
advantages  and  disadvantages  of  the  conventional  and  automated  methods  of 
calculating  PF's. 


Protection  Factor 

A  respiratory  protection  factor  is  defined  as  the  ratio  of  the  ambient 
challenge  atmosphere  concentration  external  to  the  respiratory  protective 
device  to  that  of  the  sampled  leakage  concentration  drawn  from  the  interior  of 
the  device  [1,2,8,11-20,26,32,34,37,41-43].  Formally,  this  relationship  can 
be  expressed  in  mathematical  terms: 

PF  «  ~  (2> 

where 

PF  =  protection  factor 

Ca  ®  ambient  challenge  atmosphere  concentration 
Cs  3  sampled  leakage  concentration 

Note  that  a  PF  is  a  dimensionless  quantity.  In  the  ratio,  the  units  of  con¬ 
centration  in  the  numerator  and  denominator  cancel  (assuming  that  Ca  and 
Cs  were  measured  and  appropriately  converted  to  a  consistent  set  of 
concentration  units;  e.g.,  parts  per  million,  micrograms  per  liter,  percent, 
etc.). 


Also  important  in  respirator  quantitative  fit  testing  is  the  calculation 
of  an  average  protection  factor  (PF).  This  calculation  becomes  important  when 
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the  subject  being  evaluated  performs  a  series  of  breathi™  k  j 

is  desisned  to  stress  the 
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i=l  PF-j 


where 


PF  =  average  protection  factor  for  n  exercises 
i  =  the  ith  exercise,  i  =  1,  2,  3,  ....  n 

PF  =  protection  factor  associated  with  a  particular  exercise 

gro&'Sr  ’r  ca"  be  ca'cu,ate,i  “he" 

exercise  PF's.  The  most  common  examnlp  P°rJ?nfe.are  ^signed  to  individual 
exercise  protocol  is  perfoS  fo^  a  ^  each  exercise  in  a" 

time  would  become  theweiqhtina  faetnr  f^ent  of  time;  in  this  case, 

expression  for  an  averageweJShted  PF is:  COmpleteness*  *  mathematical 
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where 


PFW  »  weighted  average  protection  factor  for  n  exercises 
i  «  the  ith  exercise,  i  =  1,  2,  3,  ....  n 
wi  B  weighting  factor  for  the  ith  exercise 
PF  =  protection  factor  associated  with  a  particular  exercise 

Conventional  Method  of  Calculating  a  Protection  Factor 
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V'>\.  \W  ws* 


a.  normal  breathing  (NB) 

b.  deep  breathing  (DB) 

c.  turning  head  side-to-side  with  deep  breathing  (TH) 

d.  moving  head  up-and-down  with  deep  breathing  (UD) 

e.  talking  (T) 

f.  facial  grimacing  (FG) 

The  analysis  of  Figure  8  begins  at  the  bottom  of  the  strip-chart  record¬ 
ing.  The  first  section  of  information  uniquely  identifies  the  particular  sub¬ 
ject  and  type  of  respirator.  The  next  section  contains  the  instrument  cali¬ 
bration  data  (a  steady-state  response  measurement  is  recorded  for  each  of  the 
serial  dilutions  and  the  5X  or  10^  sodium  chloride  challenge  concentration). 
Before  the  breathing  exercises  begin,  a  zero  calibration  or  washout  measure¬ 
ment  is  taken  to  establish  a  baseline.  The  six  exercises  follow  in  sequence, 
each  being  performed  for  a  predetermined  time  period. 

The  cyclic  nature  of  the  recorder's  trace  during  the  exercises  is  a 
direct  function  of  the  subject's  breathing  cycle.  Figures,  for  instance, 
reveals  that  the  slight  negative  pressure  created  in  the  facepiece  during 
inhalation  increases  the  penetration  of  the  challenge  atmosphere.  Exhalation, 
on  the  other  hand,  creates  a  slightly  positive  pressure,  and  acts  to  reduce 
the  penetration  of  the  challenge  atmosphere.  Because  samples  are  drawn  from 
the  visual  cavity  of  chemical  defense  respirators,  absorption  of  the  sodium 
chloride  aerosol  by  the  lungs  is  negligible  [1].  Therefore,  respirator  per¬ 
formance  is  based  on  the  average  of  the  penetration  peaks  and  valleys  for  each 
of  the  exercises.  Finally,  the  overall  respirator  performance  is  based  on  the 
arithmetic  average  of  the  six  exercise  PF's. 

The  average  of  the  penetration  peaks  and  valleys  (location  of  their  mid¬ 
point)  is  generally  deduced  visually.  By  considering  each  exercise  separate¬ 
ly,  a  line  can  be  drawn  through  the  "visual  average"  of  the  peaks  and  val¬ 
leys.  This  midpoint  line  is  then  extended  until  it  intersects  the  calibration 
curve  (Fig.  8).  A  plot  of  the  sodiun  chloride  calibration  concentrations  vs, 
their  strip-chart  recorder  displacements  reveals  an  approximate  logarithmic 
relationship  (Fig.  9).  Thus,  logarithmic  interpolation  is  used  to  identify 
the  average  penetration  value  between  adjacent  calibration  decades.  Presented 
in  Table  3  is  a  sunmary  of  the  calculations  for  the  example  in  Figure  8. 


TABLE  3.  QUANTITATIVE  FIT  TEST  PENETRATION  RECORD 


£ 

$ 


Exercise 


Average  penetration 


Normal  breathing  5.5  x  10'^ 
Deep  breathing  ?./  x  10’j? 
Turning  head  side-to-side  (deep  breathing)  7.2  x  10’^ 
Moving  head  up-and-down  (deep  breathing)  5.5  x  10’^ 
Talking  5.9  x  10’£ 
Facial  grimacing  8.3  x  10"5 
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Figure  9.  Sodium  chloride  logarithmically  scaled  calibration  data  set  used  to  interpolate  mask  leakage 

concentration. 


With  this  information,  the  PF  for  each  exercise  can  he  calculated  using 
the  following  relationship: 


PF 


10°  1 
Cs  Cs 


(5) 


That  is,  since  the  challenge  atmosphere  is  generated  from  the  5%  or  10°  sodium 
chloride  solution  concentration,  and  10°  n  i,  it  follows  that  a  PF  is  simply 
the  reciprocal  of  the  exercise  average  penetration.  To  complete  this  example. 
Table  4  depicts  the  calculated  exercise  PF's  and  the  overall  arithmetic  aver¬ 
age  PF. 


TABLE  4.  QUANTITATIVE  FIT  TEST  PF  RECORD 


Exercise  '  PF 


Normal  breathing  1.8  x  104 
Deep  breathing  1.3  x  ID4 
Tuning  head  side-to-side  (deep  breathing)  1.4  x  104 
Moving  head  up-and-down  (deep  breathing)  1 .8  x  104 
Talking  •  1.7  x  1CT 
Facial  grimacing-  1 .2  x  104 


Overall  average  PF  =  l.F  x  104 


Although  the  strip-chart  recorder  data  can  be  interpreted  without 
significant  mathematical  rigor,  this  exercise  can  be  exasperating  when  more 
than  a  dozen  subjects  are  involved.  Having  analyzed  a  series  of  13  tests, 
this  author  has  developed  an  alternative  method  that  utilizes  the  IJSAFSAM 
PDP-11/70  computer  to  perform  these  ca (culations.  This  automated  scheme 
yields  a  data  reduction  turn-around  time  of  approximately  4  min  per  subject 
vs.  40  min  per  subject  by  manual  calculation. 


USING  A  VOLTAGE- TO-FREQUENCY  CONVERTER 
CIRCUIT  TO  DO  TIME-AVERAGED  INTEGRATION 

Various  techniques  have  been  applied  to  integrate  electronic  signals 
produced  in  the  laboratory.  These  techniques  include  the:  ball  and  disk 
mechanical;  low-inertia  motor;  electrochemical ;  analog-to-digital  conversion 
followed  by  counting;  operational  amplifier;  and  voltage- to- frequency  (V/F) 
conversion  followed  by  counting  [63].  The  first  four  techniques  have  been 
evaluated  by  numerous  Investigators  for  processing  laboratory-type  recorded 
signals.  Such  Integrators,  however,  possess  certain  inherent  disadvantages; 
e.g.,  high  cost,  Indirect  readout,  Insufficient  accuracy  and  precision,  and 
Insufficient  dynamic  range. 


Operational  Amplifier  Integrators 


Operational  amplifier  integrators  have  been  used  to  solve  the  short-term 
laboratory  signal  integration  problem  [44-50,53,56,68,62-64].  A  typical  oper¬ 
ational  amplifier  integrator  circuit  is  shown  in  Figure  10  [64].  The  output 
voltage  of  the  circuit  (Eout)  is  related  to  the  input  voltage  (Ein)  by  Equa¬ 
tion  6. 


-out 


-  -  -  / 
P.C  J0 


■in 


dt  + 


If  *b  dt 
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♦  1/ 
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ip  dt 
KC 


1_ 

RC 


t 

/  E  dt 
o  os 


(6) 


The  first  term  of  Equation  6  represents  the  desired  integrated  value;  and  the 
second  through  fourth  terms  represent  the  output  error  generated  by  inte¬ 
grating  the  input  bias  current  ( i 9) »  the  current  leakage  through  the  inte¬ 
grating  capacitor  (iRc),  and  the  offset  voltage  (Eos),  respectively.  These 

errors  can  be  minimized  by  restricting  the  period  of  integration  from  1  msec 
to  100  sec  [51,52,54,55,59,70],  However,  when  the  integrating  period 
approaches  100  sec  and  an  overall  IX  accuracy  is  desired,  a  large  value,  hi qh- 
performance,  expensive  polystyrene  capacitor  is  required,  as  well  as  an 
expensive  operational  amplifier  whose  input  offset  current  and  drift  are 
negligible.  In  order  to  optimize  accuracy,  to  integrate  sfqnals  lasting  from 
milliseconds  to  hours,  and  to  keep  the  cost  of  components  within  manageable 
limits,  the  V/F  integrator  scheme  becomes  an  extremely  attractive  alternative 
[69,71-73]. 


Voltage- to- Frequency  Converter  Integrators 

The  basic  function  of  a  V/F  converter  is  to  transform  a  variable  direct- 
current  voltage  (usually  0  to  10  volts)  into  a  pulse  train  whose  repetition 
rate  (frequency)  is  a  direct  linear  function  of  the  direct- current  input  volt¬ 
age,  An  excellent  technique  for  precisely  integrating  an  analog  voltage  sig¬ 
nal  is  Simply  to  add  a  counter  stage  to  the  output  of  a  V/F  converter  and 
accumulate  the  pulse  count.  By  accumulating  the  V/F  converter  output  pulses, 
the  "area  under  the  input  voltage  curve,"  or  integral,  is  calculated.  This 
operating  principle  [68,69]  is  illustrated  in  Figure  11  In  addition,  digital 
integrators  using  V/F  *  s  have  three  primary  advantages  over  their  operational 
amplifier  counterparts- -for  it  is  very  easy  to: 

a.  simply  switch  the  counter  fco  a  “hold"  position,  and 
the  accumulated  count  (integrated  value)  is  held  in¬ 
definitely  with  absolutely  no  drift; 

b.  preset  the  digital  counter  to  any  desired  level  and 
integrate  up  (or  down)  from  that  initial  condition; 
and 

c.  adapt  standard  electro-optical  readouts  (such  as  Nixie 
tubes,  light-emitting  diodes,  etc.)  to  display  the 
integrated  count  value. 

The  attractive  features  of  the  V/F  integration  technique  stimulated  the 
development  of  the  integrator  used  with  the  USAFSAH  RQFT  sodiuo  chloride 
instrument. 
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Operating  principle  of  a  V/F  integrator 


USAFSAM  Sodium  Chloride  RQFT  Voltage- to-Frequency 
Integrator  Circuit  Design 

Shown  in  the  respective  figures  are  the  schematic  diagrams  for:  the 
USAFSAM  sodium  chloride  RQFT  instrument  (Fig.  12);  the  digital  display 
(Fig.  13);  and  the  voltmeter  (Fig.  14).  The  major  components  [66,67,74-82] 
are  listed  in  Table  5. 


TABLE  5.  MAJOR  COMPONENTS  OF  THE  V/F  CONVERTER  ELECTRONIC  INTEGRATOR 


Schematic  Diagram  Ref¬ 
erences  (Figs.  12  -  14) 

Description 

R1 

27K  ohms,  1/4  W, 

5% 

R2 

IK  ohms,  1/4  W, 

5% 

R3 

20K  ohms,  1/4  W, 

5% 

R4 

500-ohm  variable 

potentiometer 

R5 

50K-ohm  variable 

potentiometer 

R6 

100K  ohms,  1/4  W, 

5% 

R7 

100K  ohms,  1/4  W, 

5X 

R8 

50K  ohms,  1/4  W, 

5% 

R9 

100K  ohms,  1/4  W, 

5X 

RIO 

50K  ohms,  1/4  W, 

5% 

Rll 

100K  ohms,  1/4  W, 

5X 

R12 

100K  ohms,  1/4  W, 

5% 

R13 

15K  ohms,  1/4  W, 

5X 

R14 

10K  ohms,  1/4  W, 

5% 

R15 

IK  ohms,  1/4  W, 

5% 

(Cont4d.  on  facing  page) 


TABLE  5  (Cont'd.) 


Schematic  diagram  ref- 
erences  (Figs.  12  -  14) 

R16 

R 1 7 

Cl 

T1 

IC1 

IC2 

IC3 


IC4 


ICS 


1C6 


IC7 


Digital  Display  Module  - 
Integrator  Count 

Digital  Display  Module  - 
Integration  Time 


Description 


lOOK-ohm  variable  potentiometer 

10K  ohms,  1/4  watt,  5% 

0.1-yF  capacitor 

2N956  NPN  transistor 

Precision  Monolithics  Incorporated, 
operational  amplifier,  OP-7 

Analog  Devices  Incorporated,  high  per¬ 
formance  V/F  converter,  AD450J 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  decade  counter/divider, 

MC14017B 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  noninverting  hex  buffers,  MC 
14050B,  Vcc  -  Pin  1,  Vss  -  Pin  8, 
Ground  Pins  11  and  14,  Vcc  =  +5  volts 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  quad  2- input  OR  gate,  MC14071B, 
Vcc  -  Pin  14,  Vss  -  Pin  7,  Ground  Pins 
5,  6,  /,  8,  9,  12,  and  13,  Vcc  °  +5 
volts 

Texas  Instruments  Incorporated,  hex 
inverter,  SN7404,  Vcc  -  Pin  14,  Vss  - 
Pin  7,  Vcc  ■  +5  volts 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  industrial  time  base  generator, 
MC14566B 

Dial ight  LED  display  module  749-1706 


Dialight  LED  display  module  749-1704 
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Figure  12:  Part  2  (of  3).  Sodium  chloride  RQFT  instrument  integrator. 


instrument  integrator. 


Figure  13.  Sodium  chloride  RQFT  instrument  integrator  digital  display. 

(♦Delete  these  two  digits' from  the  module  for  the  4-digit  LED  time  display.) 


Figure  14.  Sodium  chloride  RQFT  instrument  integrator  voltmeter. 


-  - 


Description  of  the  Analog  Devices  AD450J 
V/F  Converter  Integrated  Circuit 

The  heart  of  the  USAFSAM  sodium  chloride  RQFT  instrument  integrator  design 
is  the  Analog  Devices  AD450J  V/F  converter.  The  electronic  specifications  for 
this  device  [66,67,78]  are  summarized  in  Table  6,  and  a  block  diagram  of  the 
converter  is  shown  in  Figure  15  [66,67,78]. 

The  AD450J  is  a  low-cost  high-performance  V/F  converter  that  provides 
exceptional  linearity  and  temperature  stability  over  a  wide  dynamic  input  sig¬ 
nal  range.  The  key  to  the  AD450J's  precise  operation  is  the  unique  charge 
balance  conversion  technique.  This  feature  means  that  an  analog  signal  is 
accurately  converted  to  a  train  of  pulses  (each  of  which  has  a  constant  width 
and  amplitude)  at  a  rate  directly  proportional  to  the  analog  signal  ampli¬ 
tude.  The  output  continuously  tracks  and  responds  directly  to  changes  in  the 
input  signal. 

The  versatile  operational  amplifier  buffer  serves  as  the  input  stage.  Its 
purpose  is  to  convert  the  applied  input  voltage  to  a  control  current  for  the 
charge  balance  conversion  circuitry. 

External  clock  synchronization  is  not  required  for  the  AD450J.  The  inter¬ 
nal  temperature  compensated  timing  reference  yields  an  accurate  square-wave 
output  pulse  train  void  of  low-frequency  cycle-to-cycle  jitter. 

Finally,  the  output  drive  circuitry  is  responsible  for  conditioning  the 
stable  output  pulse  train.  The  high  current-handling  capacity  of  this  cir¬ 
cuitry  allows  the  designer  to  interface  directly  with  any  low-cost  digital 
processing  logic  family.  These  attractive  features  motivated  the  selection  of 
the  AD450J  as  the  V/F  converter  circuit  for  the  long-term  precision  integrator 
design. 


Operation  of  the  USAFSAM  Sodium  Chloride  RQFT 
V/F  Integrator 

Operation  of  the  USAFSAM  sodium  chloride  RQFT  V/F  integrator  is  quite  sim¬ 
ple.  The  integrator's  capability  and  function  can  be  appreciated  by  analyzing 
Figure  12.  Only  one  connection  to  the  integrator  is  required.  The  analog 
voltage  signal  to  be  integrated  (the  input  signal  to  the  strip-chart  recorder) 
is  connected  directly  to  the  0-10  volt  integrator  input  terminal. 

The  integrator  is  energized  with  the  on/off  switch,  and  the  front  panel 
neon  lamp  yields  a  positive  indication  for  applied  power.  Two  direct-current 
power  supplies  (±15  volt  and  +5  volt)  provide  the  operating  voltages  for  the 
integrator's  digital  circuits. 

Because  the  output  signal  of  the  sodium  chloride  RQFT  instrument  is  biased 
with  a  low-level  noise  component,  a  Precision  Monolithics  Corporation  OP-7 
operational  amplifier  buffer  is  useu  to  compensate  for  this  unwanted  signal 
[74].  The  degree  of  compensation  is  accomplished  by  rotating  the  PMT  noise 
offset  voltage  adjustment  potentiometer  (R16).  This  adjustment  must  be 
checked  each  time  the  RQFT  instrument  is  turned  on. 


TABLE  6.  ELECTRONIC  SPECIFICATIONS  OF  THE  ANALOG  DEVICES  AD450J  V/F  CONVERTER 


Character! Stic 

Value 

TRANSFER  FUNCTION 

Voltage  Input 

font  ■  (103  jp)ei„ 

ANALOG  INPUT 

Voltage  Signal  Range  (e-jn) 

0  to  +10V  min 

Overrange 

50%  min 

Impedance  (e-j n) 

20kfl 

Max  Safe  Input  Voltage  (e-jn) 

+25 V  (-Vs) 

ACCURACY 

Warmup  Time 

1  min 

Nonl Inearl ty 

ein  =  +lmV  to  +15V 

±0.01%  max 

Full  Scale  Error 

(+0.5  +  1.5)%  max 

Gain 

vs.  Temperature  (0  to  70°C) 

±50ppm/°C  max 

vs.  Supply  Voltage 

±200ppm/%max 

vs.  Time 

±100ppm/day 

Input  Offset  Voltage 

±5mV  max 

vs.  Temperature  (0  to  70°C) 

+50pV/°C 

vs.  Supply  Voltage 

±10ppm/%max 

vs.  Time 

±10uV/day 

(Cont'd.  on  facing  page) 


TABLE  6  (Cont'd.) 
Characteristic 


Value 


RESPONSE 

Settling  Time  for  +10V  Step  Input 
Overload  Recovery  Time 

OUTPUT 
Waveform 
Pulse  Width 
Rise/Fall  Time 
Pulse  Polarity 
Logic  "1"  (High)  Level 
Logic  "0"  (Low)  Level 
Capacitive  Loading 
Fan  Out  Loading 
Impedance 

POWER  SUPPLY 

Voltage,  Rated  Performance 
Voltage,  Operating 
Current,  Quiescent 

TEMPERATURE  RANGE 
Rated  Performance 
Operating 
Storage 

CASE  SIZE 


120ps 

15ms 

Train  of  TTL/DTL  compatible  pulses 
50us 
200ns 
positive 
+2.4V  min 
+0.4V  max 
lOOOpF  max 
10  TTL  loads  min 
3.3kn 

415V  dc 

4(12  to  18)V  dc 
(+15,  -9)mA 

0  to  +70*C 
-25°C  to  +80°C 
-55°C  to  +85°C 

(1,5  x  1.5  x  0.4)  inches 
transistor-transistor  logic. 


OTL  “  diode- transistor  logic;  and  TTL 
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g  Devices  AQ45QJ  V/F  converter 


The  overall  accuracy  and  dynamic  range  of  the  AD45Q0  V/F  converter  is 
established  with  an  initial  calibration  of  the  two  trim  potentiometers  (R4  and 
R5)  shown  in  Figure  12.  This  calibration  is  accomplished  after  the  integrator 
has  reached  an  operating  equilibrium  (5-min  warmup).  A  precision  regulated 
voltage  source  of  +1.00  mV  is  connected  to  pin  6  of'ICl,  and  kn  oscilloscope 
is  connected  between  ground  and  pin  6  of  IC2.  The  zero  offset  adjustment 
potentiometer  (R5)  is  rotated  so  that  a  1  Hz  frequency  pulse  train  is  visible 
on  the  oscilloscope.  The  precision  regulated  voltage  source  is  then  set  to 
+10.000  V,  and  the  full  scale  adjustment  potentiometer  (R4)  is  rotated  so  that 
a  10  kHz  frequency  pulse  train  is  visible  on  the  oscilloscope. 

The  analog  input  voltage-to-frequency  conversion  process  is  accomplished 
with  the  Analog  Devices  AD4500  V/F  converter.  Its  output  pulse  train  is  pro¬ 
cessed  by  the  Motorola  Semiconductor  Products  MC14017B  decade  counter-divider 
integrated  circuit.  This  integrated  circuit  conditions  and  accumulates  the 
V/F  converter  output  signal  for  the  6-digit  LED  integrator  count  display.  The 
user  can  select  a  signal  dividing  constant  via  the  panel -mounted  integrator 
count  scale  switch.  Selection  of  the  di vide-by-one  position  means  a  one-uo- 
one  correspondence  exists  between  the  actual  integrator  count  (area  under  the 
analog  voltage  input  signal)  and  the  magnitude  on  the  digital  display.  Selec¬ 
tion  of  any  switch  position  2  through  9  means  that  the  magnitude  on  the  digi¬ 
tal  display  must  be  multiplied  by  the  switch  position  number  to  yield  the 
actual  integrator  count.  This  feature  permits  the  user  to  integrate  input 
signals  for  several  hours  without  overloading  the  integrator  count  display. 

An  internal  time  base  generator  is  an  inherent  part  of  the  integrator 
design,  A  Motorola  Semiconductor  Products  MC14566  integrated  circuit  is  used 
to  provide  a  clocked  output  pulse  each  second.  The  accumulated  time  (in  sec¬ 
onds)  is  available  via  the  4-digit  LED  display.  The  user  controls  the  length 
of  time  an  analog  input  signal  is  integrated.  Two  front-panel  mounted 
switches  give  the  user  the  flexibility  to  clear  both  LED  displays  to  read  zero 
and  halt  the  integration  of  a  signal  at  any  point  in  time.  To  clear  the  dis¬ 
plays,  the  user  simply  toggles  the  display  reset  switch  and  trips  the  integra¬ 
tor  function  switch  to  the  hold  position.  When  the  operator  is  ready  to  inte¬ 
grate  a  signal,  he  toggles  the  integrator  function  switch  to  the  run  position. 
Thus,  the  signal  will  be  integrated  and  the  integrator  count  and  elapsed  time 
will  be  automatically  recorded  on  the  appropriate  LEO  displays.  When  the 
integration  is  to  be  halted,  the  use’  toggles  the  integration  function  switch 
to  the  hold  position.  This  last  operation  gives  the  user  the  opportunity  to 
record  the  magnitude  of  the  integrator  count  and  elapsed  time.  When  another 
signal  is  to  be  integrated,  the  user  resets  the  displays  to  zero  and  repeats 
the  sequence  described  above. 


Data  Collection  with  the  USAFSAM  Sodium 
4  Chloride  RQKT  Instrument  V/F  Integrator 

The  collection  of  RQFT  data  for  subsequent  PF  calculations  is  a  simple 
process  when  the  integrator  is  used.  Shown  in  Figures  16  and  17  are  the  data 
sheets  used  for  this  purpose.  After  the  RQFT  instrument  has  reached  its  oper¬ 
ating  equilibrium,  and  the  integrator's  PMT  noise  offset  voltage  has  bean 
adjusted,  the  sodium  chloride  calibration  and  challenge  atmosphere  atomizer 
solutions  are  processed. 
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RESPIRATOR  QUANTITATIVE  FIT  TESTING 
--USAFSAM  SALT  FOG  INSTRUMENTATION— 


Concentration 

10  to  the  Zero 
10  to  the  Minus  One 
10  to  the  Minus  Two 
10  to  the  Minus  Three 
10  to  the  Minus  Four 
10  to  the  Minus  Five 
10  to  the  Minus  Six 


SUBJECT  NAME: 
TYPE  OF  MASK: 
DATE  TESTED: 
TIME  TESTED: 


Exercise 


Normal  Breathing  Straight  Ahead 
Deep  Breathing  Straight  Ahead 
Talking 

Side-to-Side  Head  Movements  (Deep  Breathing) 
Up-and-Down  Head  Movements  (Deep  Breathing) 
Facial  Grimacing 


Count  Time  Period 

(In  Seconds) 


Figure  16.  Sodium  chloride  RQFT  data  collection  form  No.  1. 
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RESPIRATOR  QUANTITATIVE  FIT  TESTING 
-USAFSAM  SALT  FOG  INSTRUMENTATION— 


Concentration 

10  to  the  Zero 
10  to  the  Minus  One 
10  to  the  Minus  Two 
10  to  the  Minus  Three 
10  to  the  Minus  Four 
10  to  the  Minus  Five 
10  to  the  Minus  Six 


Voltage  (In  Volts) 


SUBJECT  NAME: 
TYPE  OF  MASK: 
DATE  TESTED: 
TIME  TESTED: 


Exercise 


Count  Time  Period 

(In'  "Seconds) 


Normal  Breathing  Straight  Ahead 

Normal  Breathing  Left 

Normal  Breathing  Right 

Normal  Breathing  Down 

Normal  Breathing  Up 

Deep  Breathing  Straight  Ahead 

Deep  Breathing  Left 

Deep  Breathing  Right 

Deep  Breathing  Down 

Deep  Breathing  Up 

Talking 

Fa.ial  Grimacing 

Side-to-Side  Head  Movements  (Normal  Breathing) 
Up-and-Qown  Movements  (Normal  Breathing) 
Side-to-Side  Head  Movements  (Deep  Breathing) 
Up-and-Down  Head  Movements  (Deep  Breathing) 


Figure  17,  Sodium  chloride  RQFT  data  collection  form  No.  2. 
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The  data  collected  for  each  of  the  standards  is  the  steady-state  output  voltage 
j  displayed  on  the  integrator's  digital  voltmeter.  These  voltages  will  range 

from  a  maximum  of  4.0  volts  to  a  minimum  of  0.09  volts;  each  reading  is  re- 
j  corded  opposite  the  sodium  chloride  concentration  (Fig.  16).  The  next  block  of 

|  information  to  be  recorded  is  the  subject  and  respirator  identification  infor¬ 

mation.  After  each  exercise  is  accomplished,  the  integrator  count  on  the 
\  6-digit  LED  display  and  time  on  the  4-digit  LED  display  are  recorded.  The 

i  elapsed-time  display  is  used  to  initiate  and  terminate  each  exercise.  Before 

proceeding  to  the  next  exercise,  the  integrator  displays  are  reset  to  zero. 

.  * 

After  all  subjects  have  been  tested,  the  user  proceeds  to  a  computer  term¬ 
inal  and  enters  the  information  from  the  data  collection  forms.  The  interac¬ 
tive  curve-fitting  and  PF  calculation  program  described  in  the  next  section  is 
used  to  process  the  RQFT  data. 


LEAST  SQUARES  CURVE  FITTING  COMPUTER  PROGRAM 
TO  CALCULATE  PROTECTION  FACTORS 

The  fitting  of  empirical  data  by  formulas  or  equations  can  be  accomplished 
by  two  methods.  The  first  is  to  have  a  polynomial  that  is  satisfied  exactly  at 
the  observed  data  points;  this  is  commonly  referred  to  as  "the  polynomial 
interpolation  method"  [83 ,84,86-92] .  The  second  method,  however,  is  a  more 
desirable  way  to  analyze  data  which  have  been  gathered  from  experimental  obser¬ 
vations  that  are  biased  with  various  errors  of  measurement;  this  is  commonly 
referred  to  as  "the  least  squares  approximation  method"  [83-87,93].  Because 
the  measurement  errors  associated  with  the  RQFT  calibration  data  have  been 
empirically  determined  to  be  relatively  small  in  magnitude,  a  least  squares 
curve  fitting  algorithm  is  used  to  determine  an  interpolating  polynomial.  The 
theory  of  fitting  nonlinear  curves  by  the  method  of  least  squares  yields  an 
interpolating  polynomial  when  the  degree  of  the  polynomial  is  one  less  than  the 
number  of  data  points.  Additionally,  this  curve  fitting  polynomial  is  unique 
and  the  data  points  are  fully  parameterized  [96].  A  summary  of  the  mathemati¬ 
cal  theory  for  the  method  of  least  squares  curve  fitting  is  presented  first, 
and  is  followed  by  its  direct  application  to  process  the  RQFT  data  and  calcu¬ 
late  PF's. 


The  Method  of  Least  Squares  Curve  Fitting 

The  objective  of  the  method  of  least  sguares  curve  fitting  is  to  relate  by 
some  function,  y  «  f(x),  a  set  of  m  points  (Xj,  yj ) ,  (j  *  1,  2,  3,  ...,  m), 
which  have  been  gathered  through  some  measuring  process.  The  method  of  least 
squares  curve  fitting  assumes  that  the  function,  y  »  f(x),  can  be  written  as  a 
polynomial  of  degree  n  <  m: 


y  «  a0  +  ajx  +  a2x2  +  ...  +  =  j  a.xi 

1=0  1 


(7) 
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The  next  step  in  the  process  is  to  determine  the  value  of  the  coefficients 
a-j ,  (i  =  0,  1,  2,  n) ,  such  that  the  polynomial  described  by  Equation  7 
is  a  "good  fit"  to  the  data  (xj,  yj).  By  substituting  the  data  points 
into  the  polynomial  described  by  Equation  7,  a  set  of  m  simultaneous  equations 
are  generated: 


=  a0  +  alxl  +  a2Xi2  +  ...  +  apX^  -  yt 

o  n 

R2  =  3q  ■**  X2  +  62*2  *  •••  *  3nX2  “  y2 


(8) 


Rm  =  ao  +  aixm  +  a2xm  +  ...  +  anxm  -  ym 


These  equations  are  not  exactly  equal  to  zero,  because  the  polynomial  does  not 
necessarily  pass  through  all  the  points  (except  in  the  case  where  the  degree 
of  the  polynomial  is  one  less  than  the  number  of  data  points).  As  shown  in 
Figure  18,  the  differences  between  a  polynomial  value  and  a  data  value  can  be 
positive,  negative,  or  zero.  This  difference  is  called  a  residual.  Residuals 
are  readily  calculated  by  Equation  9. 


Rj  3  l  «ixl  -  yj  for  (j  =  1,  2,  ....  m)  (9) 
1=0  J 

Thus,  the  set  of  equations  described  by  Equation  8  can  be  referred  to  as  "the 
residual  equations."  The  principle  of  least  squares  curve  fitting  states  that 
the  best  representation  of  the  data  is  that  which  makes  the  sum  of  the  squares 
of  the  residuals  a  minimum  [83,86,87,93].  Therefore,  it  is  desirable  to  force 
the  function 


f(a0,  ai,  a2,  ....  an)  »  Ri2  ♦  R22  +  R32  +  ...  +  Rm2  (10) 


to  be  as  close  to  zero  as  possible.  The  condition  which  fulfills  this 
requirement  is  that  the  partial  derivatives  of  Equation  10  be  exactly  zero 
[85,94]. 


“5*5 r 


Substituting  the  results  of  Equation  12  into  Equation  11  yields: 


Ri  +  R2  +  R3  +  ...  Rm  =  0 


Xi Ri  +  x2^2  +  X3R3  +  ...  +  XmRm  =  0 
X!2Rl  +  x22R2  +  x32R3  +  ...  +  xm2Rm  =  0 


(13) 


XinRi  +  x2nR2  +  x3nR3  +  ...  +  xmnRm  =  0 


Replacing  the  Rj's  by  their  values,  defined  in  Equation  8,  and  collecting 
the  coefficients  of  the  (n+1)  unknowns  of  a-,-,  (i  =0,  1,  2,  3,  ...»  n), 
yields: 


mag 

+  l*i*i 

-  I x j 2  a2 

+  ...  +  Ix-j% 

-  hi  *  0 

Ixja0 

+  Ixj2a1 

+  Ixj3a2 

+  ...  +  Ixjn+1an 

l^y j  *  0 

Ixj2a0 

+  l*i\ 

+  Ixj4a2 

+  ...  +  Ix^+2an 

-  Ixj2yj  =  0 

+  +  lxjr>+2a2  +  ...  +  Ixj2%  -  I x j ny j  ■  0 

where  all  sunwations  are  from  1  to  m;  that  is, 

IXj3  «  X!3  +  X23  +  X33  +  ...  +  Xjjj3 

Ixj2yj2  0  *iay2  +  x2zy2  ♦  ♦  xm2yft1 


(14) 


(15) 


The  set  of  equations  defined  in  Equation  14  are  known  as  the  normal  equa¬ 
tions.  All  of  the  associated  summations  are  known;  $0  the  system  of  equations 


defined  by  Equation  14  is  a  system  of  (n+1)  linear  equations  in  the  (n+1) 
unknowns  of  a-j  (i  =  0,  1,  2,  ...,  n).  The  solution  of  Equation  14  yields 
the  coefficients,  a-j ,  and  thus  the  polynomial  defined  by  Equation  7  is 
determined. 

The  principle  of  least  squares  curve  fitting  is  not  limited  to  polynomi¬ 
als  defined  solely  in  terms  of  x.  The  functional  relationship  between  y  and 
x  can  be  any  known  form  (i.e.,  ex,  log10x,  sin  x,  1/x,  etc.)  as  long  as  the 
functional  form  is  defined  and  the  resulting  normal  equations  can  be  solved. 

The  simplest  case  by  far,  however,  is  for  y  defined  directly  in  terms  of  x 
[83-87,93]. 

The  coefficients  for  Equation  7  are  found  by  calculating,  from  the  mea¬ 
sured  data,  all  the  sums  defined  by  Equation  14.  These  sums  are  then  substi¬ 
tuted  into  this  system  of  equations,  and  the  coefficients  aj  (i  =  0,  1,  2, 

....  n)  are  thus  determined. 

Application  of  the  Method  of  Least  Squares 
Curve  Fitting  to  Calculate  RQFT  PF's 

The  method  of  least  squares  curve  fitting  was  applied  to  process  the  cal¬ 
ibration  and  exercise  integrator  count  data;  the  final  product  is  the  calcula¬ 
tion  of  PF's.  The  author  wrote  a  Fortran  computer  program  that  accepted  the 
calibration  data  and  calculated  a  least  squares  polynomial  curve  fit  equa¬ 
tion.  The  integrator  count  data  for  each  exercise  was  then  substituted  into 
the  polynomial  equation,  and  a  corresponding  mask  leakage  penetration  (concen¬ 
tration)  was  calculated.  Finally,  PF's  were  calculated  using  Equation  5.  The 
subsequent  paragraphs  of  this  report  outline  the  analysis  used  to  accomplish 
this  method  of  calculating  PF's. 

The  objective  of  writing  a  least  squares  curve  fitting  computer  program 
was  to  take  advantage  of  the  intrinsic  importance  of  the  calibration  data  and 
relate  this  set  of  seven  points  (xj,  yj),  (j  =  1,  2,  3,  ....  7)  by  some 
function  y  =  f(x).  In  this  particular  case,  the  x  coordinate  is  the  digital 
voltmeter  reading  (in  volts),  and  the  y  coordinate  is  the  corresponding 
sodium  chloride  calibration  or  challenge  concentration.  (This  determination 
was  made  to  avoid  introducing  additional  error  in  the  calculation  of  a  PF  by 
having  to  implement  a  complex  inverse  iteration  algorithm.)  By  referring  to 
Figure  16,  the  reader  can  observe  that  these  seven  ordered  pairs  correspond  to 
the  data  at  the  top  of  the  RQFT  form.  A  typical  set  of  calibration  data  are 
*hown  in  Table  7. 

After  several  sets  of  calibration  data  had  been  examined,  and  various 
functional  definitions  of  the  variable  had  been  tried,  it  was  determined  that 
a  polynomial  of  the  following  form  yielded  the  best  and  most  stable  fit: 

y  =  a0  +  ajex  +  a2(ex)2  +  a3(ex)3  +  a4(ex)H  +  a5(ex}5  +  a6(ex)6  (16) 

The  "goodness"  of  fit  was  checked  utilizing  the  conditional  Eq.  10  to  verify 
that  the  polynomial  indeed  passed  through  the  RQFT  calibration  data  points. 

The  computer  program  has  calculated  the  "goodness"  of  fit;  and  the  results  for 
a  typical  run,  using  the  data  In  Table  7,  are  presented  in  Table  8.  In 
addition,  a  computer- generated  (Calcomp)  plot  (semilogarithmic)  is  shown  In 
Figure  19. 
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TABLE  7.  TYPICAL  SODIUM  CHLORIDE  RQFT  CALIBRATION  DATA  USED 
FOR  LEAST  SQUARES  CURVE  FITTING 


X-coordinate 
(voltage  reading) 

[in  volts] 

Y-coordinate 
(sodium  chloride 
concentration) 

3.380 

1.0 

2.915 

0.1 

2.310 

0.01 

1.500 

0.001 

0.545 

0.0001 

0.165 

0.00001 

0.105 

0.000001 

TABLE  o.  LEAST  SQUARES  CURVE  FIT  CALCULATIONS 
FOR  THE  DATA  CONTAINED  IN  TABLE  7 


The  order  of  the  desired  polynomial  =  6 

The  polynomial  functional  definition  of  the  variable  (X)  is  in  terms  of: 
(Exponential  (X)) 


Coefficient 

Number 

1 

= 

-4.17904E-05 

Coefficient 

Number 

2 

S 

-7.31487E-05 

Coefficient 

Number 

3 

as 

1.25863E-04 

Coefficient 

Number 

4 

= 

-2.68095E-05 

Coefficient 

Number 

5 

= 

3.80138E-06 

Coefficient 

Number 

6 

s 

-1.71541E-07 

Coefficient  Number 

7 

= 

3.88406E-09 

i  residuals  are 

calculated  by  the  following  equation: 

[Y (DATA )- Y(ESTIMATED)] 

SAMPLE  NUMBER 

Y(DATA) 

Y(ESTIMATED) 

RESIDUAL 

1 

1.00000 

1.00000 

-0.483548E-07 

2 

0.100000 

0.100000 

-0.189111E-07 

3 

0.100000E-01 

0.100000E-01 

-0.711342E-08 

4 

0.100000E-02 

0.100000E-02 

-0.401929E-08 

5 

Q.100000E-03 

0.100004E-03 

-0.390294E-08 

6 

0.100000E-04 

0.100034E-04 

-0.340602E-08 

7 

0.100000E-05 

0.100433E-05 

-0.432502E-08 

Sum  of  the  square  residuals  =  2.80811E-15 
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From  the  information  in  Table  8,  the  polynomial  curve  fit  equation  is: 

y  =  (-4.17904  x  10'5  )  +  (-7.31487  x  10"5)ex  +  (1.35863  x  lO’4^)2  + 
(-2.68095  x  10-5)(ex)3  +  (3.80138  x  10‘6)(ex)4  + 

(-1.71541  x  10-7)(ex)5  +  (3.88406  x  10"9)(ex)6  (17) 


The  residuals  were  calculated  to  measure  the  "goodness  of  fit"  for  the  poly¬ 
nomial  curve  fit  Equation  17.  Since  the  RQFT  data  points  are  fully  parameter¬ 
ized  by  Equation  17,  the  residuals  and  sum  of  the  square  residuals  should  be 
exactly  zero.  However,  these  values  are  not  exactly  equal  to  zero  because  of 
computer  calculation  round-off  errors.  The  computer  program  made  these  calcu¬ 
lations  by  substituting  the  set  of  sodium  chloride  calibration/challenge  con¬ 
centration  (Y-coordinate  values  of  Table  7)  values  into  the  polynomial  curve 
fit  Equation  17;  the  corresponding  calculated  Y  values  are  called  the  Y-esti- 
mated  values.  One  can  readily  observe  that  the  sum  of  the  square  residuals  is 
2.80811  x  10“ 1 5 ,  a  very  small  number  indeed.  Interestingly,  the  worst-case 
residual  occurs  for  the  10'6  sodium  chloride  calibration  standard;  this 
finding  implies  that  the  accuracy  of  predicting  a  particular  exercise  mask 
leakage  concentration  will  not  be  in  error  by  more  than  0.433%  [(0.100433  x 
lO'5  -  0.100000  x  10-5)/(0. 100000  x  10-5)  =  0.00433  x  100].  The  magnitude  of 
this  error  is  representative  of  that  found  for  several  sets  of  calibration 
data  that  have  been  analyzed. 

The  following  explanation  is  the  key  to  understanding  how  a  best  fit 
polynomial  curve  (relating  sodium  chloride  calibration/challenge  concentration 
to  a  digital  voltmeter  response)  can  be  used  to  calculate  a  PF.  Two  concepts 
are  involved  in  the  analysis. 

First:  the  Y-axis  is,  in  reality,  an  exact  scale  for  the  time-averaged 
mask  leakage  penetration  (concentration).  When  the  user  calibrates  the  RQFT 
instrument,  the  strip-chart  recorder  displays  (on  its  vertical  or  Y-axis)  the 
PMT's  response  for  a  known  diluted  sodium  chloride  concentration  sample.  It 
is  important  to  recognize  that  the  response  for  the  diluted  concentration  is 
equivalent  to  that  which  would  be  measured  if  one  were  to  evaluate  the  fit  of 
a  mask  with  the  10°  sodium  chloride  challenge  concentration,  and  have  the  mask 
fit  be  perfect,  except  for  a  "calibrated"  leak  that  would  permit  a  correspond¬ 
ing  known  diluted  concentration  of  the  10°  challenge  to  penetrate  the  mask. 

The  essence  of  the  polynomial  curve  fit  equation  is  that  it  permits  "interpo¬ 
lation"  between  the  responses  for  adjacent  sodium  chloride  calibration  con¬ 
centrations. 

Second:  the  integrator  count--the  source  of  data  used  to  calculate  a 
mask  leakage  penetration  (concentration)--is,  in  reality,  a  time-averaged 
voltage  response.  This  fact  can  be  derived  through  the  following  analyses: 

a.  Integrator  sensitivity  is  1000  counts-per-volt*sec. 

b.  Each  exercise  is  performed  for  a  predetermined  length  of  time;  for 
example,  10  sec. 

c.  The  Integrator  count  (IC)  value  recorded  for  a  particular  exercise 
Is  actually  the  time-averaged  area  under  the  strip-chart  recorder 
response  (refer  to  Figs.  8  and  11). 


Therefore, 


in  which  V  volts  is  the  time-averaged  voltage  for  a  particular  exercise. 

The  computer  program,  written  to  accept  the  integrator  exercise  count 
data  and  the  time  duration  for  each  exercise,  calculates  corresponding  exer¬ 
cise  time-averaged  voltages.  Fach  time-averaged  voltage  is  substituted  into 
the  polynomial  curve  fit  Equation  17,  and  a  corresponding  mask  leakage  pene¬ 
tration  (concentration)  is  calculated.  The  final  set  of  calculations  per¬ 
formed  by  the  program  are  the  individual  exercise  PF‘s,  an  arithmetic  average 
PF,  and  a  time-weighted  average  PF.  Equation  5  is  used  to  calculate  individ¬ 
ual  PF's,  and  Equations  3  and  4  are  used  to  calculate  the  average  PF's.  The 
results  of  the  data  in  Tables  7  and  8  are  given  in  Table  9 

Discussion  of  the  Computer  Programs  Used  to  Process  the 
USAFSAM  Sodium  Chloride  RQFT  Integrator  Data 


Two  computer  programs  are  used  to  process  the  integrator  exercise  data. 
The  first  and  primary  program  is  called  NACIRQFT.FTN;  the  second, 
NACIGRAPH.FTN. 


The  purpose  of  NACIRQFT.FTN  is  to  use  the  RQFT  information  collected  on 
the  RQFT  data  sheet  (Fig.  16)  and  calculate  a  set  of  PF's.  The  results  of 
this  program  are  stored  on  three  disk  files: 

1.  OATA.XXX  contains  the  initial  calibration  data,  the 
test  identification  data,  the  time  period  for  each 
exercise,  and  a  listing  of  the  exercises  performed  and 
their  associated  integrator  count  values. 


CALCX.XXX  contains  the  polynomial  curve  fit  equation 
coefficients,  the  calculated  residuals,  the  sum  of  the 
square  residuals,  a  composite  listing  of  the  identifi¬ 
cation  data,  the  exercises  performed  and  their  corre¬ 
sponding  PF's,  and  the  average  PF's. 

GRPHX.XXX  contains  an  array  of  401  X-axis  values  and 
401  Y-axis  values.  These  values  were  generated  using 
the  polynomial  curve  fit  Equation  15.  This  array  Is 
used  by  the  NACLGRAPH.FTN  program  to  generate  a 
CALCOHP  calibration  curve  for  the  soditia  chloride  RQFT 
cal i brat ion  data. 
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TABLE  9.  PROTECTION  FACTOR  CALCULATIONS  FOR  THE  DATA  CONTAINED  IN 
TABLE  7 


SUBJECT  NAME:  Captain  Edward  S.  Kclesar,  Jr. 
TYPE  OF  MASK:  USA:  H17  -  Medium  (no  glasses) 
DATE  TESTED:  9  April  1980 
TIME  TESTED:  1330  hours 

EXERCISE  INTEGRATOR  COUNT  DATA 


Exercise 

Integrator 

count 

Time  period 

~ T 

Protection  factor 

Normal  breathing  straight  ahead 

3904 

10 

1 .8E+04 

Deep  breathing  straight  ahead 

4751 

10 

1.3E+04 

Talking 

4628 

10 

1 ,3E+04 

Side-to-side  head  movements 
(deep  breathing) 

3976 

p 

1./T+04 

Up-and-down  head  movements 
(deep  breathing) 

40  i' 6 

m 

1 .7E+04 

Facial  grimacing 

4937 

10 

1.2E+04 

Overall  arithmetic  average  protection  factor  (PF)  for  all  categories  of 
exercises  actually  performed  *  1.5E+04. 

Overall  time-weighted  average  protection  factor  (PF)  for  all  categories  of 
exercises  actually  performed  •  1.5E+04. 


The  NACLGRAPH.FTN  program  has  been  written  as  a  separate  program  because 
the  USAFSAM  POP- 11/70  CALCOMP  plotter  is  an  off-line  device.  That  is, 
NACLGRAPH.FTN  uses  the  GRPHX.XXX  disk  file  as  input  data  and  sequences  through 
the  CALCOMP  plot  subroutine  library  files,  and  produces  a  user-named  output 
disk  file  (for  example,  PLOT.SCD}.  The  NACLGRAPH.f  "N-generated  output  file  is 
then  transferred  to  a  magnetic  tape  by  one  of  the  computer-room  operators. 

The  magnetic  tape  is  then  mounted  on  the  CALCOMP  terminal  and  the  plot  is 
generated  [95]. 

Each  of  the  programs  discussed  in  this  section  is  documented  with  com¬ 
ments  that  define  the  variables  and  explain  the  operations  performed.  There¬ 
fore,  a  line-by-line  analysis  of  the  code  will  not  be  done.  For  the  inter¬ 
ested  reader,  however,  the  following  information  is  available  in  eight  appen¬ 
dixes  (A  -  H) 

Appendix  A:  NACLRQFT.FTN  Fortran  listing 

Appendix  B:  DATA.XXXX  file  contents  for  data  in  Table  7 

Appendix  C:  CALCX.XXX  file  contents  for  information 
in  Tables  8  and  9 

Appendix  0:  GRPHX.XXX  file  contents  for  use  with 
NACLGRAPH.FTN  program 

Appendix  E:  NACLGRAPH.FTN  fortran  listing 

Appendix  F:  CALCOMP  generated  plot  (semi logarithmic)  of  the 
GRPHX.XXX  data. 

Appendix  G:  User's  guide  for  the  NACLRQFT.FTN  computer  program 

Appendix  H:  User's  guide  for  the  NACLGRAPH.FTN  computer  program 
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c 

C  THIS  PROGRAM  CALCULATES  PROTECTION  FACTORS  FOR  DATA  COLLECTED  ON 
C  THE  USAFSAM/VNL  SODIUM  CHLORIDE  RQFT  INSTRUMENT. 

C 

C  CALCULATIONS  ARE  BASED  ON  INTEGRATING  THE  PHOTOMULTIPLIER 
C  TUBE'S  OUTPUT  OVER  A  VARIABLE  TIME  PERIOD  (10-20  SECONDS)  FOR 

C  EACH  EXERCISE  OF  THE  RQFT  PROTOCOL.  SINCE  THE  INTEGRATOR'S  j 

C  SENSITIVITY  IS  1000  COUNTS  PER  (VOLT-SECOND),  AND  THE  TIME  PERIOD 
C  FOR  EACH  EXERCISE  IS  VARIABLE  (10-20  SECONDS),  A  SCALED  OR 

C  'AVERAGE'  VOLTAGE  FOR  THE  PHOTOMULTIPLIER  OUTPUT  CAN  BE 

C  CALCULATED  BY  THE  FOLLOWING  RELATIONSHIP: 

C  I 

C  APTVO=  (IC)  *  (1/IS)  *  (1/VTPEE) 

C  i 

C  WHERE,  APTVO=AVFRAGE  PHOTOMULTIPLIER  TUBE  VOLTAGE  OUTPUT  j 

C 

C  IC= INTEGRATOR  COUNT 

C 

C  IS = INTEGRATOR  SENSITIVITY 

C 

C  VTPEE= VARIABLE  TIME  PERIOD  FOR  EACH  EXERCISE 

C 

C  NEXT,  A  LEAST-SQUARES  CURVE  FITTING  METHOD  IS  USED  TO  GENERATE  A 
C  POLYNOMIAL  FUNCTION  OF  THE  FORM:  Y-F(X).  THE  DATA  USED  TO 
C  GENERATE  THIS  FUNCTION  ARE  THE  ORDERED  PAIRS  OF  VALUES  OF  THE  FORM 
C  (X , Y) ,  TABULATED  WHEN  CALIBRATING  THE  RQFT  INSTRUMENT  WITH  THE 

C  SODIUM  CHLORIDE  CALIBRATION  STANDARDS.  IN  THIS  CASE; 

C 

C  X= PHOTOMULT I PL I ER  DETECTOR  TUBE  VOLTAGE  OUTPUT  (IN  VOLTS) 

C 

C  YaSODIUM  CHLORIDE  CAI IBRATION  STANDARD  CONCENTRATION 

C 

C  THE  RESULT  OP  THE  CURVE  FITTING  PROCEDURE  WILL  BE  THE  CAPABILITY 

C  TO  DIRECTLY  CALCULATE  THE  MASK  LEAK  CONCENTRATION  FOR  AN  EXERCISE. 

C  THIS  IS  SO  BECAUSE  THE  POLYNOMIAL  CURVE  FITTING  RELATIONSHIP  OF 
C  THE  FORM:  Y«F(X),  DIRECTLY  RELATES  THE  PHOTOMULTIPLIER  TUBE'S 
C  AVERAGE  OUTPUT  VOLTAGE  TO  CONCENTRATION  (THAT  IS,  THE  INTEGRATED 
C  COUNT  IS  DIRECTLY  RELATED  TO  MASK  LEAKAGE  CONCENTRATION  FOR  A 

C  PARTICULAR  EXERCISE)  .  j 

C  i 

C  THUS,  FOR  A  PARTICULAR  EXERCISE,  A  MASK  LEAKAGE  CONCENTRATION 

C  IS  CALCULATED  BY  SUBSTITUTING  THE  'ADJUSTED'  INTEGRATOR  COUNT,  .J 

C  (THAT  IS,  THE  APTVO  VALUE),  INTO  THE  POLYNOMIAL  CURVE  FITTING  i 

C  FUNCTION.  ! 

c  i 

C  FINALLY,  A  MASK  LEAK  CONCENTRATION  IS  CONVERTED  TO  A  PROTECTION 
C  FACTOR  BY  THE  FOLLOWING  RELATIONSHIP: 

C  1 

C  PF"  (CC)  /  (ML) 

C 

C  WHERE,  PF-PROTECTION  FACTOR  j 

C  j 

C  CC-CHALLENGE  CONCENTRATION  (FOR  THIS  SYSTEM,  10  TO  THE  | 

C  ZERO  OR  1.0)  { 

C  | 

C  ML-MASK  LEAK  CONCENTRATION  FOR  A  PARTICULAR  EXERCISE  j 

C 


i 

j 
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C  SINCE  THE  LEAST-SQUARES  CURVE  FITTING  PROCEDURE  IS  THE  KEY  TO 

C  THIS  METHOD  OF  PROCESSING  THE  RQFT  DATA,  THE  FOLLOWING  COMMENTS 

C  WILL  DESCRIBE  ITS  IMPLEMENTATION  IN  THE  PROGRAM. 

C 

C  THIS  IS  A  PROGRAM  WHICH  UTILIZES  A  LEAST-SQUARES  METHOD  OF 
C  APPROXIMATION  TO  DO  CURVE  FITTING. 

C 

C  THIS  PROGRAM  WILL  TAKE  A  SET  OF  (N)  NUMBER  OF  DATA  POINTS  AND 

C  WILL  FIT  A  POLYNOMIAL  (UP  TO  DEGREE  <N-1) 3  TO  THE  SET  OF  POINTS. 

C 

C  THIS  PROGRAM  DETERMINES  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  BY 
C  SOLVING  A  SYSTEM  OF  NORMAL  EQUATIONS  WHICH  ARE  DERIVED  BY 

C  CONSIDERING  SOME  QUANTITY  (Y)  AS  A  POLYNOMIAL  FUNCTION  OF  ANOTHER 

C  QUANTITY  (X)  . 

C 

C  TO  FIND  THE  BEST  FITTING  CURVE  FOR  A  GIVEN  SET  OF  DATA,  THE 

C  PROGRAM  USER  IS  ONLY  REQUIRED  TO  ENTER  THE  OBSERVED  DATA  POINTS 

C  AND  THE  DEGREE  OF  THE  POLYNOMIAL  APPROXIMATION. 

C 

C  THE  FIRST  PART  OF  THE  PROGRAM  (STATEMENT  LABEL  NUMBERS:  5101  THRU 
C  403)  ACCEPTS  THE  DATA,  STORES  IT,  AND  ORGANIZES  IT  FOR  USE  IN 
C  THE  REMAINDER  OF  THE  PROGRAM. 

C 

C  THE  SECOND  PART  OF  THE  PROGRAM  (STATEMENT  LABEL  NUMBERS:  400  THRU 

C  500)  ESTABLISH  THE  NORMAL  EQUATIONS.  STATEMENT  LABEL  NUMBERS:  30 

C  THRU  90  SOLVES  THE  NORMAL  EQUATION  SYSTEM.  IN  THE  FINAL  SECTION 

C  OF  THE  PROGRAM  (STATEMENT  NUMBERS:  16  THRU  11)  ,  THE  ERRORS 

C  (DIFFERENCES  BETWEEN  THE  ACTUAL  AND  THE  ESTIMATED  VALUES)  ARE 
C  CALCULATED.  IN  ADDITION  THE  SUM  OF  THE  SQUARE  ERRORS  IS 
C  CALCULATED.  THE  POLYNOMIAL  FUNCTION  PRODUCING  THE  SMALLEST 

C  SUM  OF  SQUARE  ERRORS  YIELDS  THE  BEST  CURVE  FIT.  OBSERVATION 
C  OF  THE  INDIVIDUAL  ERRORS,  ON  THE  OTHER  HAND,  GIVES  A  CLUE 

C  AS  TO  THEIR  RELATIVE  DISTRIBUTION  ABOUT  THE  BEST  FIT 

C  POLYNOMIAL  CURVE. 

C 

C  THE  USER  IS  GIVEN  THE  OPTION  OF  SELECTING  A  FUNCTIONAL 
C  DEFINITION  OF  THE  VARIABLE  (X)  SO  THAT  APPROXIMATIONS 
C  CAN  BE  MADE  TO  CURVES  THAT  ARE  NOT  NECESSARILY 

C  POLYNOMIALS  IN  (X).  FOR  EXAMPLE:  Y-EXP(X)  OR  Y-l/ (X)  . 

C 

C* ******************************************************************** 


c*****  ***** 

c***»*  IF  YOU  HAVE  ANY  QUESTIONS  CONCERNING  THIS  PROGRAN  CALL  ***** 

C*  *  *  *  *  ***** 

C*****  CAPTAIN  EDWARD  S.  KOLESAR,  JR.  ***** 

C*****  ***** 

C*****  USAFSAM/VNL  BROOKS  AFB,  TX  ***** 

C***»*  ***** 

*****  AUTOVON  240-2154  COMMERCIAL  (512)536-2154  ***** 

*****  ***** 


ARRAYS  ARE  DOUBLE  PRECISION  TO  INSURE  ACCURACY  OF  THE 
CALCULATIONS. 

IMPLICIT  INTEGER*4  (I-N) 
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INTEGER  FILTK 

DOUBLE  PRECISION  D< 8, 9 )  ,W(  0  : 13 )  , 2 ( 0 : 13 ) ,B ( 8, 9) ,E ( 8, 9) 

DOUBLE  PRECISION  XAS<7) ,A(14) ,X(13) ,Y(13) .THEORY , RESID,SRESD, XXX 
DOUBLE  PRECISION  DEL.XMIN ,XMAX ,XEST , BEGIN, XXMIN,XXMAS,XXAS(7 ) 
DOUBLE  PRECISION  PF (17 ) , IC ( 17 > , ICTP<17 ) ,XXXX 
DIMENSION  GRPH1 (14) ,GRPH2(802) ,ICCBS(2) 

DIMENSION  XX(401) ,YX(401) 

BYTE  SECN(9) 

BYTE  P ( 9)  ,C ( 1 0)  ,G (1 0) 

BYTE  ANS,COMMEN ,REP , COM 
BYTE  SELECT (7) .GROUP  1, GROUP  2 
BYTE  YES, NO 

BYTE  NAME (45) ,MASK(45) ,DATE(45) ,TIME(45) 

DATA  SECN/ ' 1 ' , 1 2 ' , ' 3 ' , 1 4 ' , ' 5 ' , ' 6  *  # ‘7 ' » 1 8 1 , * 9 ' / 

YES” 1 Y  * 

NO” '  N 1 


ESTABLISH  A  FILE  COUNTER  AND  DECURE  THE  FILE  NAMES. 
S101  FILTK*1 

VARIABLES  ARE  SET  EQUAL  TO  ZERO  SO  THAT  ITERATIVE  RONS 
ON  THE  SAME  DATA  CAN  BE  READILY  PERFORMED. 


H-0 

NP-0 

ICCBS II ) ”72*256+27 
ICCBS  (2)  -74 *2 56 *27 
4137  FORHATdH  ,2A2) 
IC1-0.0 
IC2-0.0 
IC3-0.0 
IC4-0.0 
ICS-0.0 
IC6-0.O 
IC7-0.0 
IC8-0.0 
IC9-0.0 
ICIO-O.^S) 

ICU-0.0 
IC1 2-0.0 
IC13-0.0 

104- 0,0 

105- 0,0 
108-0.0 
XO-O.O 
XC2-0.0 
XC3-0.0 
XC4-0.0 
XCS-0.0 
XC6-C.0 
XC7-0.0 
XCB-0.0 
XC9-0.0 

xo  o-o.o 

XC 11-0.0 
XO  2-0.0 
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XC.I 3=0 . 0 
XC14-0.0 
XC15-0. 0 
XC16-0 . 0 
CALCPF-0 . 0 
XEST-0 . 0 
PFEST-0.0 
DEL=Q . 0 
THEORY-0. 0 
RESID-0.0 
SRESD-0.0 
XXX-0.0 
BEGIN-0.0 
DO  9333  1-1,17 
IC ( I ) -0 . 0 
PF  ( I )  -  0 . 0 
ICTP ( I ) -0 . 0 

9333  CONTINUE 

DO  4153  1-1,7 
XAS ( I ) -0 . 0 
XXAS ' I ) -0 . 0 

4153  CONTINUE 

THE  FILES  THAT  HOLD  VARIOUS  SEGMENTS  OF  DATA  ARE  NAMED. 

DO  9334  1-1,13 
XtII-0.0 
YU1-0.Q 

9334  CONTINUE 
P( 1 ) • *  D’ 

PC2) -• A1 
P  C  3 )  » •  T  * 

Pt4)-*A‘ 

PtS) »•  .  * 

P 1 9  > -0 

etn-'c 

etzi-’A1 
C (3 !  ■ '  L 1 
Ct4)-'C 

051 -8KCN  <  FILTK ) 

06>«*  .  * 

0101-0 
GUI-'C 
OS)«'R' 

GtJi-'S’* 

Gt4!-*(l‘ 

OS)  -SECN  I F I LTK ) 

OS)«* 

0101-0 

THE  FILE  CALLED  DATA. XXX,  CONTAINS  THE  SODIUM  CHLORIDE 
CAUHRAT10N  STANDARD  CONCENTRATION  DATA  TO  HE  PITTED 
WITH  THE  POLYNOMIAL  FUNCTION. 

THE  FILE  CALLED  CALCX.XXX,  CONTAINS  THE  CALCULATED  POLYNOMIAL 
COEFFICIENTS,  RESIDUALS,  AND  OTHER  DESCRIPTIVE  INFORMATION. 
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THE  FILE  CALLED  GRPHX.XXX  CONTAINS  THE  FOLLOWING  ARRAYS: 

[GRPH1 1  SODIUM  CHLORIDE  CALIBRATION  5TANDARD 
CONCENTRATION  DATA. 

IGRPH2)  ESTIMATED  VALUES  DERIVED  FROM  EXERCISING 
THE  POLYNOMIAL  FUNCTION. 

IF  A  PLOT  OF  THE  ACTUAL  DATA  AND  ESTIMATE*  VALUES  IS  DESIRED, 
THE  USER  CAN  USE  THE  FILE  CALLED  GRPHX.XXX  FOR  THIS  PURPOSE. 


NUMBER  THE  FILES  SEQUENTIALLY  SO  THAT  THEY  CAN  BE  EASILY  RETRIEVED 
FOP  PRINTING  AND  PLOTTING. 

TYPE  2006 

2006  FORMATUX,’  '///> 

1002  FORMAT (1  All 

PROGRAM  USER  WILL  UNIQUELY  IDENTIFY  EACH  DATA  FILE  WITH  A 
SEQUENTIAL  NUMBERING  SYSTEM  FOR  EASE  OF  RECALL. 

TYPE  4137, ICCBS ( 1 1 , ICCBS ( 2) 

TYPE  2006 
TYPE  2002 

2002  FORMATUX. ’USER  ATTENTION:  IN  ORDER  TO  KEEP  TRACK  OF  THE  DATA’/ 

C'  SETS  BEING  ANALYZED,  IT  IS  RECOMMENDED  THAT  THEY  BE  '/ 

C‘  SEQUENTIALLY  NUMBERED.*/) 

TYPF  306 
TYPE  2003 

2003  FORMATUX, ’ENTER  THE  FOLLOWING:  QQ1  FOR  THE  FIRST  DATA  SETj  */ 

C‘  002  FOR  THE  SECOND  DATA  SET:  003  FOR  THE  THIRD  DATA  SET.  ETC.*/) 
TYPE  306 
TYPE  2004 

2004  FORMATUX,  'ENTRY*  *,S) 

ACCEPT  2005, P(6) ,P(7) ,P)8) 

2005  FORNATUA1) 

TYPE  2006 

TYPE  4137, teCHS U) , ICCBS  12) 

TYPE  2006 
TYPE  5100 

5100  FORMAT! IX,  MISER  ATTENTION:  IN  ORDER  TO  KEEP  TRACK  OF  THE*/ 

C  POLYNOMIAL  CURVE  FITTING  COEFFICIENTS  AND  X ESI DUALS,  IT*/ 

C*  IS  RECOMMENDED  THAT  THEY  BE  SEQUENTIALLY  NUMBERED.*/) 

TV  PE  306 
TYPE  5200 

5200  FORMAT* IX. 'ENTER  THE  FOLLOWING:  001  FOR  THE  FIRST  RESIDUAL  SET: */ 
G*  002  FOR  THE  SECOND  RESIDUAL  SET:  003  FOR  THE  THIRD,  ETC.*/) 

TYPE  306 
TYPE  5300 

3300  FORMATUX, ’ENTRY-  '  ,  f) 

ACCEPT  5400,0?)  ,C(8)  ,CC9> 

5400  FORMAT)  3AD 
TYPE  2006 

TYPE  41J3,iccasU)  ,ICCOS<2> 

TYPE  2006 
TYPE  4001 

4001  FORMATUX, 'USER  ATTENTION:  IN  ORDER  TO  KEEP  TRACK  OF  THE  GRAPH'/ 
C*  SETS  UBINC  ANALYZED,  It  IS  RECOMMENDED  TUAT  THEY  Ofi'/ 
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C'  SEQUENTIALLY  NUMBERED1/) 

TYPE  306 
TYPE  4004 

1004  FORMATUX,  'ENTER  THE  FOLLOWING:  001  FOR  THE  FIRST  GRAPH  SET;'/ 
C'  002  FOR  THE  SECOND  GRAPH  SET;  003  FOR  THE  THIRD,  ETC.1/) 

TYPE  306 
TYPE  2004 

ACCEPT  2005, G( 7) ,G(8) ,G!9> 

TYPE  2006 

TYPE  41  37  ,  ICCBS  (1)  ,ICCBS(2) 

TYPE  2006 
REP-'N1 
GO  TO  3081 
6000  FILTK=FILTK+1 

CIS) =SECN(FILTK> 

GC  5) »SECN I  FI LTK) 

CLOSE  (UNIT*  3) 

CLOSE! UNIT* 1) 

OPEN (UNIT* 3, NANE-C, DISPOSE-1 SAVE' .TYPE* ‘NEW1) 

OPEN  ( UNI  T*  1 ,  NANE  »G ,  DI S  PO  SE- 1  SAVE 1  ,TY  PE*  1  NEW 1 ) 

REWIND  2 
GO  TQ  6001 

3081  OPEN (UNIT*2,NAME*P, DISPOSE* ‘SAVE1 .TYPE* ‘NEW1) 

OPEN ( UNIT*  3 , NAME«C , DI SPOSE- 1  SAVE 1 .TYPE- 1  NEW 1 ) 

OPEN ( UN  IT- 1 .  NAHE-G, DISPOSE- 1  SAVE 1 ,TY  PE- 1  NEW 1 ) 

C 

C  ENTER  THE  DATA  TO  BE  PROCESSED.  IN  THIS  SITUATION.  THE  VOLTAGE 

C  MEASUREMENTS  ARE  THE  X-AXIS  DATA  AND  THE  SODIUM  CHLORIDE 

C  STANDARD  CALIBRATION  CONCENTRATIONS  THE  Y-AXIS  DATA. 

C 

1003  TYPE  1000 

1UO0  FORMATUX, ‘ENTER  THE  NUMBER  OF  S0OIUM  CHLORIDE  CALIBRATION*/ 

C1  CONCENTRATION  STANDARDS',/! 

TYPE  2004 
ACCEPT  2000. NP 
2000  FORMAT t 112) 

C 

C  OUTPUT  THE  INFORMATION  USED  TO  PERFORM  THE  CALCULATIONS  IN 
C  THE  DA TAX. XXX  FILE. 

C 

TYPE  2006 

TYPE  41  37.ICCDSU)  »ICCUS<3> 

TYPE  2006 
KRITEiS, 20061 
200?  FORMATUX.//////) 

NRITEI2.  J03DNP 

3001  FORMAT! 1X.24X, ‘THE  NUMBER  OF  SODIUM  CHLORIDE 
C  CALIBRATION  CONCENTRATION  STANDARDS-  M2) 

WRITE 12. 2006) 

TYPE  3000 

3000  FORMATUX, 'DEPRESS  RETURN  KEY  AFTER  ENTERING  a  VOLTAGE  *  / 

C1  MEASUREMENT'//) 

TYPE  2006 

C 

C  THE  FOLLOWING  INFORMATION  tS  AVAILABLE  FROM  THE  DATA  SHEET 
C  USED  DURING  AN  ROFT  EVALUATION. 

C 
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TYPE  4000 

4000  FORMAT (IX, 'ENTER  THE  DATA  POINTS'/ 

C'  SAMPLE  NUMBER  VOLTAGE  SODIUM  CHLORIDE  / 

r '  MEASUREMENT  CALIBRATION'/ 

r<  IN  VOLTS  CONCENTRATION'/ 

Ct  (x  DATA)  (I  DATA)  '/) 

WRITE( 2 ,3003) 

3003  FORM AT (IX, 2 4X, 'SAMPLE  NUMBER' , 18X ,' VOLTAGE ', 23X ,' SODIUM  CHLORIDE') 
WRITE ( 2 , 3082 ) 

3082  FORMAT (IX , 24X , 13X , 18X , ' MEASUREMENT  * ,19X , ' CALIBRATION' ) 

VJRITE  (2 , 3  083) 

3083  FORMAT (IX , 24X , 13X , 18X , ' (IN  VOLTS) ', 20X , 'CONCENTRATION ' ) 

WRITE( 2 , 3084 ) 

3084  FORMAT ( IX , 24X , 13X , 18X , ' (  X  DATA  ) ' ,20X. ' (  Y  DATA  )  /) 

LINES-0 

DO  5001  1=1, NP 
TYPE  4500,  I 

4500  FORMAT (1H+,T7,H2,T20,1H  ,$) 

ACCEPT "5000 ,  X ( I ) 

CALL  CLEAR(LINES) 

TYPE  4501,  I  ,X ( I ) 

4501  FORMAT( '+' ,T7, 112, T21,1G13.6, '  ',?) 

ACCEPT  5000,  Y ( I > 

TYPE  305 

5000  FORMAT (1G1 3 .6) 

XAS  ( I)  =X  ( I) 

XXAS (I) =X(I) 

WRITE(2, 5002)  1,  X  (I )  ,Y(I) 

5002  FORMAT ( IX, 29X,1 12, 23X,1G13 .6 , 19X ,  1G13 .6) 

5001  CONTINUE 

C  ENTER  THE  DESCRIPTIVE  INFORMATION  CONCERNING  THE  SUBJECT,  MASK, 

C  DATE,  AND  TIME  TESTED. 

C 

TYPE  2006 

TYPE  4137,  ICCBS  (1)  ,  ICCBS  (2) 

TYPE  2006 
3  077  TYPE  3  085 

3085  FORMAT (IX, 'SUBJECT  NAME:',2X,S) 

ACCEPT  3 08 6  , NAME 

3086  FORMAT  ( 4  5A1) 

TYPE  3087 

3  087  FORMAT  (IX, 'TYPE  OF  MASK  s  '  ,  2X,  $> 

ACCEPT  3 086, MASK 
TYPE  3088 

3088  FORMAT (IX, 'DATE  TESTED t ' ,2X,$) 

ACCEPT  3 08 9, DATE 

3089  FORMAT  ( 4 5A1) 

TYPE  3  090 

3090  FORMAT  (IX, 'TIME  TESTED  i  ',  2X ,  $) 

ACCEPT  3089, TIME 

TYPE  2006 
3092  FORMAT (12) 

TYPE  2006 

C  ENTER  INTEGRATION  EXERCISE  COUNT  DATA  FOR  THE  TEST  PROTOCOL. 

C 
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c 

C  THIS  INFORMATION  IS  AVAILABLE  FROM  THE  DATA  SHEET 
C  USED  DURING  AN  RQFT  EVALUATION. 

C 

C 

C  SELECT  THE  EXERCISE  PROTOCOL. 

C 

TYPE  2006 

TYPE  4137 , ICCBS (1 ) , ICCBS (2) 

TYPE  2006 
4139  TYPE  3127 

3127  FORMATdX, 'THE  USER  IS  FREE  TO  SELECT  ONE  OF  TWO  GROUPS  OF'/ 

C*  EXERCISE  PROTOCOLS.'//) 

TYPE  3128 

3128  FORMATdX, 'THE  (GROUP  1]  EXERCISE  PROTOCOL  CONSISTS  OFi '// 

C'  til  NORMAL  BREATHING  STRAIGHT  AHEAD'/ 

C'  (21  DEEP  BREATHING  STRAIGHT  AHEAD'/ 

C*  (31  TALKING'/ 

C'  ( 4 J  SIDE-TO-SIDE  HEAD  MOVEMENTS  (DEEP  BREATHING)'/ 

C'  (5)  UP-AND-DOWN  HEAD  MOVEMENTS  (DEEP  BREATHING) '/ 

C*  (61  FACIAL  GRIMACING'///) 

TYPE  2006 

TYPE  4137,  ICCBS  (1 )  ,  ICCBS  (2) 

TYPE  2006 
TYPE  3129 

3129  FORMATdX, 'THE  (GROUP  2)  EXERCISE  PROTOCOL  CONSISTS  OFi'// 

C*  (1)  NORMAL  BREATHING  STRAIGHT  AHEAD'/ 

C  (21  NORMAL  BREATHING  LEFT'/ 

C'  (31  NORMAL  BREATHING  RIGHT'/ 

C'  (41  NORMAL  BREATHING  DOWN*/ 

C*  (51  NORMAL  BREATHING  UP'/ 

C'  161  DEEP  BREATHING  STRAIGHT  AHEAD'/ 

C'  (7)  DEEP  BREATHING  LEFT'/ 

C'  (81  DEEP  BREATHING  RIGHT’) 

TYPE  3130 

3130  FORMATdX,  '  (91  DEEP  BREATHING  DOWN'/ 

C'  (101  DEEP  BREATHING  UP'/ 

C  (11)  TALXING '/ 

C*  (125  FACIAL  GRIMACING'/ 

C'  (131  SIDE-TO-SIDE  HEAD  MOVEMENTS  (NORMAL  BREATHING)'/ 

C‘  (141  UP-AND-DOWN  HEAD  MOVEMENTS  (NORMAL  BREATHING) '/ 

C  (15)  SIDE-TO-SIDE  HEAD  MOVEMENTS  (DEEP  BREATHING)'/ 

C*  1161  UP-AND-DOWN  HEAD  MOVEMENTS  (DEEP  BREATHING)'//) 

4138  TYPE  2006 

TYPE  4137,  ICCDSd ) ,  ICCBS ( 2) 

TYPE  2006 
TYPE  3131 

3131  FORMATdX, 'TO  SPECIFY  THE  EXERCISE  PROTOCOL  GROUP  OP  INTEREST,'/ 

C'  TYPE  EITHER:  GROUP  1  OR  GROUP  2  ') 

TYPE  3  06 
TYPE  3199 

3199  FORMATdX, 'ENTRY  -  *,S) 

ACCEPT  3122, SELECT 
3122  FORMAT (7A1) 

IF  t SELECT (7) .NS. ' 1 ' .AND. SELECT  (?) .NS. *2 ' )  GO  TO  4138 
TYPE  2006 

TYPE  4137, ICCBSti ), ICCBS (2) 
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TYPE  2006 
TYPE  3021 

3021  FORMAT (IX, 'USER  ATTENTION:  IF  NO  EXERCISE  COUNT  DATA  WAS'/ 

C*  COLLECTED  FOR  A  PARTICUIAR  EXERCISE,  TYPE:  000001.  ALSO,'/ 

C'  FOR  EACH  TYPED  EXERCISE  COUNT  DATA  ENTRY,  SIX  DIGITS  MUST') 

TYPE  3022 

3022  FORMAT (IX, 'BE  TYPED,  THAT  IS,  IF  YOU  HAVE  A  SIX  DIGIT  NUMBER,'/ 

C'  TYPE  ALL  SIX  DIGITS.  IF  YOU  HAVE  A  FIVE  DIGIT  NUMBER,  TYPE'/ 

C'  ONE  LEADING  2ERO  AND  THEN  THE  FIVE  DIGITS.  IF  YOU  HAVE  A') 

TYPE  3033 

3033  FORMAT ( '  FOUR  DIGIT  NUMBER,  TVPE  TWO  LEADING  ZEROS  AND  THEN  THE'/ 
C'  FOUR  DIGITS,  ETC.  SEVERAL  EXAMPLES  FOLLOW  AS  AN  ILLUSTRATION') 

TYPE  3034 

3034  FORMAT (IX, 'FOR  EXAMPLE:  COUNT  DATA=743182  TYPED  ENTRY«743 182 ' ) 
TYPE  3035 

3035  FORMAT (IX, 'FOR  EXAMPLE:  COUNT  DATA=18726  TYPED  ENTRY-018726') 
TYPE  3036 

3036  FORMAT (IX, 'FOR  EXAMPLE:  COUNT  DATA-64 12  TYPED  ENTRY-006412') 
TYPE  2006 

TYPE  41 37  ,  ICCBS  ( 1)  ,  ICCBS  (2) 

TYPE  2006 
TYPE  3119 

3119  FORMATdX, 'DEPRESS  THE  RETURN  KEY  AFTER  ENTERING  AN'/ 

C  INTEGRATOR  COUNT'/! 

TYPE  2006 
TYPE  3037 

3037  FORMATdX, 'EXERCISE  COUNT  DATA:'/) 

TYPE  3038 

3038  FORMATdX,  'EXERCISE*  ,29X,  'INTEGRATOR1  #7X,  'TIME  PERIOD') 

TYPE  3126 

3126  FORMATdX, 40X, 'COUNT' ,9X,  '  (IN  SECONDS)'/) 

IF ( SELECT (7)  .BQ.  *2')  OO  TO  3133 
TYPE  3039 

3039  FORMATdX, 'NORMAL  BREATHING  STRAIGHT  AHEAD  *,S) 

ACCEPT  3040.IC1 

CALL  CLEAR (LINES) 

3040  FORMAT (16) 

TYPE  3112, IC1 

3112  FORMATUH+,  'NORMAL  BREATHING  STRAIGHT  AHEAD '  ,T39 , 16 ,T59, 1H  ,$) 
ACCEPT  3113  ,ICTP1 

3113  FORMAT) 12) 

TYPE  3041 

3041  FORMAT (IX, 'DEEP  BREATHING  STRAIGHT  AHEAD  ',*> 

ACCEPT  3040  , IC2 

CALL  CLEAR i LIMES) 

TYPE  3114, IC2 

3114  FORNATOM*,  'DEEP  BREATHING  STRAIGHT  AHEAD' ,T39, 16, T59, 

•111  ,5) 

ACCEPT  3113 ,ICTP2 
TYPE  3042 

3042  FORMATdX, 'TALKING  ',*) 

ACCEPT  3 040, ICS 

CALL  CLEAR (LINES) 

TYPE  3115, IC3 

3115  FORMAT  dH+*  'TALKING '  ,T39 , 16,T59,  111  ,  S) 

ACCEPT  3113, ICTP3 

TYPE  3043 
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3043  FORMATtlX, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C*  (DEEP  BREATHING) ',20X,1H  ,$> 

ACCEPT  3  040  ,  IC4 
CALL  CLEAR (LINES! 

TYPE  3116, IC4 

3116  FORMAT ( 1H+,  '  (DEEP  BREATHING!  '  ,T39 , 16  ,T59, 1H  ,  $) 

ACCEPT  3113 ,  ICTP4 

TYPE  3044 

3044  FORMAT (IX, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING) ' ,20X,1H  ,$) 

ACCEPT  3  040  ,  IC5 
CALL  CLEAR (LINES) 

TYPE  3116, ICS 
ACCEPT  3113 , ICTP5 
TYPE  3045 

3045  FORMATUX,  'FACIAL  GRIMACING *, 2 OX ,  1H  ,$) 

ACCEPT  3  040  , IC6 

CALL  CLEAR ( LINES) 

TYPE  3117, IC6 

3117  FORMAT (1H+, 'FACIAL  GRIMACING ' ,T39 , 16 ,T59 , 1H  , $) 

ACCEPT  3113 , ICTP6 

TYPE  2006 

TYPE  41  37  ,  ICCBS  ( 1)  ,  ICCBS  ( 2) 

TYPE  2006 
GO  TO  3153 

3133  CONTINUE 
TYPE  3039 
ACCEPT  3040  , IC1 
CALL  CLEAR(LINES) 

TYPE  3112, IC1 
ACCEPT  3113 , ICTP1 
TYPE  3134 

3134  FORMATtlX,  'NORMAL  BREATHING  LEFT  ,S) 

ACCEPT  3040.IC2 

CALL  CLEAR (LINES) 

TYPE  3135, IC2 

3135  FORMAT till*, 'NORMAL  BREATHING  LEFT ‘ ,T39 , 16 ,T59, 1H  , S) 
ACCEPT  3U3.ICTP2 

TYPE  3136 

3136  FORMATtlX, 'NORMAL  BREATHING  RIGHT  ’,S) 

ACCEPT  3040, IC3 

CALL  CLEAR (LINES) 

TYPE  3 137  ,  IC3 

3137  FORMAT!  Ill*  ,  'NORMAL  BREATHING  RIGHT*  ,T39, 16  ,TS9, 1H  ,5) 
ACCEPT  3113, ICTP3 

TYPE  3130 

3130  FORMATtlX, 'NORMAL  BREATHING  DOWN  ',5) 

ACCEPT  3040  ,  IC4 
CALL  CLEAR (LINES) 

TYPE  3133, IC4 

3139  FORMAT  till*,  'NORMAL  BREATHING  DOWN*  ,T39 , 16 ,T59, 11)  ,S) 
ACCEPT  m3,ICTP4 

TYPE  3140 

3140  FORMATtlX, 'NORMAL  BREATHING  UP  \S> 

ACCEPT  3040, ICS 

CALL  CLEAR (LINES) 

TYPE  3 141, ICS 
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3141  FORMAT (1H+, 'NORMAL  BREATHING  UP ' ,T39 , 16 #T59 , 1H  , $) 
ACCEPT  3113, ICTP5 

TYPE  3  041 
ACCEPT  3040, IC6 
CALL  CLEAR (LINES) 

TYPE  3114,  IC6 
ACCEPT  3113, ICTP6 
TYPE  3142 

3142  FORMAT (IX, 'DEEP  BREATHING  LEFT  ',$) 

ACCEPT  3040, IC7 

CALL  CLEAR  <  LINES) 

TYPE  3 143  ,  IC7 

3143  FORMAT (1H+, 'DEEP  BREATHING  LEFT ' ,T39 , 16 ,T59 , 1H  , $) 
ACCEPT  3113,ICTP7 

TYPE  3144 

3144  FORMAT (IX, 'DEEP  BREATHING  RIGHT  ',$> 

ACCEPT  3 040 , IC8 

CALL  CLEAR (LINES ) 

TYPE  3145  ,  IC8 

3145  FORMATUH+,  'DEEP  BREATHING  RIGHT ' ,T39 , 16 ,T59, 1H  ,$) 
ACCEPT  3113 , ICTP8 

TYPE  3146 

3146  FORMAT (IX. 'DEEP  BREATHING  DOWN  ',$) 

ACCEPT  3040.IC9 

CALL  CLEAR(LINES) 

TYPE  3147  ,IC9 

3147  F0RMAT(.1H+,  'DEEP  BREATHING  DOWN*  ,T39 , 16  ,T59 , 1H  ,  $) 
ACCEPT  3113,ICTP9 

TYPE  3148 

3148  FORMAT (IX, 'DEEP  BREATHING  UP  * , S) 

ACCEPT  3040, IC10 

CALL  CLEAR(LINES) 

TYPE  3149  IC10 

3149  FORMAT(IH+, ' DEEP  BREATHING  UP ' ,T39 , 16 ,T59 , 1H  , S) 
ACCEPT  3113.ICTP10 

TYPE  3042 
ACCEPT  3040, IC11 
CALL  CLEAR (LINES) 

TYPE  3115, IC11 
ACCEPT  3113,ICTV11 
TYPE  3045 
ACCEPT  3040, IC12 
CALL  CLEAR (LINES) 

TYPE  3117, IC12 
ACCEPT  3113.ICTP12 
TYPE  3150 

3150  FOPNATtlX, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING)  ',S> 

ACCEPT  3040, IC1 3 
CALI.  CLEAR  (LINES) 

TYPE  3151, IC13 

3151  FORMATUH+,'  (NORMAL  BREATHING)  '  ,T39 , 16 ,T59, III  ,S> 
ACCEPT  3113,ICTP13 

TYPE  3152 

3152  FORKATUX.'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  t NORMAL  BREATHING)  ',*) 

ACCEPT  3040, IC14 
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CALL  CLEAR(LINES) 

TYPE  3151, XC14 
ACCEPT  3113 , ICTP14 
TYPE  3043 
ACCEPT  3040  ,IC15 
CALL  CLEAR (LINES) 

TYPE  3116, IC15 
ACCEPT  3113 ,ICTP15 
TYPE  3044 
ACCEPT  3040  ,IC16 
CALL  CLEAR (LINES) 

TYPE  3116, IC16 
ACCEPT  3113 ,ICTP16 
TYPE  2006 

TYPE  4137 , ICCBS (1 ) ,ICCBS(2) 

3153  TYPE  2006 

WRITE(2, 2006) 

WRITE(2, 3046)  NAME 

3046  FORMAT(6X, 'SUBJECT  NAME! ' ,2X,45A1) 

WRITE (2, 3 047) MASK 

3047  FORMAT(6X, 'TYPE  OF  MASK !  ’ , 2X ,  4 5A1) 

WRITE(2 , 3048)  DATE 

3048  FORMAT(6X, 'DATE  TESTED  !',2X,45A1) 

WRITE (2, 3 049) TIME 

3049  FORM AT (6X, 'TIME  TESTED  t',2X,45Al) 

WRITE(2 ,2006) 

WRITE(2,9006) 

9006  FORHAT(6X, 'EXERCISE  INTEGRATOR  COUNT  DATAt '/) 

WRITE(2 ,3094) 

3094  FORMAT(6X, 'EXERCISE' ,29X, 'INTEGRATOR  COUNT' ,7X» 'TIME  PERIOD') 
WRITE(2 , 3120) 

3120  FORMAT(6X,8X,29X,16X,7X, ' (IN  SECONDS)',/) 

IF(SELECT(7)  .EQ.  ' 2  * )  GO  TO  3154 

XC1-IC1 

XC2-IC2 

XC3-IC3 

XC4-IC4 

XC5-IC5 

XC6-IC6 

IC(1)»XC1 

IC ( 2) »XC2 

IC(3)«XC3 

IC(4)-XC4 

IC(5)«XC5 

IC(6)-XC6 

icTPdxem 

ICTP  (2)  ■ICTP2 

ICTP(3) "ICTP3 

ICTP(4)-ICTP4 

ICTP(5) "ICTP5 

ICTP(6) "ICTP6 

WRITE(2,3051> IC1 , ICTP1 

3051  F0RMAT(6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD' ,1 IX. 16, 17X, 12) 
WRITE (2, 3 052) IC2 , ICTP2 

WRITE(2, 3053) IC3,ICTP3 

3052  FORHAT(6X, 'DEEP  BREATHING  STRAIGHT  AHEAD' ,13X,I6,17X,I2) 

3053  FORHAT(6X, 'TALKING', 35X, 16, 17X, 12) 
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WRITER  2 , 3  054) IC4,ICTP4 

3054  FORMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING) ’ ,26X, 16, 17X, 12) 

WRITE(2,3055) IC5,ICTP5 

3055  FORMAT(6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING)  '  ,26X,  16, 17X,  12) 

WRITE(2,3056) IC6,ICTP6 

3056  FORMAT(6X, 'FACIAL  GRIMACING ' ,26X , 16 ,17X , 12//) 
GO  TO  3155 

3154  CONTINUE 
XC1=IC1 
XC2-IC2 
XC3=IC3 
XC4»IC4 
XC5*=IC5 
XC6-IC6 
XC7-IC7 
XC8-IC8 
XC9-IC9 
XC10-IC10 
XC11-ICH 
XC12-IC12 
XC13-IC13 
XC14-IC14 
XC15»IC15 
XC16-IC16 
IC(1)«XC1 
IC(2)»XC2 
IC(3)»XC3 
IC ( 4) »XC4 
IC(5)«XC5 
IC(6) »XC6 
IC(7).XC7 
IC(8)«XC8 
IC(9)-XC9 
IC(10>-XC10 
lCtll)«XCU 
IC  1 12) *XC1 2 
IC(13)«XC13 
IC(14)-XC14 
IC(1S)-XC15 
IC(16)-XC16 
ICTPtU-ICTPl 
iCTP  ( 2)  ■  ICTP2 
ICTP  (3)  “ICTP3 
ICTP < 4)  “ICTP4 
ICTP(S) "ICTPS 
ICTP<6) ■ICTP6 
ICTP  (7)  «ICTP7 
ICTP(8)«ICTP8 
ICTP  (9)  ■ICTP9 
ICTP ( 10) “ICTPIO 
ICTP  Ul)  "ICTP  11 
ICTP  (12) -ICTP  12 
ICTP  (13) -ICTP  13 
ICTP (14) “ICTP 14 
ICTPUSI-ICTPIS 
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ICTP (16) =ICTPI6 
WRITE(2/3051) IC1,ICTP1 
WRITE(2,3156) IC2,ICTP2 

3156  F0RMAT(6X, 'NORMAL  BREATHING  LEFT ' ,21X , 16 , 17X, 12) 

WRITE (2,3157) IC3,ICTP3 

3157  FORMAT(6X, 'NORMAL  BREATHING  RIGHT ' ,20X , 16 ,17X, 12) 

WRITE  (2, 3158)IC4, ICTP4 

3158  FORMAT ( 6 X, 'NORMAL  BREATHING  DOWN* ,21X, 16 , 17X , 12) 

3159  F0RMAT(6X, 'NORMAL  BREATHING  UP' , 23X, I6,17X, 12) 

WRITE (2/3159)105/1 CTP5 

WRITE (2,3052) IC6/ICTP6 
WRITE(2 , 3160) IC7/ICTP7 

3160  FORMAT(6X, 'DEEP  BREATHING  LEFT ' , 23X , 16 , 17X , 12) 

WRITE(2 , 3161) IC8, ICTPO 

3161  FORMAT(6X, 'DEEP  BREATHING  RIGHT ' , 22X , 16 , 17X, 12) 

WRITE(2, 3162)109, ICTP9 

3162  FORMAT(6X, 'DEEP  BREATHING  DOWN' ,23X, I6,17X, 12) 

WR1TE<2,3163) IC10/ICTP10 

3163  F0RMAT(6X, 'DEEP  BREATHING  UP ' ,25X , Ifi ,17X , 12) 

WRITE (2 ,3053 ) IC11 , ICTP11 

WRITE(2,3164) IC12,ICTP12 

3164  FORM AT (6X, 'FACIAL  GRIMACING ' ,26X , 16 , 1 7X , 12) 

WRITE (2, 3165)1013, ICTP13 

3165  F0RMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C'  (NORMAL  BREATHING) *,24X, 16, 17X, 12) 

WRITE(2,3166) IC14 , ICTP14 

3166  FORM AT < 6 X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (NORMAL  BREATHING)  *,24X, 16, 17X, 12) 

WRITE (2, 3 054) IC15,ICTP15 
WRITE (2, 3055) IC16,ICTP16 
3155  REP* ' N ' 

3888  IF ( REP. EQ.NO)  GO  TO  3999 
GO  TO  9983 

C 

C  CALCULATE  THE  AVERAGE  PHOTOMULTIPLIER  TUBE  OUTPUT  VOLTAGE  FOR 
C  EACH  EXERCISE 
C 

3999  CONTINUE 

IF (SELECT (7)  . EQ.  * 2 * >  GO  TO  3167 
DO  3093  1*1,6 

IC(J)*ICtI)/(l 000 .0* ICTP ( I ) ) 

3093  CONTINUE 
GO  TO  3168 

3167  CONTINUE 

DO  3169  1*1,16 

IC( I )*1C( I ) / (1000.0*ICTP( I ) ) 

3169  CONTINUE 

3168  CONTINUE 
END  FILE  2 

9983  CONTINUE 
REWIND  2 

6001  READ(2,2006> 

READ (2, 3001) HP 
3002  FORMAT(lX«52X, 12) 

C 

C  AFTER  ALL  PERTINENT  INFORMATION  IS  ENTERED  INTO  THE  DATA. XXX  FILE, 

C  THE  FILE  IS  REWOUND  AND  THE  POINTER  MOVED  TO  READ  THE  DATA  NEEDED 
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C  TO  CALCULATE  THE  LEAST-SQUARES  POLYNOMIAL  FUNCTION. 

C 

READ (2, 2006) 

READ (2,3003) 

READ ( 2 ,3082 ) 

READ (2,3083) 

READ (2,3084) 

C 

C  TO  MAKE  THE  LEAST-SQUARES  POLYNOMIAL  CURVE  FIT  ROUTINE  AS 

C  FLEXIBLE  AS  POSSIBLE,  THE  USER  CAN  SELECT  A  VARIETY  OF 

C  FUNCTIONAL  DEFINITIONS  FOR  THE  VARIABLE  (X). 

C 

2013  TYPE  2006 

TYPE  4137.ICCBSU  >  ,ICCBS(2) 

TYPE  2006 
TYPE  2  009 

2009  FORMAT (IX, 'THE  USER  IS  FREE  TO  SELECT  ANY  ONE  OF  THE  FOLLOWING'/ 
C'  FUNCTIONAL  DEFINITIONS  OF  THE  VARIABLE  (X)./> 

TYPE  306 
TYPE  2010 

2010  FORMATdX,  '  tlJ  (X)  ’/ 


C' 

(21 

1/ (X)  '/ 

C’ 

(3) 

EXPONENTIAL 

(X)  '/ 

C* 

141 

EXPONENTIAL 

(-X)  '/ 

C' 

(51 

EXPONENTIAL 

d/x) '/ 

C' 

(61 

EXPONENTIAL 

(-1/X)  •/ 

C' 

(71 

NATURAL  LOG 

(X)  '/ 

C' 

(81 

LOG  BASE  TEN 

(X!  '/ 

C' 

(91 

SIN  (X)  '/) 

TYPE  306 
TYPE  2  011 

2011  FORMATdX, 'TO  SPECIFY  A  FUNCTIONAL  DEFINITION  OF  (X),  SELECT'/ 
C'  THE  CORRESPONDING  NUMBER  INSIDE  THE  BRACKETS'/) 

TYPE  306 
TYPE  2C04 
ACCEPT  20)2, IAX 

2012  FORMATd ID 
TYPE  2006 

TYPE  4137,  ICCBS  (D  ,  JCCOS  (2) 

TYPE  2006 

IPtIAX.GT.0. AND.IAX.LT. 10)  GO  TO  2014 
GO  TO  2013 

C 

C  DOUBLE  PRECISION  CALCULATIONS  ARE  HADE  TO  RESOLVE  PROBLEMS 
C  ENCOUNTERED  WHEN  DEALING  WITH  DATA  THAT  RANGES  OVER  SEVERAL 
C  ORDERS  OF  MAGNITUDE. 

C 

C 

C  THE  FOLLOWING  'DO  LOOP'  CALCUIATE5  AN  EOCI VALENT  VALUE 

C  FOR  THE  (X)  VARIABLES  BAS UD  ON  THE  SS,SCTa>  FUNCTIONAL 

C  DEFINITION. 

C 

2014  DO  300  t«l,NP 

K-l 

READ ( 2,S002!K,X(K) ,¥UO 
IFdAX.EO.l)  X(I)«X(I) 

IFdAX.EQ.2)  X(I)«1.0/XtI> 
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X  ( I ) =DEXP ( 1 . 0/X ( I ) ) 

X ( I >  «DEXP (-1 .  0/X ( I ) ) 


NAME 

MASK 

DATE 

TIME 


IF ( IAX.EQ. 3)  X(I)»DEXP(X(I) ) 

IF ( IAX.EQ. 8)  X ( I )  *>DL0G1 0 ( X  ( I ) ) 

IF ( IAX. EQ.7)  X ( I ) =DLOG  (X  l I ) ) 

IF ( IAX. EQ. 9)  X ( I > »DSIN(X  < I ) ) 

IF ( IAX. EQ. 4)  X ( I ) -DEXP ( -X ( I ) ) 

IF ( IAX. EQ. 5) 

IF ( IAX. EQ. 6) 

300  CONTINUE 

IF(REP.EQ.YES)  GO  TO  9777 
GO  TO  9888 
9777  READ < 2 * 2006) 

READ ( 2 ,3046 ) 

READ ( 2,3047 ) 

READ ( 2,3048) 

READ  ( 2 ,3049) 

READ (2,2006) 

READ(2,9006) 

READ(2,3094) 

READ  ( 2 ,3120) 

IF(SEL£CT (7)  .EQ.  '2')  GO  TO  3170 

READ (2, 3 051)  ICl.ICTPl 

READ< 2 ,3052)  IC2,ICTP2 

READ ( 2 ,3053 )  IC3,ICTP3 

READ ( 2 ,3054)  IC4.ICTP4 

READ(2,3055)  IC5.ICTP5 

READ (2, 3056)  IC6,ICTP6 

XC1-IC1 

XC2-IC2 

XC3-IC3 

XC4-TC4 

XCS-IC5 

XC6-IC6 

IC(l)«XCl 

IC(2)-XC2 

ICC3J-XC3 

IC(4)«XC4 

IC(5)-XC5 

:C(6)*XC6 

ICTP(1)«ICTP1 

ICTP (2) -ICTP2 

ICTP(3) -ICTP! 

1CTP(4)-ICTP4 
ICTP(S) -ICTPS 
ICTP (6) ■ICTP6 
DO  9417  1-1,6 

tCil)-IC<l>/(1000.0*lCTP(X> ) 

9417  CONTINUE 
GO  TO  3171 

1170  READ (2, 1051) IC1»1CTP1 
R£AD(2,31S6) IC2,ICTP2 
READ(2,3157) ICl.ICTPl 
R£AD(2,11S8) IC4, ICTP4 
READ (2, 1159) ICS, ICTPS 
READ (2, 1052) ICS.IC^PO 
READ (2, 11 60! IC7,ICTP7 
READ (2,1161) IC8 » ICTP8 
READ (2, 1162)  IC9,  ICTP9 
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READ ( 2 ,3163) ICIO.ICTPIO 

READ (2,3053)IC11,ICTP11 

READ (2,3164)1012, ICTP 12 

READ (2, 3 165) IC13.ICTP13 

READ (2, 3 166)  IC14  ,  ICTP14 

READ (2, 3 054) IC15,ICTP15 

READ (2, 3055) IC16,ICTP16 

XC1«IC1 

XC2-IC2 

XC3-IC3 

XC4-IC4 

XC5-IC5 

XC6»IC6 

XC7-IC7 

XC8-IC8 

XC9-IC9 

XClOaXCIO 

XC11-1C11 

XC12-XC1-2 

XC13-IC13 

XC14-IC14 

XCIS-ICIS 

XC16-IC16 

ICTP ( 1 ) «ICTP1 

ICTP (2) "ICTP2 

ICTP (3 ) "ICTP3 

ICTP(4) ■ICTP4 

ICTP  C  5) •ICTPS 

ICTP (6) "ICTP6 

ICTP ( 7) ■ ICTP’ 

ICTP  t 8) »ICTP8 

ICTP(9)«ICTP9 

ICTP<10)»ICTP10 

ICTP (11) -ICTP 11 

ICTP(12)*XCTP12 

ICTP(13)-ICTP13 

ICTP(14)*ICTP14 

ICTP(1S)*ICTP1S 

ICTP(16) -ICTP16 

icm-xci 

1C(2)«XC2 

icn )  »xcs 

ICUI-XC4 

IC(S)«XC$ 

IC«6)-XC6 
IC ( 7) •XC7 
IC 1 8) «XC8 
IC(9)«XC9 
ICU0)«XC10 

ic  ( 1 1)  «xei  i 

IC(12)*XC12 

1C  (13)  -XC11 

tCU4)«XC14 

IC(15)»XC1S 

IC(  16)  »XC16 

00  3172  1*1,16 

IC(  I)*ICtl)/(1000.0*lCTP( I) ) 
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3172  CONTINUE 
3171  CONTINUE 
REWIND  2 
9888  TYPE  2006 
306  FORMAT (IX, 1  '/) 

THE  USER  NOW  SELECTS  THE  ORDER  OF  THE  POLYNOMIAL  CURVE 
FIT  EQUATION.  THE  MAXIMUM  ORDER  IS  ONE  LESS  THAN  THE 
NUMBER  OF  (X)  DATA  VALUES. 

TYPE  2006 

TYPE  4137,ICCBS(1> ,ICCBS(2) 

TYPE  2006 
TYPE  301 

301  FORMAT (IX, ‘ENTER  THE  ORDER  OF  THE  DESIRED  POLYNOMIAL’/) 

TYPE  306 
TYPE  302 

302  FORMAT ( IX, ' (MAXIMUM  -  THE  NUMBER  OF  SCCIUM  CHLORIDE'/ 

C‘  CALIBRATION  CONCENTRATION  STANDARDS  -  1) V) 

TYPE  306 
TYPE  2004 
ACCEPT  303.M 

303  FORMAT(I2) 

TYPE  2006 

TYPE  4137 , ICCBS(l) , ICCBS(«) 

TYPE  2006 
DO  403  JK-1,14 
W(JK-l) -0.0 
J(JK-1)»0.0 
403  CONTINUE 
W(0)-NP 
DO  401  I-l,NP 
DO  400  4-1, 2*M 
WU)-W(4)  ♦*(!>  *»J 
S(J)«2(J)*Y(X)*X(I)»*J 

400  CONTINUE 

J ( 0!  (0! *Y( I ) 

401  CONTINUE 
DO  300  I- 1 ,M*1 
DO  S00  4-1, M*1 
0(1 »4)«Wt4“2*l) 

B(  I  ,M*2)»®(I“1) 

S0C  CONTINUE 

C  THE  FOLLOW  I  NO  SERIES  OF  NEST  (2)  DO  LOOPS  ACCOMPLISH  THE  FOLLOWING! 

C  THE  FIRST  INCREMENTS  THE  NESTED  DO  LOOP  SYSTEM  AND  THUS  SOLVES 

C  FOR  EACH  POLYNOMIAL  COEFFICIENT  ON  A  SUCCESSIVE  OASIS, 

C 

C  THE  SECOND  INCREMENTS  THE  NESTED  DO  LOOP  SYSTEM  AND  THUS  DECREASES 
C  THE  MATRIX  SUE  OY  ONE  ROW  AND  COLUMN  EACH  ITERATON. 

C 

C  THE  THIRD  AMD  FOURTH  TAKE  THE  STH  ROW  AND  DIVIDE  IT  BY  ITS  NTH 

C  COEFFICIENT  AND  THEN  MULTIPLIES  THE  ITH  MOW  BY  THE  NTH  COEFFICIENT 

C  OF  THE  (I*1)TH  ROW.  FINALLY  IT  SUBTRACTS  THE  (I-1ITH  ROW  FROM  THE 
C  ITH  ROW. 

C 

C  THE  LA5T  TWO  DO  LOOPS  RESTRUCTURE  MATRIX  O  SO  THAT  IT  CONTAINS  THE 
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E  MATRIX  MINUS  THE  NTH  COLUMN  AND  THE  LAST  RCM. 

DO  50  N=1,M+1 
DO  30  11=1, M+l 
DO  30  Jl=l,M+2 
D(I1,J1)=B(I1,J1) 

30  CONTINUE 

DO  40  K=1,M 
IF ( N-K) 35,35,33 
33  L=  1 

GO  TO  37 
35  L*2 

37  DO  31  1  =  1  ,M+1-K 
DO  31  J“1 ,M*3-K 

E(I,J)=(D(I.J>/D(I,L) >*D(1-H,L) -DO  +  l.J) 

31  CONTINUE 

THE  COEFFICIENTS  FOR  THE  POLYNOMIAL  LEAST-SQUARES  CURVE  FIT  EQUATION 
ARE  CONTAINED  IN  THE  'A(N) *  ARRAY. 

DO  40  1*1 ,H+1-K 
DO  4Q.  J“i ,M+2-E 
IF(J  .GE.  L)  GO  TO  45 
D ( 1 , J ) *E ( I , J) 

G©  TO  40 

45  D(I,J)-E(I,J+1> 

«Q  CONTINUE 

A ( N)  =D  ( 1 , 2)  /D  ( 1 , 1 ) 

50  CONTINUE 
WRITEO  ,2007) 

WRITEO, SUM 

OUTPUT  THE  CALCUIATQ3  ^FORMATION  INTO  THE  CALC X.XXX  FILE. 

51  FORMAT! IX, 24X,  ’THE  ORDER  OP  THE  DESIRED  POi-VNGHIAL*  MU 
WRITE! J. 2006) 

WRITE! 1 ,20061 
IP! IAX.EQ.il  GO  TO  52 
IF! IAX.FQ.25  GO  TO  S3 
IFUAX.KQ.J)  GO  TO  54 
lFtIAX.EQ.45  GO  TO  55 
tmAX.EO.55  GO  TO  56 
tEtlAX.EQ.65  GO  TO  57 
Sft IAX.EQ.7)  GO  TO  56 
If  I IAX.EQ.8)  GO  TO  59 
IF  tlAX.8Q.95  CO  TO  $0 

52  WRITE! 3,615 

61  FORMAT t IX , I4X ,  'THE  POLYNOMIAL  FUNCTIONAL  DEFINITION  OP  THE 
e  VAR  I  ABLE  t  SO  IS  IN  TERMS  OF  t  130*///) 

GO  TO  99 
57  WRITEO  ,62  5 

62  PORHATUX.IAX,  'THE  POLYNOMIAL  FUNCTIONAL  DEFINITION  OF  THE 
C  VARIABLE  <X>  IS  IN  TERNS  OFi  C1/X5V//5 

GO  TO  89 
54  WRITEO, 635 

63  FORM AT  1 1 X , 2  4X , ’ THE  POLYHCH  I AL  FUNCTIONAL  DEFINITION  OF  THE 
C  VARIABLE  !XI  IS  IN  TERMS  OP:  <  EXPONENT  SAL  IX5IV//J 
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GO  TO  89 

55  WRITE! 3  ,  64 ) 

64  F0RMAT(1X,24X ,  'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

56  WR I T  E ( 3 , 6  5 ) 

65  FORMAT! IX , 24X , 'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

57  WRITE13 ,66) 

66  FORMAT! IX, 24X, 'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

58  WRITE(3 ,67) 

67  FORMAT(iX,24X, 'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

59  WRITE* 3 ,68) 

68  F0RMAT!1X,24X, 'THE  POLYNOMIAL 


FUNCTIONAL  DEFINITION  Cc  THE 
(EXPONENTIAL  <-X) )'///) 


FUNCTIONAL  DEFINITION  OF  THE 
(EXPONENTIAL  U/X))V//> 


FUNCTIONAL  DEFINITION  OF  THE 
(EXPONENTIAL  (-1/X) )*///> 


FUNCTIONAL  DEFINITION  OF  THE 
{NATURAL  LOG  (X)  '///> 


FUNCTIONAL  DEFINITION  OF  THE 


C  VARIABLE  (X)  IS  IN  TERMS  OF:  (LOG  BASE  TEN  (X)'///) 

GO  TO  89 
60  WRITE(3 ,69) 

69  F0RMAT(1X,24X, 'THE  POLYNOMIAL  FUNCTIONAL  DEFINITION  OF  THE 
C  VARIABLE  (X)  IS  IN  TERMS  OF:  (SIN  (X) )'///) 

GO  TO  89 

89  CONTINUE 

00  90  I-I  ,M+l 
TYPE  TO ,  I  ,A ( I ) 

71  FORMAT! IX , 24X , ’COEFFICIENT  NUMBER', 13,'-  '1PE1S.5) 

70  FORMAT!!*,  'COEFFICIENT  NUMBER* ,  1 3,  '•  ‘.1PE15.5) 

WR7TFQ  , 71)  I  ,A(  I ) 

90  CONTINUE 
TV  ft  2006 

TYPE  4137, ICCBSt I) , ICCBSt 2> 

TYPE  2006 
TYPE  2006 
IS  TYPE  19 

19  FOSMA?US*'DO  YOU  WANT  A  LIST  OP  RESIDUALS  ?'/ 

C*  (AN&iSR  YES  OR  NO)  *,2X,S> 

ACCEPT  5  v!)  2,  COM  MEN 
TYPE  2036 

IPfvOMHEN  ,SQ.  YES)  GO  TO  18 
IFteOMM&N  NO)  GO  TO  17 

SPUrOKMSN  .tit.  n.S  .AND.  COMMEN  ,Nt.  NO)  GO  TO  16 
TYPE  2006 

TYPE  4137 , iCCBstl) ,ICCBS(2) 

TYPE  *006 
18  wpmo.200?) 

WRITE tl, 23) 

23  FORMAT {1X>24P, 'THE  RESIDUALS  ARE  CALCULATED  BY  THE  FOLLOWING') 
WRITES)! ,  )  3) 

13  FORM  AT  (IX,  2  4X,  'EQUATION:  tY(DATA)-Y(ESTIMATED) )  V/) 
WaiTEn,2G06> 

WS1TE(3.24) 

24  FORMAT! IX, I BX, 'SAMPLE  NUMBER* , 17X, 'YOATA' ,23X, 'Y (ESTIMATED) ' , 
C»9X, 'RESIDUAL*/) 
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THE  RESIDUALS  ARE  CALCULATED  USING  THE  POLYNOMIAL  CURVE 
FIT  EQUATION  AND  THE  ORIGINAL  DATA. 

SUMRES=0 . 0 
DO  15  K=1,NP 
THEORY=0 . 0 
DO  21  1=1, M+I 

IF(I  .EQ.  1  .AND.  X(K)  .EQ.  0.0)  GO  TO  20 
THEORY=A ( I ) *X(K>  **(I-1) +THEORY 
GO  TO  21 

20  THEORY=A(I) *1.0** (1-1) +THEORY 

21  CONTINUE 

CALCULATE  THE  SUM  OF  THE  SQUARES  RESIDUAL. 

RESID=Y(I<)  -THEORY 
SRESD=R!oIC ‘RESID 
SUMRES=i.  .Eld  f  SUMRES 
TYPE  10, K,-.  .SID 

10  FORMAT (IX, 'SAMPLE  NUMBER  ' , 13 , 12X , ' RESIDUAL-  ' , 1PE15 . 5) 
WRITE(3,22)K,Y(K) , THEORY, RESID 

22  FORMAT(1X,23X,I2,20X,1G13.6,18X,1G13.6,16X,1G13.6) 

15  CONTINUE 

TYPE  2  006 
TYPE  7 2, SUMRES 

72  FORMAT (IX, 'SUM  OF  THE  SQUARE  RESIDUALS*  ' rlPElS.5> 

TYPE  2007 

WRITE(3 ,2007) 

WRITE(3  ,73)  SUMRES 

73  FORMAT (IX, 2 4X  ,  'SUM  OF  THE  SQUARE  RESIDUALS’*  MPE15.5) 

17  CONTINUE 

CALCULATE  INDIVIDUAL  EXERCISE  PROTECTION  FACTORS 


NOTE:  SINCE  THE  LEAK  MEASURING  SENSITIVITY  OF  THE  SODIUM 
CHLORIDE  RQFT  INSTRUMENT  IS  ONE  PART  IN  TEN  TO  THE 
SIXTH,  ANY  EXERCISE  SCALED  INTEGRATOR  COUNT  VALUE 
YIELDING  A  PF  GREATER  THAN  1.0E+06,  WILL  BE 
REPORTED  AS  1.0E+06.  REPORTING  A  MASK  LEAKAGE  GREATER 
THAN  1.0E+06  WOULD  BE  ERRONEOUS.  ANY  EXERCISE  SCALED 
INTEGRATOR  COUNT  VALUE  YIELDING  A  PROTECTION  FACTOR 
GREATER  THAN  1.0  WILL  BE  REPORTED  AS  1.0E+00. 

TYPE  2006 

TYPE  41 37 , ICCBS ( 1 ) ,ICCBS(2) 

TYPE  2006 

IF (SELECT (7)  .EQ.  '2')  GO  TO  3173 
IDL=6 
IDLP=7 
GO  TO  3174 

3173  CONTINUE 
IDL=16 
IDLP=17 

3174  CONTIN'  E 
PFEST<  J . 0 

DO  9173  1=1, IDL 
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PF ( I ) =0 . 0 
9173  CONTINUE 

CALCPF=0 . 0 
DO  3  057  1  =  1,  IDL 

IF(IC(I)  .LE.XXAS  (7)  )  GO  TO  9227 
IF(IC(I) .GE.XXAS (1) )  GO  TO  4136 
IFdAX.EQ.l)  PFEST  =  IC(I) 

IF(IAX.EQ.2)  PFEST=1 .0/10(1) 

IFdAX.EQ.3)  PFEST=DEXP  (IC  (I)  ) 

IF ( IAX. EQ. 4)  PFEST=DEXP(-IC(I) ) 

IF(IAX.EQ.5)  PF EST=DEXP  (1 . 0/ IC  ( I ) 

IF ( IAX. EQ. 6)  PFEST=DEXP(-1.0/IC(I) ) 

IF (IAX. EQ. 7)  PFEST=DLOG(IC(I> ) 

IF (IAX. EQ. 8)  PFEST=DLOG10( IC (I ) ) 

IF (IAX. EQ. 9)  PFEST=DSIN(IC(I) ) 

CALCPF=0 . 0 
GO  TO  9228 

9227  CALCPF=0. 000001 
GO  TO  3058 

4136  CALCPF=0 . 0 
GO  TO  3058 

9228  DO  3058  J=1,M+1 

IF ( J. EQ. 1 .AND. PFEST. EQ. 0. 0)  GO  TO  3059 
CALCPF=A ( J) *PFEST*  * ( J  —  1 > +CALCPF 
GO  TO  3058 

3059  CALCPF=A (J)  *1.0** (J-l)  +CALCPF 
3058  CONTINUE 

C 

C  BY  DEFINITION,  THE  PROTECTION  FACTOR  (PF)  IS: 

C 

C  PF  =  (CHALLENGE  ATMOSPHERE  CONCENTRATION) / (MASK  LEAK  CONCENTRATION) 

C 

C 

IF (CALCPF , NE , 0 . 0)  GO  TO  4135 
PF  (I )  =0 . 0 
GO  TO  3  057 

4135  PF  (I )  =1 . 0/ CALCPF 
PF(I) “1.0/ CALCPF 
3057  CONTINUF 
C 

C  CALCULATE  AN  OVERALL  ARITHMETIC  AVERAGE  PROTECTION  FACTOR  FOR  ALL 
C  EXERCISES. 

C 

KOUNT=0 
PFSUM=0 . 0 
DO  3060  MT= 1 , IDL 
KOUNT=KOUOT  + 1 
PFSUM=PF:7'!'-’.  +  PF(MT) 

3060  CONTINUE 

PF ( IDLP) =PFSUM/KOUNT 
C 
C 

C  CALCULATE  AN  OVERALL  TIME  WEIGHTED  AVERAGE  PROTECTION  FACTOR 
C  FOR  ALL  EXERCISES. 

C 

C 

WPF=0 . 0 
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KKOUMT=0 
PPSUM=0 . 0 
DO  3194  1MT=1,IDL 
KKOUNT=KKOUNr  +  ICTP(IMT) 

PPSUM=PPSUM  +  (  PF  ( I WTS  * ICTP ( IMT) ) 

3194  CONTINUE 

WPF=PPSUM/KKOUNT 
WRITE(3 ,3061) 

3061  FORMAT(lHl) 

WRITEO  ,3062) 

3062  FORMAT(6X, 'THE  DESCRIPTIVE  AND  PROTECTION  FACTOR  CALCULATIONS:'/) 
WRITEO  ,9542) 

9542  FORMAT(6X, 'NOTE:  ANY  PROTECTION  FACTOR  THAT  IS  LISTED  AS'/ 

C'  1.0E+06  HAS  BEEN  ASSIGNEE  THIS  VALUE  BY  DEFAULT'/ 

C'  BECAUSE  THE  SENSITIVITY  OF  THIS  RQFT  INSTRUMENT  IS'/ 

C'  AT  MOST  ONE  PART  IN  TEN  TO  THE  SIXTH.  THE  INTEGRATOR'/ 

C'  COUNT  VALUE  FOR  A  PARTICULAR  EXERCISE  IN  QUESTION’/ 

C'  IS  MERELY  REPRESENTATIVE  OF  INTEGRATING  THE  ELECTRICAL'/ 

C'  NOISE  AND  THE  TRUE  PROTECTION  FACTOR  IS  INDEED  LESS') 

WRITE(3 ,9543 ) 

9543  FORMAT(6X, 'THAN  1.0E+06.  ANY  EXERCISE  SCALED  INTEGRATOR'/ 

C'  COUNT  VALUE  YIELDING  A  PROTECTION  FACTOR  GREATER'/ 

C'  THAN  1.0  WILL  BE  REPORTED  AS  0.0E-01.'//) 

WRITE(3 , 3176) 

3176  FORMAT (IX, '  '//) 

TYPE  3  062  / 

TYPE  2006 
TYPE  9542 
TYPE  9  54"* 

TYPE  2005 

TYPE  41 37 , ICCBS (1)  ,ICCBS(2) 

TYPE  2006 
WRITE(3 ,3046)  NAME 
WRITE (3, 3 047) MASK 
WRITE (3 , 3048) DATE 
WRITE(3,3049)TIME 
WRITE  (3 ,3176) 

TYPE  3 046  , NAME 
TYPE  3047, MASK 
TYPE  3 04 8, DATE 
TYPE  3049, TIME 
TYPE  2006 

TYPE  41 37 , ICCBS ( 1 ) ,ICCBS(2) 

TYPE  2006 
WRITE(3 ,3063) 

TYPE  3063 

3063  FORMAT(6X, 'EXERCISE' ,29X, 'PROTECTION  FACTOR'/) 

IF ( SELECT (7)  .EQ.  '2')  GO  TO  3175 

WRITE(3 ,3064)  PF (1 ) 

TYPE  3  064  ,  PF  ( ) ) 

3064  FO RM AT (6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD ' ,6X , 1PE1 2 . 1) 

WRITE!? ,3065)  PF  < 2 > 

TYPE  3  065  ,  PF  (  2) 

3065  FORM  AT  (6X, 'DEEP  BREATHING  STRAIGHT  AHEAD '  ,  8X ,  1PE1 2 . 1) 

WRITE ( 3  , 3 066)  PF(3) 

TYPE  3  066  ,PF  !3) 

3  066  FORMAT(6X,  'TALKING  ' , 3 OX ,  1PE1 2 . 1) 
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WRITE(3 ,3067)  PF  <  4 ) 

TYPE  3067, PF(4) 

3067  FORMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS’/ 

C'  (DEEP  BREATHING) ' ,  21X , 1PE12 . 1) 

WRITE  (3 ,3  068)  PF(5) 

TYPE  3068,  PF  (  5) 

3068  FORM AT <6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING)  '  ,2 IX ,  1PE12 . 1) 

WRITE(3 ,3069)  PF(6) 

TYPE  3069, PF(6) 

3069  FORMAT(6X, 'FACIAL  CRI MACING ', 21X , 1PE1 2 . 1//// ) 

WRITE(3 ,3070)  PF  (7) 

TYPE  3  070  ,PF  ( 7) 

3070  FORMAT(lHO,5X, 'OVERALL  ARITHMETIC  AVERAGE  PROTECTION  FACTOR  '/ 

C'  FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMED  =  ' 

C, 1PE8 .1) 

TYPE  2006 
WRITE(3 ,2006) 

TYPE  3195, WPF 
WRITE (3,3195)  WPF 

3195  FORMAT(lHO,5X, 'OVERALL  TIME  WEIGHTED  AVERAGE  PROTECTION  FACTOR  '/ 
C'  FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMED  =  ' 

C, 1PE8 . 1 ) 

GO  TO  3177 
3175  CONTINUE 

WRITE(3,3178)PF(1) 

TYPE  3178, PF(1) 

3178  FORMAT(6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD 8X, 1PE12 . 1) 
WRITE(3,3179)PF(2) 

TYPE  3179, PF  (  2) 

3179  FORMAT(6X, 'NORMAL  BREATHING  LEFT ' , 18X , 1PE12 .1) 

WRITE (3 , 3180) PF (3 ) 

TYPE  3180, PF<3) 

3180  FORMAT (6X, 'NORMAL  BREATHING  RIGHT ’, 17X, 1PE12 . 1) 

WRITE(3,3181)PF<4) 

TYPE  3181,  PF  (  4) 

3181  FORMAT (6X, 'NORMAL  BREATHING  DOWN' ,18X , 1PE1 2 . 1) 

WRITE(3,3182)PF(5) 

TYPE  3182, PF(  5) 

3182  FORMAT(6X, 'NORMAL  BREATHING  UP' ,2 OX , 1PE12 . 1) 

WRITE (3 , 3183 ) PF (6) 

TYPE  3183, PF  (6) 

3183  FORM AT ( 6 X, 'DEEP  BREATHING  STRAIGHT  AHEAD ' , 1  OX , 1PE1 2 . 1) 
WRITE(3,3184)PF(7) 

TYPE  3184, PF(7) 

3184  FORMAT(6X, 'DEEP  BREATHING  LEFT ' , 2 OX , 1PE1 2 . 1) 

WRITE(3,3185)PF(8) 

TYPE  3185, PF  ( 8) 

3185  FORMAT(6X, 'DEEP  BREATHING  RIGHT ', 19X , 1PE12 . 1) 

WRITE(3,3186)PF(9) 

TYPE  3186, PF(9! 

3186  FORMAT(6X, 'DEEP  BREATHING  DOWN' ,20X , 1PE12 . 1) 

WRITE(3,3187)  PF(10) 

TYPE  3 187  ,  PF  ( 1 0) 

3187  FORMAT (6X, 'DEEP  BREATHING  UP' ,22X, 1PE12.1) 

WRITE  (3 , 3 188)  PF  ( 11) 

TYPE  3188,  PF(ll) 


•** 
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3188  F0RMAT16X, ‘TALKING* ,32X, 1PE12.1) 

WRITE  ( 3 , 3189)  PF  ( 1 2) 

TYPE  3189 ,  PF  ( 12) 

3189  FORMAT(6X, 'FACIAL  GRIMACING ' ,23X , 1PE12 . 1) 

WRITE <3 ,3190)  PF  (13) 

TYPE  3190  ,PF  (13) 

3190  FORMAT(6X, ‘SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C*  (NORMAL  BREATHING) ' ,21X,1PE12.1) 

WRITE (3 , 3191 )  PF  (14) 

TYPE  3191  ,PF  (14) 

3191  FO RM AT (6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (NORMAL  BREATHING)  ' ,21X , 1PE12 .1) 

WRITE  (3, 3192)  PF  (15) 

TYPE  3192, PF  (15) 

3192  FORMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING)  '  ,23X ,  1PE12 .1) 

WRITE  (3 , 3193)  PF  ( 16) 

TYPE  3193  ,PF  (16) 

3193  FORMAT(6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C'  (DEEP  BREATHING) ' ,23X , 1PE12 .1////) 

WRITE  (3 , 3  07  0)  PF  ( IDLP) 

TYPE  3070  ,PF  ( IDLP) 

TYPE  2006 
WRITE(3 ,3176) 

TYPE  3195, WPF 
WRITE (3 ,3195)WPF 
3177  CONTINUE 
TYPE  2006 

TYPE  41 37  , ICCBS (1 )  ,  ICCBS  (2) 

TYPE  2006 

FORM  THE  ARRAYS  FOR  PLOTTING  PURPOSES  AND  LOAD  THEM  INTO  THE 
FILE  CALLED  GRPHX.XXX. 


FIND  THE  MINIMUM  AND  MAXIMUM  OF  THE  (X)  DATA  VALUES 

AND  FORM  A  DELTA  ELEMENT  SO  THAT  INTERMEDIATE  VALUES 

CAN  BE  CALCULATED  USING  THE  POLYNOMIAL  CURVE  FIT  EQUATION. 

TYPE  3061 
DO  2  1*1  ,NP 
GRPH1 ( I ) »XX AS ( I ) 

2  CONTINUE 
Nl-NP-1 
DO  8  1*1, N1 
11*1  +  1 
DO  8  J*I1 ,NP 

IF (XAS ( J) . GE , XAS ( I ) )  GO  TO  8 
XSAVE-XAS ( I ) 

XAS ( I ) *XAS( J) 

XAS ( J) *XSAVE 
8  CONTINUE 
XMIN*XAS(1> 

XMAX-XAS(NP) 

DEL* (XMAX-XMIN) /400 . 0 

XXX-XMIN-DEL 

DO  7  KA*1 ,401 
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XXX=XXX  +  DEL 
XX (KA) =XXX 
XXXX=XXX 

IF(IAX.EQ.l)  XEST*XXXX 
IF ( IAX. EQ . 2)  XEST-l.O/XXXX 
IFUAX.EQ.3'  XEST=DEXP (XXXX) 

IF ( IAX. EQ. 4)  XEST=DEXP(-XXXX) 

IF ( IAX. EQ . 5)  XEST=DEXP (1 . O/XXXX) 

IF(IAX.EQ,6)  XEST=DEXP(-1. O/XXXX) 

IF ( IAX. EQ.7)  XEST*DLOG(XXXX) 

IF ( IAX. EQ. 8)  XEST=DLOG10(XXXX) 

IF(IAX.EQ.9)  XEST=DSIN(XXXX) 

BEGIN=0 . 0 
DO  6  KB«1,M+1 

IF (KB. EQ. 1 .AND .  XEST .  EQ.  0 . 0)  GO  TO  5 
BEGIN=A (KB) *XEST** (KB-1) +BEGIN 
GO  TO  6 

5  BEGIN=A(KB) *1 . 0** (KB-1) +BEGIN 

6  CONTINUE 
YY (KA) =BEGIN 

7  CONTINUE 
DO  4  1=1,401 
GRPH2 ( I ) =XX ( I ) 

4  CONTINUE 
DO  3  1=1/401 
GRPH2  (1  +  401)  =YY  ( I ) 

3  CONTINUE 
DO  1  I— 1 ,  NP 
GRPH1(I+NP)»Y(I) 

1  CONTINUE 

OUTPUT  THE  INFORMATION  THAT  CAN  BE  USED  FOR  GENERATING  A 
CALIBRATION  PLOT  INTO  THE  FILE  CALLED  GRPHX.XXX. 

WRITE(1 .2007) 

WRITE (3  ,9) 

9  FORMAT (IX, 2 4X, 'THE  ARRAYS  FOR  PLOTTING'///) 

WRITE (1 , 2007 ) 

WRITE(1 ,3006) 

3006  FORMAT (IX, 2 4X, 'GRPHl  ARRAY  CONTAINS  THE  SODIUM  CHLORIDE 
C  CALIBRATION  STANDARD  DATA'///) 

WRITE ( 1 ,3004 ) 

3004  F0RMAT(1X,24X, 'SAMPLE  NUMBER' , 19X ,' VALUE ' /) 

DO  34  1=1 ,NP*2 

WRITEd, 3005)1, GRPHl(I) 

3005  F0RMAT(1X,28X, 13 ,22X, 1G13 .6! 

34  CONTINUE 

WRITEd, 2007) 

WRITEd, 3007) 

3007  FORMAT ( IX , 2 4X , ' GRPH2  ARRAY  CONTAINS  THE  ESTIMATED  VALUES'///) 
WRITE (1 ,3004 ) 

DO  32  1=1,802 
WRITE (1,3005)  1 ,GRPH2 ( 1 5 
32  CONTINUE 
29  TYPE  2006 

TYPE  4137, ICCBS ( 1) , ICCBS (2) 

TYPE  2006 
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TYPE  27 

27  FORMAT (IX, *DO  YOU  WISH  TO  USE  THE  SAME  SODIUM  CHLORIDE  STANDARD'/ 
C'  CALIBRATION  VOLTAGE  MEASUREMENTS  AND  EXERCISE  INTEGRATOR  COUNT1/ 
C*  DATA  BUT  CALCULATE  A  DIFFERENT  DEGREE  OR  FORM  OF  THE'/ 

C'  POLYNOMIAL  CURVE  FITTING  FUNCTION?  (ANSWER  YES  OR  NO)',2X,$> 
ACCEPT  1002, REP 
IF ( REP . EQ . YES)  GO  TO  6000 
IF ( REP . EQ . NO)  GO  TO  28 
IF(REP. NE.YES.AND.REP.NE.NO)  GO  TO  29 

28  CONTINUE 
CLOSE <UNIT=2) 

CLOSE (UNIT=3) 

CLOSE <UNIT=1) 

9674  TYPE  2  007 

TYPE  4137, ICCBS ( 1 ) , ICCBS (2) 

TYPE  2006 
TYPE  14 

14  FORMATdX, 'DO  YOU  WISH  TO  CALCUIATE  PROTECTION  FACTORS  FOR'/ 

C'  ANY  OR  ALL  OF  THE  FOLLOWING  CONDITIONS:  1)  A  DIFFERENT'/ 

C’  SUBJECT  2)  A  DIFFERENT  SET  OF  SODIUM  CHLORIDE  CALIBRATION'/ 

C'  STANDARD  VOLTAGE  MEASUREMENTS  3)  A  DIFFERENT  SET  OF'/ 

C'  EXERCISE  INTEGRATOR  COUNT  DATA.'/) 

TYPE  9999 

9999  FORMATdX,  '  (ANSWER  YES  OR  NO)',2X,$) 

ACCEPT  1002  tCOM 

IF (COM. EQ. YES)  GO  TO  5101 

IF (COM.EQ.NO)  GO  TO  9672 

IFtCOM. NE.YES.AND.COM.NE.no)  GO  TO  9674 

9672  CONTINUE 
TYPE  2  007 

TYPE  4137  ,  ICCBS ( 1 )  ,  ICCBS (2) 

TYPE  2007 
TYPE  9599 

9599  FORMATdX, 12X, 'JOB  SUCCESSFULLY  COMPLETED',//////) 

STOP 

END 


SUBROUTINE  CLEAR  IS  CALLED  WHEN  ENTERRING  DATA  IN  THE 
SODIUM  CHLORIDE  CALIBRATION  CONCENTRATION  STANDARDS 
SECTION  AND  EXERCISE  INTEGRATOR  COUNT  SECTION.  BY 
CALLING  THIS  SUBROUTINE  THE  OPERATOR  CAN  MAKE  TWO 
DATA  ENTRIES  ON  THE  SAME  LINE.  SUBROUTINE  CLEAR  ERASES 
THE  LINE  INWHICH  THE  FIRST  DATA  ENTRY  WAS  MADE  AND 
RETYPES  THAT  LINE,  INCLUDING  THE  FIRST  DATA  ENTRY; 

THIS  ALLOWS  A  SECOND  ENTRY  TO  BE  MADE  ON  THAT  LINE. 


SUBROUTINE  CLEAR! LINES) 
BYTE  A ( 3 > 

A(l)  -27 
A (2) -65 
A(3) -75 

IF  (LINES. EQ.O)  LINES-1 
DO  1  1-1, LINES 
1  TYPE  4 , A ( i ) ,A(2) ,A(1) ,A(3) 
LINES-Q 
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RETURN 

4  FORMAT  <1H+,4A1,$) 
END 


APPENDIX  B: 


DATA. XXX  File  Contents  for  Data  in  Table  7 
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tiie  number  r»F  sodium  chloride  calibration  concentration  standards*  i 


SAMPLE  NUMBER 


VOLTAGE 
MEASUREMENT 
(Itl  VOLTS ) 

(  X  DATA  ) 


SODIUM  CIIIORIOE 
CALIBRATION 
CONCENTRATION 
(  Y  DATA  ) 


1 

2 

3 

4 

5 

6 
7 


3 ,30000 
2.91500 
2.31000 
1.50000 
0.545000 
0.165000 
0.105000 


1.00000 
0.100000 
0 . 1 OOOOOE-Ol 
0 . 1 00000E-02 
0 . 1 OOOOOE-03 
0 . 1 00000E-04 
0. 100000E-05 


SUBJECT  NAME )  CAPTAIN  EDWARD  S.  tCOLKSAR,  JR. 
TYPE  Of  NARKS  USA  s  M17  -  MEDIUM  -  MO  (.MASSES 
DATE  TESTED  t  9  APR  1 1,  1980 
TIME  TESTED  *  1330  HOURS 


EXE  IK*  tSE  INTEGRATOR  CODIN'  DATA: 


:  EE |K*  IRE 

INTEGRATOR  COtllM* 

TIME  PERIOD 
UN  SECONDS) 

M‘ 'RMAI.  r.REATIUNG  STRAIGHT  AHEAD 

3904 

10 

Dt:t:p  BREATHING  STRAIGHT  AHEAD 

47  51 

10 

TALKING 

S1DE-TO-SIDE  HEAD  MOVEMENTS 

4620 

10 

tuEEp  BREATHING) 

"P- AND-OOUN  HEAD  MOVEMENTS 

3976 

10 

tDRi.p  BREATH  lilt!) 

4016 

10 

!'A‘*t At.  GRIMACING 

4937 

10 

Hm  bummkv  num 


r*r 


\ 

» 

i 


/ 


s  r:mmna  tuM  blok-ma  nuu> 
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Tin:  ouni'i'.  or  tiu:  desired  1‘oi.yiioiit.m.=  c 

the  roi-YikaiTAf,  ruiKTiotiAf.  iutjijitioh  of  the  vap.i.w'm:  oo  is  in  terms  op:  ( expoueutial  tx>  > 


cni:rm:t!:iiT  nunnur. 

ro!;mcn::iT  nuincr. 
coefficient  tnmnEU 

I'onmcimiT  imiinER 
cormcirm’  uofinEi: 
roiTFicir.tiT  tiunitEi: 
mr.Fncir.tr  tutunr.r. 


1=  -4 .  1790  “IE-05 
?.=  -7.31407E-0S 

7=  1.25P.C3E-04 

4=  -2.6fi095K-05 

5=  3.«0I3fiK-0G 

6=  -1.7IS4U-.-07 

7=  3.fl3'!C6c-n9 


tm:  re.*:  i  oiiAi.0  ,m:i:  c.M.ctKA-ir.n  r.y  tiu:  folumiuk; 

EOltjV-'IOtl:  1 Y  (DATA)  -Y  (  ESV 1 1  SAT  ID)  ) 


✓ 


■i.r  siu:i.':i:i: 


Ytl.V'-A 


y (estimated) 


RESIDUAL 


t 

i 


1 .ooooo 
(1.)  00000 
o.tooooor-oi 
0,  looocfir.-02 
o.iouotuif:-oj 
o.  iaot)0or-o.i 
O.IOOOOOE-OS 


1.00000 
o.  100000 
o.iaooooc-oi 
0.1000006-02 
O.1OO0Q4K-D3 
0 , 1 000346-04 
0.1004136-05 


-o.40354er:-07 
-o.ioomc-o? 
-0.71 1342E-00 
-0.4019296-01: 
-0.390294E-0r 
-0.3406026-01*. 
-0.4325026-0E 


nu.*:  ritr.  sou/.w:  residual::*  ?,ooiuifc-l3 


pwrffffNO  mm  M  fur-hit  nuts 
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THi.  nnrcni  -’-tvr  and  protection  pactor.  calcuiatione: 


A —  pi' ION  pr.C’OR  'M! ,  .p 

’.o::+0'.  !■;.;•  ree"  aeoigned  this  value  dy  orr 

i  "v"~  cir.’oj  *pj\ij*yy  Qr  pnrm  t  '*  f!”  *1!!"  *. 


*•*'  dap.t  in  ten  to  the  '“ix’*’;* .  t.e  i n~’ lcnn\t r. 

con”?  value  rcr.  ••.  parti  cl-:'  rc  i~  f.  i  ::  cunmc:: 

MCRCLY  REPRESENTATIVE  Of  INTEGRATING  T!!E  ELCCTP.I CAL 

r t:  apt  the  true  ppotdctioi.  .'actor  is  indee?  lees 

THAN  1.0E+(1f..  ANY  EXEP.CISE  SCALED  INTEGRATOR 
COUITT  VALUE  YIELDING  A  PROTECTION  FACTOR  GREATER 
THAN  1.0  '..’ILL  EE  REPORTED  AS  0.0E-01. 


SUBJECT  NAME:  CAPTAIN  EDWARD  S.  KOLESAR,  JR. 
TYPE  OP  MASK:  USA  :  Ml 7  -  MEDIUM  -  NO  GLASSES 
DATE  TESTED  :  3  APRIL  1980 

TIME  TESTED  :  1330  HOURS 


EXERCISE  PROTECTION  EACTOP. 

NORMAL  E-REATHING  STRAIGHT  AHEAD  1.8E*0< 

DEEP  BREATHING  STRAIGHT  AHEAD  1.3E  +  0C 

TALKING  1.3E+0? 

SIDE-T^-SIDE  HEAD  MOVEMENTS 

(DEES'  RRSATtUNC)  1  .?£♦<>.* 

t’P-AMP-DOUN  !(EAD  MOVEMENTS 

(DEEP  BREATHING)  I.7*>3< 

FACIAL  GRIMACING  1.2E*0i 


OVERALL  ARITHMETIC  AVERAGE  PROTECTION  FACTOR 

FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMS?  »  1 .SE+04 


OVERALL,  TIME  WEIGHTED  AVERAGE  PROTECTION  FACTOR 
FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMS) 


l.S£»04 


APPENDIX  0; 

GRPHX.XXX  File  Contents  for  Use  with  NACLGRAPH .FTN  Program 
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TRE  ARRAYS  FOR  PLOTTING 


GPPiil  ARRAY  CONTAINS  THE  SODIOtl  CHLORIDE  CALIBRATION  STANDARD  DATA 


S.V'.PLF  ;x"!FRF 


VALUE 


1 


e 

7 

3 

9 

10 

n 

n 

i  :> 


5.1P.C0C 

2.91500 

c.nooo 

1.50000 
0.545000 
0.165000 
0 . 105000 
1.00000 
0.100000 
O.IOOQOOS-Ol 
0 . 1C0C0QE-Q2 
0 , 100000E-03 
0 . 100000E-04 
0.100000E-05 


GRPH2  ARRAY  CONTAINS  THF.  ESTIMATED  VALUES 


SAMPLE  NUMBER 


VALUE 


1 

I 

3 

4 

5 

6 
7 
(1 
9 

10 


0.105000 
0.11318? 
0.121375 
0.129562 
0.137750 
0.145930 
0 . 1541 25 
0.162313 
0.170500 
0.178607 


f 

«1 


momm  jusi umumwujm 
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11 

0 .  186C75 

12 

0.195063 

13 

0.203250 

14 

0.211438 

15 

0.219625 

16 

0.227812 

17 

0.236000 

16 

0.244188 

19 

0.252375 

20 

0.260563 

21 

0.268750 

22 

0.276930 

23 

C. 235125 

24 

0.293313 

25 

0.301500 

26 

0.309687 

27 

0.317875 

28 

0.326063 

29 

0.334250 

30 

0.342438 

31 

0.350625 

32 

0.358813 

33 

0.367000 

34 

0.375188 

35 

0.383375 

36 

0.391563 

37 

0.399750 

38 

0.407937 

39 

0.416125 

40 

0.424313 

41 

0.432500 

42 

0.440688 

43 

0.448875 

44 

0.457063 

45 

0.465250 

46 

0.473438 

47 

0.481625 

48 

0.489813 

49 

0.498000 

50 

0.506187 

51 

0.514375 

52 

0.522563 

53 

0.530750 

54 

0.538938 

55 

0.547125 

56 

0.555313 

57 

0.563500 

58 

0.571688 

59 

0.579875 

60 

0.588063 

61 

0.596250 

62 

0.604438 

63 

0.612625 

64 

0.620813 

65 

0.629000 

66 

0.637188 

67 

0.645375 

68 

0.653561 

69 

0.661750 

70 

0.669937 

71 

0.678125 

102 
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72 

72 

74 

75 

76 

77 

78 

79 

80 
31 
82 

83 

84 

85 

86 
87 

8e 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 


0.686312 
0 .G9450C 
0.702688 
0.710875 
0.719063 
0.727250 
0.735438 
0.743625 
0.753813 
0.760000 
0.768188 
0.776375 
0.784563 
0.792750 
0.800938 
0.809125 
0.817313 
0.825500 
0.833688 
0.841875 
0.850063 
0.858250 
0.866438 
0.874625 
0.882812 
0.891000 
0.899188 
0.907375 
0.915563 
0.923750 
0.931938 
0.940125 
0.948313 
0.956500 
0.964688 
0.972875 
0.981063 
0.989250 
0.997438 
1.00563 
1.01381 
1.02200 
1.03019 
1.03838 
1.04656 
1.05475 
1.06294 
1.07113 
1.07931 
1.08750 
1.09569 
1.10388 
1.11206 
1.12025 
1.12844 
1.13663 
1.14481 
1.15300 
1.16119 
1.16938 
1.17755 
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133 

1.18575 

134 

1.19394 

135 

1.20213 

136 

1.21031 

137 

1.21850 

133 

1.22669 

139 

1.23488 

140 

1.24306 

141 

1.25125 

142 

1.25944 

143 

1.26763 

144 

1.27581 

145 

1.28400 

146 

1.29219 

147 

1.30037 

148 

1.30856 

149 

1.31675 

150 

1.32494 

151 

1.33312 

152 

1.34131 

153 

1.34950 

154 

1.35769 

155 

1.36588 

156 

1.37406 

157 

1.38225 

158 

1.39044 

159 

1.39863 

160 

1.40681 

161 

1.41500 

162 

1.42319 

163 

1.43138 

164 

1.43956 

165 

1.44775 

166 

1.45594 

167 

1.46413 

168 

1.47231 

169 

1.48050 

170 

1.48869 

171 

1.49688 

172 

1.50506 

173 

1.51325 

174 

1.52144 

175 

1.52963 

176 

1.53781 

177 

1.54600 

178 

1.55419 

179 

1.56238 

180 

1.57056 

181 

1.57875 

182 

1.58694 

183 

1.59513 

184 

1.60331 

185 

1.61150 

186 

1.61969 

187 

1.62788 

188 

1.63606 

189 

1.64425 

190 

1.65244 

191 

1.66063 

192 

1.66881 

193 

1.67700 
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194 


1.68519 


195 

196 

197 
190 

199 

200 
201 
202 

203 

204 

205 

206 
207 
200 

209 

210 
211 
212 


1.69330 

1.70156 

1.70975 

1.71794 

1.72613 

1.73431 

1.74250 

1.75069 

1.75888 

1.76706 

1.77525 

1.78344 

1.79163 

1.79981 

1.80800 

1.81619 

1.82438 

1.83256 


213 

214 

215 

216 
217 
21S 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
23* 

239 

240 

241 

242 

243 

244 

245 

246 

247 
241 

249 

250 

251 

252 

253 

254 


1.84075 
1.84894 
1.85713 
1.36531 
1,87350 
1.88169 
i.88988 
1.89806 
1.90625 
1.91444 
1.92263 
1.93081 
1.93900 
1.94719 
1 .95538 
1.96356 
1.97175 
1,97994 
1.98813 
1.99631 
2.00450 
2.01269 
2.02087 
2.02906 
2.03725 
2.04S44 
2.05363 
2.06181 
2.07000 
2.07919 
2.0863? 
2.09456 
2.10275 
2.11094 
2.11913 
2.12731 
2.13SSO 
2.14169 
2. 15188 
2.16006 
2.16825 
2.176*4 
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:s5 
2  55 
25“ 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
232 

283 

284 

285 

286 
287 
28B 
209 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 
3  OS 

306 

307 

308 

309 

310 

311 

312 


3  IS 


2.18463 

2.19281 

2.20100 

2.20919 

2.21738 


2.22556 

2.23375 

2.24194 

2.25013 

2.25831 

2.26650 

2.27469 

2.28288 

2.29106 

2.29925 

2.30744 

2.31563 

2.32381 

2.33200 

2.34019 

2.34338 

2.35656 

2.36475 

2.37294 

2.38112 

2.38931 

2.39750 

2.40569 

2.41388 

2.42206 

2.43025 

2.43844 

2.44662 

2.45481 

2.46300 

2.47119 

2.47930 

2.48756 

2.49575 

2.50394 

2.51213 

2.S2031 

2.52850 

2.53669 

2.54488 

2.55306 

2.56125 

2.56944 

2,57763 

2.58S61 

2.59400 

2.60219 

2.61038 

2.61856 

2.62675 

2.63494 

2.64313 

2.65131 

2.65950 

2.66769 

2.67580 
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316 

2.68406 

317 

2.69225 

318 

2.70044 

319 

2.70863 

320 

2.71681 

321 

2.72500 

322 

2.73319 

323 

2.74138 

324 

2.74956 

325 

2.75775 

326 

2.76594 

327 

2.77413 

328 

2.78231 

329 

2.79050 

330 

2.79869 

333 

2.80687 

332 

2.81506 

333 

2.82325 

3  34 

2.83144 

335 

2.83963 

336 

2.84781 

337 

2.85600 

338 

2.86419 

3  39 

2.87238 

340 

2.880  56 

341 

2.88675 

342 

2.89694 

343 

2.90513 

344 

2.91331 

345 

2.92150 

346 

2.92969 

347 

2.93768 

348 

2.94606 

349 

2.95425 

350 

2.96244 

351 

2.97063 

352 

2.97881 

353 

2.98700 

354 

2.99519 

355 

3.00336 

356 

3.01156 

357 

3,01975 

356 

3.02794 

359 

3.03613 

360 

3.04431 

361 

3.05250 

362 

3.06069 

363 

3.06866 

364 

3.07706 

365 

3.06525 

366 

3.09344 

367 

3.10163 

366 

3.10981 

369 

3.11800 

370 

3.12619 

371 

3.13438 

372 

3.14256 

373 

3.15075 

374 

3.15894 

375 

3.16713 

376 

3.17531 
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377 


3.18350 


378  3.19162 

379  3.19988 


380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 


3.20806 

3.21625 

3.22444 

3.23263 

3.24081 

3.24900 

3.25719 

3.26538 

3.27356 

3.28175 

3.28994 

3.29813 

3.30631 

3.31450 

3.32269 

3.33088 


396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

4"7 

W 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 


3.33906 
3.34725 
3.35544 
3.36363 
3.37181 
3.38000 
0.100432E-05 
0.216730E-05 
0.3350336-05 
0.455367E-05 
0. 5777636-05 
0.702248E-05 
0.8288526-05 
0.9576055-Q5 
0.10G854E-04 
0.122160E-04 
0.135706E-04 
0 .14947  IE-04 
0.163466E-04 
0.177694E-04 
0.192160E-04 
0.206865E-04 
0.2218136*04 
0.237007E-04 
0.2524526*04 
0.2681506-04 
0.2841046*04 
0.3003196*04 
0.3167976-04 
0.3335436-04 
0.3505606*04 
0.3678516*04 
0.3854216-04 
0.4032736-04 
0.4214116-04 
0.4398396-04 
0.4S85606-04 
0.4775796-04 
0.4969006-04 
0. 516527 E-04 
0.5364636*04 
0.556714E-04 
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433 
4  3? 

4  40 

441 

442 
4  43 

444 

445 

446 

447 

448 

449 
4  50 

451 

452 

453 

454 

455 
4  56 
4  57 
4  58 

459 

460 

461 

462 

463 

464 

465 

466 

467 
4f  8 
469 
<70 

471 

472 

473 

474 

475 
4.6 

477 

478 

479 

48v> 

4tl 

462 

46j 

464 

465 

466 

467 
403 
469 

490 

491 

492 

493 

494 

495 

496 

497 
4*C 


0 . 57  7283E-G4 

0 . 5981 74E-04 

0.619392E-04 

0 . 64  09  42E-04 

0.662827E-04 

0 .685053  E-04 

0.707623E-04 

0.730542E-04 

G.753816E-04 

0.777449E-04 

0,801445E-04 

0.825810E-04 

0.850548E-04 

0 .875665E-04 

0.901165E-04 

0.927055E-04 

0  .953  338E-04 

0.980021E-04 

0 . 1007 1  IE-03 

C.103461E-03 

0 . 1 06252E-03 

0.1 Q9Q86E-Q3 

0.111962E-03 

0.114881S-03 

Q.117845E-03 

0.120C53E-03 

Q.1239G6E-Q3 

Q.127005K-03 

1.1301506-03 

0.133342E-Q3 

0.136582E-03 

Q.139B70E-03 

0,1432076-03 

0.1465946-01 

0.150031E-03 

0.1S3519E-03 

0.1570596-03 

0.16QG92E-Q3 

0,1642936-03 

0 .1679986-03 

0.1717536*03 

0.1755446*63 

0.179430E-03 

0.1633556-03 

0.1673376-03 

0.1913766-03 

0.19  54606-03 

0.1996416-03 

0. 2036656-03 

0,2061516-03 

0.2 1250  l.E-03 

0.2169156-03 

0.2213956-03 

0.2259416-03 

0. 2305556-03 

0. 135 2376-03 

0.2399696-03 

0./44I’ 2E-03 

0,2497066-03 

0.2546746-03 

0.2597166-03 
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499 

500 

501 

502 

503 

504 

505 

506 

507 
SOS 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 
S  20 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 
515 
536 
517 

538 

539 

540 

541 

542 
S4J 

544 

545 

546 

547 

548 

549 
5  SO 

551 

552 

553 

554 

555 

556 

557 

558 

559 


0 . 264834E-03 
0 . 270028E-03 
0.275300E-Q3 
0.280652E-03 
0.286084E-03 
0.291598E-03 
0.297196E-03 
0.302878B-03 
0.308647E-03 
0.314504E-03 
0.320450E-03 
0.326488E-03 
0.332617E-03 
0.338842E-03 
0.345162E-U3 
0.351580E-03 
0.358098E-03 
0. 364718E-Q3 
0.371441E-03 
Q.37826SE-03 
O.3652Q5E-03 
0 . 392251E-Q3 
0.399409E-03 


0.406680E-03 

0 . 414068E-03 

0.42157SE-03 

0,4392036-03 

0.436954E-03 

0. 44483  16-03 

0.452838E-03 

0,4609766-03 

0,4692466-01 

0,47 765 86-Q3 

0.486308E-03 

0.  **49016-03 

C.S01741E-03 

0,5127316-03 

0.5218746-03 

0.531174E-03 

0.5406146-03 

0.5502586-03 

0,5600496*03 

0.5700136-03 

O.S60 1516-03 

0.5904726-03 

0. 6009  7  66-03 

0.611C69E-03 

0.623SSSE-03 

0.633640E-03 

0 ,64  49286-03 

0. 6564246-03 

0.660) 34E-03 

0.6800636-03 

0,6922176-03 

0.7046016-03 

0.7172216-03 

0.7300646-03 

0.7431966-03 


0.7545636-03 

0.7701936-03 


0. 7640926-01 
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560 

561 

562 
56  3 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 
509 

590 

591 

592 
591 

594 

595 
594 

597 

598 

599 

600 
601 
602 

603 

604 

605 
C0« 
607 
601 

609 

610 
611 
612 
623 

614 

615 

616 
617 
611 

619 

620 


0.798268E-03 
0.81 2729E-03 

0 .8274  81 E-03 
0.842533E-03 
0.857894E-03 
U.873571E-03 
0 . 8  89574E-03 
0.905912E-03 
0.922595E-03 
0.939631E-03 
0.957032E-03 
0.974808E-O3 
0.992969E-03 
0.1 01153E-02 
0.103049E-02 
0.104988B-02 
0.106970E-02 
0.108996E-02 
0.1110686-02 
0.113187E-02 
0.1153SSE-02 
Q.117573E-Q2 
0.119843E-Q2 
0.122165E-02 
0.124543E-02 
0.126976E-02 
0.129468E-02 
0.132019E-02 
0.134633E-02 
0.1373096-02 
0.140051E-02 
0. 14286  IE -92 
0.145740E-02 
0.148692E-02 
0 . 1517175-02 
0,1548186-02 
0, 157998E-02 
0 . 1612596*02 
0.1646046-02 
0 . 1680356-02 
0.1T1554E-02 
0.1751666-02 
0.178*736-02 
0 . 1628766-02 
0.1865816-02 
0.1905906-02 
0.1947066-02 
0.1919326-02 
0.2032736-02 
0.2077316-02 
0.2123106-02 
0.2170156-02 
0.2211506-02 
0.2368176-02 
0.2319236-02 
0.2371706-02 
0,2425656-02 
0.2411106-02 


0.2538126-02 


Q.2596 756-02 


0.2657056-02 
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621 

622 

623 

624 

625 

626 

627 

628 

629 

630 
6  31 

632 

633 

634 
6  35 

636 

637 

638 
6  39 

640 

641 

642 

643 

644 
64  5 

646 

647 

648 

649 

650 

651 

652 
652 
6*4 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 
667 
661 

669 

670 

671 

672 

673 

674 

675 

676 

677 
676 
679 
610 
6#S 


0 .27 1906E-02 

0.27  8284E-02 

0.284846E-02 

0.2915D6E-02 

0.298541E-02 

0.305687E-02 

0.313041E-02 

0.320608E-02 

0.328398E-02 

0.336415E-02 

0.344668E-02 

0.353164E-92 

0.361911B-02 

Q.370918E-02 

0.380 19 2E-Q2 

0.3B9742S-02 

0.399577E-02 

0.409706E-02 

0.420139E-03 

0 .43QG85E-C 1 

0.441955E-02 

0 .4  S3  3S9E-02 

0.465108E-02 

0,47 7JX3E-02 

0.48968SE-02 

O.S02S37E-02 

0.515780E-02 

0.S29427E-02 

0 .S43492E-02 

0.557987E-0I 

0 .5729276-02 

0 .5883266-02 

0.6041995-02 

0.620562E-03 

0.6374295-02 

0. 6548186-02 

0.6727456-02 

0.6912295-02 


0.7102666-02 
0.72993SE-02 
0.7S0197E-02 
0.7710906-02 
0.7926365-32 
0. 6148556-02 
0.6377695-02 
0.1614025-02 
0.6657765-02 
0.9109165-02 
0.9356415-02 
0.96 35965-02 
0.9911666-02 
0.1016656-01 
0.1069016-01 
0.1079315-01 
0. 1110565-01 
0.1142615-01 
0.1176065-01 
0.1216416-01 


0.1243646-01 

0.1262396-01 


0. 1320126-01 
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682 

0 . 135905E-01 

6S3 

0 . 139923E-01 

684 

0 . 14  4070E-01 

685 

0 . 14  8351E-01 

686 

0 . 1527  69E-01 

687 

0 . 157330E-01 

688 

0.162038E-01 

689 

0 . 166899E-01 

690 

0 . 171917E-G1 

691 

0.177099E-01 

692 

0 . 182449E-01 

693 

0.187974E-O-1 

694 

0 . 193679E-01 

695 

0 . 199570E-01 

696 

0.205656E-01 

697 

0 .2119  41E-01 

698 

0.218433E-01 

699 

0 . 225140F-01 

700 

0.232068E-01 

701 

0 . 239227E-  01 

702 

0 . 246  623E-01 

703 

0 . 254267E-01 

704 

0.262165E-01 

705 

0.270329E-01 

706 

0 . 2787  67E-01 

707 

0.287489E-01 

708 

0.29S506E-01 

709 

0.305829E-01 

710 

0.315469E-01 

711 

0.325438E-01 

712 

0.335748E-01 

713 

0.346412E-01 

714 

0.357443E-01 

715 

0.368856E-01 

716 

0.380665E-01 

717 

0.392885E-01 

719 

0.405533E-01 

719 

0 . 418624E-01 

720 

0.432177E-01 

721 

0.446209E-01 

722 

0.460739E-01 

723 

0.475787E-01 

724 

0.491375E-01 

725 

0.507523E-01 

726 

0.524255E-01 

727 

0.541594E-01 

728 

0.559565E-01 

729 

0.578195E-01 

730 

0.597510E-01 

731 

0.617540E-01 

732 

0.638314E-01 

733 

0.659864E-01 

734 

0.682223E-01 

735 

0.705426E-01 

736 

0.729509E-01 

737 

0.754510E-01 

738 

0.780469E-01 

739 

0.807429E-01 

740 

0.835433E-01 

741 

0 . 864528E-01 

742 

0.894763E-01 
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743 

0 .926 189E-01 

744 

0 . 9  58  860E-01 

745 

0 . 9  92332E-01 

746 

0.102816 

747 

0.106492 

•748 

0.110317 

749 

0.114297 

750 

0.118441 

751 

0.122755 

752 

0.127248 

753 

0.131929 

754 

0.136806 

755 

0.141889 

756 

0.147187 

757 

0.152711 

758 

0.158473 

759 

0.164483 

760 

0.170755 

761 

0.177300 

7  62 

0.184133 

763 

0.191267 

7  64 

0.198718 

765 

0.206502 

766 

0.214636 

767 

0.223137 

768 

0.232024 

769 

0.241316 

770 

0.251035 

771 

0.261202 

772 

0.271840 

773 

0.282975 

774 

0.294631 

775 

0.306835 

77G 

0.319618 

777 

0.333007 

778 

0.347037 

779 

0.361740 

780 

0.377152 

781 

0.393310 

782 

0.410255 

783 

0.428028 

784 

0.446673 

785 

0.466237 

786 

0.486770 

787 

0.508323 

788 

0.530952 

789 

0.554714 

790 

0.579671 

791 

0.605887 

792 

0.633431 

793 

0.662376 

794 

0.692797 

795 

0.724775 

796 

0.758395 

797 

0.793747 

798 

0.830926 

799 

0.870032 

800 

0.911171 

801 

0.954454 

802 

1.00000 
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DIMENSION  XYK50)  ,XY2(810)  ,TTL4(5)  ,NOE(2) ,IPTTL(4> 

DIMENSION  IPCRV(2) ,KTYPE(2) ,TTL3(4> ,OTTL(4) ,TSIZE(4) , RYTTL ( 4 ) 
BYTE  XTL ( 40)  ,YTL(40)  ,TTL(240)  ,DNAME(130)  ,PNAME(13) 

BYTE  YES ,NO,ANS,SYMB 

BYT E  RED , GR EEN , BLACK , PSCL , PTL , PCRV ,  ESC , HOME , CLR 
COMMON/GRPH/YTLRX , XTLRY 
DATA  YES/ 'Y' / 

DATA  NO/'N'/ 

DATA  RED/ 'R'/ 

DATA  GREEN/ 'G'/ 

DATA  BLACK/ 'L  / 

DATA  TTL3/;«  AC','TUAL','  DAT' , 'A'/ 

DATA  TiTi4/  '«  ES  ',  'TIMA' ,  'TED  ' ,  'VALU* ,  'ES  '/ 

ESC-27 
HOME- 7 2 
CLR*74 
SYMD* 1 
ISCL-0 

CALL  ERRSET(64, . /RUE. , .TRUE. , .TRUE. , .FALSE. ,3 0) 

CALL  ERRSET(29,.  TRUE.,  .TRUE.  ,  .TRUE. ,  .FALSE.  ,30) 

CALL  ERRSET(43,.TR0E., .TRUE., .FALSE., .FALSE. ,30) 

CALL  ERRSET (59,. TRUE., .FALSE., .TRUL., .FALSE., 30) 

CALL  ERRSET (3 0> .TRUE. , . FALSE. , .TRUE. , .FALSE. , 30) 

11  CONTINUE 

WRITE(5f 555)  ESC, HOME, ESC, CLR 
WRITE(5, 415) 

1  READ (5 , 405 ,ERR-2)  IRW 

IF ( IRW . EQ . 1 , OR. IRW . EQ . 2 . OR . IRH . EQ , 3 )  GO  *0  3 

2  WRITE(5,155) 

GO  TO  1 

3  IF (IRW. EQ. 2. OR. IRW. EQ. 3)  IRLUN-5 
IF(IRW.EO.l)  IRLON-1 

WRITE  C 5  65) 

10  f READ (5, 55)  LENP , (PNAME(I) ,1-1,12) 

IF (LENP.GT.12)  GO  TO  20 
PNAME (LENP+1) -0 
CALL  ERRTST(43,K) 

CALL  ASSIGN (2, PNAME) 

CALL  ASSIGNU,  'GHAPH1.TXT') 

IF(IRW.NE.3)  CALL  ASSIGN (1 , 'GRAPH1 .COM ' ) 

CALL  ERRTST(43,K) 

IF(K.EQ,2)  GO  TO  21 

20  WRITE(5,115) 

GO  TO  10 

21  WRTTPfi  SPl 

WRI1’E(4,305)  (PNAME  (I)  ,1-1, LENP) 

30  ISCL-ISCL+1 

IXMM-0 
IYMM-0 

WRITE(5,555)  ESC, HOME, ESC, CLR 
NDIG-2 

IF(ISCL.LT.IO)  NDIG-1 

IF(IRW.EQ.2.0R.IRW.EQ.3>  WRITE<5,355)  ISCL 
IF(IRW.EQ.2.0R.IRW,EQ,3)  WRITE(5,465) 

40  REAL (IRLUN, 105)  LEN, ANS 
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IF(LEN.EQ.O)  ANS-'N' 

IF ( ANS.  EQ. YES. OR. ANS. EQ. NO)  GO  TO  50 
WRITE (5, 15 5) 

GO  TO  40 

50  IF(ANS.EQ.YES)  IXMM-1 

IF(IXMM.EQ.O)  GO  TO  110 
IF < IRW . EQ . 2 . OR.  IRW . EQ . 3)  WRITEd,  485) 
60  READdRLUN,* , ERR*70)  XMIN 

GO  TO  80 

70  WRITEd, 155) 

GO  TO  60 


80 

90 


GO  TO  110 


IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITEd,  495) 
READ  < IRLUN, * , ERR*10  0)  XMAX 


100 


WRITEd,  155) 


GO  TO  90 

110  WRITEd,  245) 

WRITE(4, 355)  ISCL 
WRITEd,  465) 

WRITEd, 265)  ANS 
IF(IRW.EQ.2)  WRITEd  ,585)  ANS 
IF (ANS. EQ. NO)  GO  TO  111 
WRITE (4 ,485) 

WRITE(4,385)  XMIN 
IF(IRW.EQ.2)  WRITEd ,435)  XMIN 
WRITE(4, 495) 

WRITE(4,385)  XMAX 
IF(IRW.EQ.2)  WRITEd, 435)  XMAX 

111  IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,475) 

120  READ (IRLUN,  105)  LEN, ANS 

IF (LEN ,EQ. 0)  ANS«'N 1 

IF(ANS. EQ.YES.OR.ANS.EQ.NO)  GO  TO  130 

WRITEd,  155) 

GO  TO  120 

130  IF(ANS.EQ.YES)  IYMM-1 

IF(IYMM.EQ.O)  GO  TO  170 
IFURW.EQ.2.0R.IRW.EQ.3)  WRITE(5,485) 

140  READ (IRLUN,*, ERR* 150)  ,YMIN 

GO  TO  160 

150  WRITE(5 , 155) 

GO  TO  140 

160  IF ( IRW . EQ . 2 . OR. IRW . EQ . 3)  WRITEd, 495) 

READdRLUN**)  YMAX 

170  WRITE(4, 475) 

WRITE(4, 265)  ANS 
IF(IRW.BQ.2)  WRITEd, 585!  ANS 
IF(ANS.EQ.NO)  GO  TO  171 
WRITE(4,485) 

WRITE(4, 385)  YMIN 
IF(IRW,EQ.2)  WRITEd, 435)  YMIN 
WRITE(4, 495) 

WRITEd, 385)  YMAX 

IF  (IRW. EQ, 2)  WRITEd, 435)  YMAX 

171  IF ( IRW , EQ . 2 . OR • IRW . EQ . 3)  WRITE(5,555) 

IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITE(5,125) 


ESC, HOME, ESC, CLR 
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350  IF(ANS.EQ.YES)  ISCX-1 

WRITE (4/245) 

WRITE(4,355)  ISCL 
VmiTE(4,255) 

WRITE<4,195) 

WRITE (4, 275)  XLEN 

IF ( IRW. EQ. 2)  WRITE(1,425)  XLEN 

WRITE (4 , 205) 

WRITE(4, 275)‘  XBLINT 
IF(IRW.EQ.2)  WRITEtl ,425)  XBLINT 
WRITE(4 , 185) 

WRITE(4 , 265 )  ANS 
IF(IRW.EQ.2>  WRITE(1,585)  ANS 

IF ( IRW . EQ . 2 . OR . IRW . EQ . 3 )  WRITE(5,555)  ESC , HOME , ESC , CLR 
IF<IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,355)  ISCL 
IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITE(5,375> 

IF ( IRW . EQ . 2 . OR . IRW . EQ . 3)  WRITE(5,215) 

360  READ ( IRLUN, * , ERR“37  0)  YLEN 

IF(YLEN.GE.1..AND.YLEN.LE.10.)  GO  TO  380 
370  WRITE (5, 155) 

GO  TO  360 

380  IF ( IRW . EQ . 2 . OR, IRW . EQ . 3)  WRITE(5,205) 

390  READ < IRLUN, *,ERR«400>  YBLINT 

IF(YBLINT.LE.O. .OR.YBLINT.GT.YLEN)  GO  TO  400 
GO  TO  410 

400  WRITE(5,155) 

GO  TO  390 

410  IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,185) 

420  READ (IRLUN, 105)  LEN, ANS 

IF(LEN.EQ.O)  ANS-'N1 
IF(ANS. BQ.YES.OR.ANS.EQ.NO)  GO  TO  430 
WRITE(5 , 155) 

GO  TO  420 

430  IF(ANS.EQ.YES)  ISCY»1 

WRITE(4,245) 

WRITE(4 , 355)  ISCL 
WRITE(4,375) 

WRITE (4, 2 15) 

WRITE (4, 27 5)  YLEN 
IF(IRW,EQ.2)  WRITE(1,425)  YLEN 
WRITE (4, 2 05) 

WRITE(4,275)  YBLINT 
IF(IRW,EQ,2)  WRITE<1 ,425)  YBLINT 
WRITE(4, 185) 

WRITE<4 , 265 )  ANS 
IF(IRW.EQ.2)  WRITEtl , 585)  ANS 

IF ( IRW . EQ . 2 .OR. IRW . EQ . 3)  WRITE(5,555)  ESC, HOME, BSC, CLR 
IF ( IRW . EQ . 2 . OR. IRW , EQ . 3 )  WRITE(5,355)  ISCL 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,95) 

IF ( IRW . EQ . 2 • OR. IRW . EQ . 3 )  WRITE(5,225) 

440  READ ( IRLUN, 4 55, ERR«45 0)  LEN, IPSCL 

IF(LEN.EQ.O)  IPSCL-1 

IFdPSCL.GT.O. AND, IPSCL.LT. 4)  GO  TO  460 
450  WRZTE(5,155) 

GO  TO  440 


} 


m  ■:  •  ■  • 
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IF(IRW.EQ.2.0R.IRW.EQ.3>  SJRITE(5,365) 

180  READ ( IRLUN, 4  05, ERR*  19  0)  IC 

IF(IC.EQ.l .0R.IC.EQ.2.0R.IC.EQ.3)  GO  TO  200 
190  WRITE(5 ,155) 

GO  TO  180 

200  IPUC.NE.l)  GO  TO  210 

NTYPE-0 
GO  TO  250 

210  IF ( IRW . EQ . 2 .OR. IRW , EQ , 3 )  WRITE(5,395) 

220  READ (IRLUN,4  05,  ERR* 23  0)  NCYCLE 

IF ( NCYCLE. GE.l)  GO  TO  240 
230  WRITE(5,155) 

GO  TO  220 

240  IF(IC.EQ.2>  NTYPE-NCYCLE 

IF(IC.EQ.3)  STY  PE— NCYCLE 

250  WRITE(4,245) 

WRITE(4,125) 

WRITE (4, 3 65) 

WRITE(4,285)  IC 
IF(IRW.EQ.2)  WRITE(1,405)  IC 
IFdC.EQ.l)  GO  TO  251 
WRITE<4,395) 

WRITE(4, 285)  NCYCLE 
IF(IRW.EQ.2)  WRITE(1,405)  NCYCLE 

251  IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,165) 

260  -  READ (IRLUN, 105)  LEN , AN S 

IF(LEN.EQ.O)  ANS-'N' 

IF(ANS. EQ.YES.OR.ANS.EQ.NO)  GO  TO  270 
WRITE(5,155) 

GO  TO  260 

270  IF(ANS.EQ.YES)  IRT-5 

IF(ANS.EQ.NO)  IRT-1 
WRITE(4f 165) 

WRITE(4f 265)  ANS 
IF(IRW.EQ.2)  WRITE(1 , 585)  ANS 

IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,555>  ESC, HOME, ESC, CLR 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,355)  ISCL 
IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITE(5,255) 

IF ( IRW . EQ . 2 . OR . IRW . EQ . 3 )  WRITE(5,195> 

280  READ ( IRLUN, * , ERR" 29 0)  XLEN 

IF (XLEN.GE.l.)  GO  TO  300 
290  WRITE(5,155) 

GO  TO  280 

300  IF ( IRW . EQ . 2 . OR . IRW • EQ , 3)  WRITE(5,205) 

310  READ (IRLUN, *,ERR«320>  XBLlNT 

IF (XBLXNT , LE . 0. .OR.XBLINT .GT.XLEN)  GO  TO  320 
GO  TO  330 

320  WRITE<5,155) 

GO  TO  310 

330  IF(IRW.EQ.2.0R.IRW.Eg.3)  WRITE(5,185) 

340  READ (IRLUN, 105)  LEN, ANS 

IF(LEN.EQ.O)  ANS-'N' 

IF (ANS. EQ.YES.OR.ANS.EQ.NO)  GO  TO  350 
WRITE(5, 155) 

GO  TO  340 
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460  IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITEQ, 565) 

470  READ (IRLUN, 575, ERR-480)  LEN, ASIZE 

IF(LEN.EQ.O)  ASIZE-, 1 

IF(ASIZE.LT. .07. OR. ASIZE. GT.. 5)  GO  TO  480 
GO  TO  490 

480  WRITEQ, 155) 

GO  TO  470 

490  WRITE<4,245> 

WRITE<4,355)  ISCL 
WRITE(4,95) 

WRITE(4, 225) 

WRITE<4,285)  IPSCL 

IF (IRW.EQ.2)  WRITEQ  , 405)  IPSCL 

WRITE(4,565) 

WRITE(4, 275)  ASIZE 

XF( IRW.EQ.2)  WRITEC1 , 595)  ASIZE 

IF ( IRW.EQ.2. OR. IRW.EQ. 3)  WRITEC5,555)  ESC , HOME, ESC, CLR 
IF (IRW.EQ.2 .OR. IRW.EQ. 3)  WRITE(5,355)  ISCL 
IF (IRW.EQ. 2. OR. IRW.EQ. 3)  WRITE(5,445) 

IF (IRW.EQ. 2. OR. IRW.EQ. 3)  WRITE(5,315) 

500  READ (IRLUN, 325)  NXTL,XTL 

IF (NXTL.GT. 40)  GO  TO  510 
IF  ( 1RW . EQ . 2 . AND . NXTL . EQ . 0)  WRITEQ  ,.585) 

IF (NXTL. EQ . 0)  GO  TO  589 
GO  TO  520 

510  IF (IRW.EQ. 2. OR. IRW.EQ, 3)  WRITE(5,15) 

GO  TO  500 

520  IF ( IRW . EQ , 2 . OR . IRW . EQ . 3 )  WRITE(5,565> 

530  READ (IRLUN, 575, ERR-540)  LEN, XSIZE 

IF(LEN.EQ.O)  XSIZE-. 2 

IF(XSXZE.LT. ,07 .OR.XSIZE.GT. .5)  GO  TO  540 
GO  TO  550 

540  WRITE(5,155> 

GO  TO  530 

550  IF  (IRW.EQ, 2 .OR. IRW.EQ. 3)  WRITEQ,  225) 

560  READ ( IRLUN, 455, ERR- 57  0)  LEM, IPXTL 

IF(LEN.EQ.O)  IPXTL-1 

IFQPXTL.GT.O.AND. IPXTL.LT. 4)  GO  TO  580 
570  WRITE(5,155) 

GO  TO  560 

580  '  WRITE (4, 245) 

WRITE'4, 445) 

WPTT’Pfi  ne) 

WRITE(4Q05)  (XTL(I) ,1-1, NXTL) 

IF  (IRW.EQ,  2)  WRITEQ, 585)  (XTL  (I)  ,  >1  ,NXTL) 

WRITE(4, 565) 

WRITE(4, 275)  XSIZE 

IF  (IRW.EQ.  2)  WRITEQ, 595)  XSIZE 

WRITE(4,225) 

WRITE(4,285)  IPXTL 

IF  (IRW.EQ,  2)  WRITEQ, 405)  IPXTL 

589  IF ( IRW . EQ , 2 . OR . IRW , EQ , 3)  WRITEQ, 335) 

590  READ (IRLON, 325)  NXTL, XTL 
IF (NXTL.GT. 40)  GO  TO  600 

IF ( IRW . EQ . 2 .AND . NXTL . EQ . 0)  WRITEQ, 585) 
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IF(NYTL.EQ.O)  GO  TO  676 
GO  TO  610 

600  IF(IRW.EQ..2.0R.IRW.EQ,3)  WRITE(5,15) 

GO  TO  590 

610  IF (IRW.  EQ .  2.0R.IRW  .EQ  ,3)  WRITE(5,565> 

620  READ  (IRLtJN,  575,  ERR-63  0)  LEN, YSIZE 

IF(LEN.EQ.O)  YSIZE*  ,2 

IF(YSIZE.LT,.07.OR.YSI ZE.GT. . 5)  GO  TO  630 
GO  TO  640 

630  WRITE(5, 155) 

GO  TO  620 

640  IF (IRW . EQ . 2 . OR . IRW . EQ . 3)  WRITE(5,225) 

650  READ(IRLUN,455,ERR-66  0)  LEN,  I  PYX’D 

IF (LEN .EQ. 0)  IPYTL*1 

IFdPYTL.GT.O.AND.IPYTL.LT.4)  GO  TO  670 
660  WRITE(5,155) 

GO  TO  650 

670  WRITE<4,335) 

WRITE(4, 305)  (YTL(I) ,  I»1,NYTL) 

IF(IRW.EQ.2)  WRXTE(1,585)  (YTL(I) ,I-1,NYTL) 

WRITE(4,565) 

WRITE(4 , 275)  YSIZE 
IF(IRW.EQ.2)  WRITE(1,595>  YSIZE 
WRITE(4,225) 

WRITE(4, 285)  IPYTL 

IF ( IRW. EQ, 2)  WRITE(1,405)  IPYTL 

676  IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITE(5,555)  ESC, HOME, ESC, CLR 

WRITE(4, 245) 

DO  780  >1,4 

IF (IRW. EQ. 2. OR. IRW. EQ. 3)  WRITE(5,345)  I 

ISTRT* (1-1)  *60+1 

IEND-ISTRT+59 

680  READ (IRLUN, 325)  NTTL (I) , (TTL(J) ,J-ISTRT,IEND) 

IF(NTTLd)  .GT.60)  GO  TO  690 
IF(IRW.EQ.2.AND.Mm(I).EQ.0)  WRITEU,585> 

IF(NTTLd)  .EQ.O)  GOTO  780 
GO  TO  720 

690  IF (IRW. EQ. 2. OR, IRW. EQ. 3)  WRITE(5,15) 

GO  TO  680 

720  IF(IRW.EQ.2.0R.IRW.EQ.3>  WRITE(5,235) 

730  READ ( IRLON ,575, ERR*7 4  0)  LEN,TSIZE(I> 

IF  (LEN.  EQ.O)  TSIZE(I)  **,2 

IF(TSIZEd)  .LT..07  .OR.TSIZE(I)  .GT..5)  GO  TO  740 
GO  TO  750 

740  WRITE(5,155) 

GO  TO  730 

750  IF(IRW.EQ.2,OR.IRW.EQ.3)  WRITE(5,175) 

760  READ  (IRLUN, 455, ERR-77  0)  LEN,IPTTL(I) 

IF  (LEN.  EQ.O)  IPTTL(D-) 

IFdPTTL(I)  .GT. O.AND. IPTTL(I)  .LT.4)  GO  TO  771 

770  WRITE (5, 155) 

GO  TO  760 

771  WRITE(4, 345)  I 
JEND*ISTRT+NTTL(I) -1 

IF(NTTLd)  .GT.O)  WRITEU,295)  (TTL(J) , J»ISTRT,JEND) 
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IF(IRW.EQ.2)  WRITEd  ,  585)  (TTL(  J)  ,J-ISTRT,  JEND) 
WRITE<4,235) 

WRITE(4,275)  TSIZE(I) 

IF(IRW.£Q.2)  WRITEd, 595)  TSIZE(I) 

WRITE(4,175) 

WRITE(4,285)  IPTTL(I) 

IF(IRW.EQ.2)  WRITEd, 405)  IPTTL(I) 

780  CONTINUE 

WRITE(5 , 555)  ESC, HOME, ESC, CLR 
NDIG1-2 
NDIG2-2 

IFdSCL.LT. 10)  NDIG1-1 
IFCICRV.LT. 10)  NDIG2-1 
WRITE(5,355)  ISCL 
WRITEC5 , 45)  ISCL 
IFNS-(ISCL-l)  *12+1 
IFNE-ISCL+12 

790  READ (5,55)  LEND, (DNAME(I) ,I-IFNS,IFNE) 

LENDE-IFNS+LEND 
DNAME (LENDE) *0 
CALL  ERRTST(29,K) 

OPEN  ( UNIT* 3 , NAME-DNAME ( IFNS) , READONLY, TYPE ■'OLD1 ,ERR»800) 
CALL  ERRTST(29,K) 

IF (K.EQ.2)  GO  TO  810 

800  WRITEC5 , 155) 

GO  TO  790 

810  WRITE(4,245) 

WRITE(4,355)  ISCL 
WRITEC4, 45)  ISCL 

WRITE(4,305)  (DNAMECI) ,I-IFNS,IFNE) 

DO  820  1-1,24 
READ <3 , 55)  LEN,T 

820  CONTINUE 

NOECD-O 
1-1 

830  READ (3,5)  L£N,XY1(I) 

IFCLEN.EQ.O)  GO  TO  840 
NOE (1) -NOE(l) +1 
>1+1 
GO  TO  830 

840  NOE(l) -NOEC1) -1 

IF C NOE ( 1 ) .LE.1000)  GO  TO  850 
IFCIRW.EQ.2.0R.IRW.EQ.3)  WRITEC5,75) 

GO  TO  9999 

850  DO  860  >1,11 

READ (3 ,55)  LEN,T 

860  CONTINUE 

NOE  (2)  -1 

>1 

870  READ (3,5, QiD-SBO)  LEN,XY2(I) 

NOE (2) ■NOE(2) +1 
>1+1 
GO  TO  870 

880  •  NOE(2)-NOE(2)-l 

REWIND  3 
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DO  890  1-1,24 
READ  (3 , 55)  LEN ,T 
890  CONTINUE 

DO  900  I-l.NOECl) 

READ  (3  ,  *)  LEN , XY 1 ( I ) 

900  CONTINUE 

DO  910  1-1,12 
READ (3, 55)  LEN ,T 
910  CONTINUE 

DO  920  I-l,NOE(2) 

READ (3 ,  *)  LEN,XY2< I) 

920  CONTINUE 

CLOSE  (UNIT-3) 

NOE(l) -NOE(l) /2 
NOE(2)-NOE(2)/2 
XMINT-1 .E27 
YHINT-1.E27 
XMAXT-1  .E-27 
YMAXT-1  .E-27 
DO  930  1-1 ,NOE(l) 

IF(XYKI).LT.XMINT)  XHINT-XYKI) 

IF<XY1(I)  .GT.XMAXT)  XNAXT-XY1 ( I) 

IF(XYl(NOE(l) +1) .LT.YHINT)  YHINT-XYKNOSCl) +1) 
IF(XYKNOEd)  +1)  .ST.YNAXT)  YNAXT-XYKNOE(l) +1) 

930  CONTINUE 

DO  940  I-l ,NOE(2) 

IF ( XY 2(1) .LT.XMINT)  XMINT-XY2<I> 

IF<XY2(I) .QT.XNAXT)  XMAXT-XY2U) 

IF(XY2(NOE(2)  •*■!>. LT.YNINT)  YHINT-XY2 (NOB  2)*1 
IF(XY2(NOE(2) +1) .OT.YNAXT)  YNAXT-XY2(NOE(2) +1) 

940  CONTINUE 

IFUXHM.EQ.l)  GO  TO  950 

XKIN-XMINT 

XHAX-XMAXT 

950  IF(IYNK.EQ.l)  GO  TO  9«0 

YNIN-YNINT 
YHAX-YHAXT 

960  WRITE(4,245) 

DO  990  ICRV-l , 2 

IFdCRV.BQ.l  .AND.  tIRV.EQ.2.0R.TRW,EQ.3)  )  «dTE(5,135) 
IF(ICRV.EQ.2.AND. (IRW.EQ.2.0R.IRW.EQ.3) )  WRITE(5,145> 
IF( IRM.EQ.2.0R.IRN.EQ.3)  WRITE<5,225) 

970  READ ( IRLUN, 455, ERR-98 0)  LEN,  IPCRV dCRV) 

IP(LEN.EQ.O)  IPCRV (ICRV) -1  _ 

IFUPCRVUCRV)  .GT.O. AND. IPCRV (XCRV)  ,LT.4»  GO  TO  981 

980  WRITE(5,155) 

GO  TO  970 

981  IF(ICRV.EQ.l)  WRITE(4,135) 

IF( ICRV.EQ.2)  WRITE(4, 145) 

WRITE (4, 225) 

WRITE(4 , 285)  IPORV(ICRV) 

IF( IRW.EQ.2)  WRITEd ,405)  IPOWdCRV) 

990  CONTINUE 

IF ( IRW.EQ.2 .OR. IRW.EQ.3)  WRITE(5,55S)  ESC, HONE, ESC, CLR 
XQIBB-XBLINT/.25 
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YQIBB-YBLINT/.25 
CALL  Ba}STP(OTYPE,-2,OHIT.0. 1) 
rVr  }SCX'£0‘0)  GO  TO  1000 
10nnt'L  ™ •  WM.XLEN.XQibb) 

10°°  „  XPtISCY.BQ.0)  GO  TO  1010 

CALL  BOPTMM(VNIN,yMAX,yLEN,TOIBB) 
0X11  "WP8IKIWCU 
DO  1020  JPI-1.IRT 
FG»PLOAT(JPX-l) 
TFA-«.2831853E0*FJ/5. 
AX«.008*COS(TPA) 

AY-.0Q8*SIN(TFA> 

IP ( IRT . EQ . 1 )  AX-0.0 
IPdRT.EQ.l)  AY-0.0 
CALL  PLOT(AX,AY,-3) 

law 

)/2.+1.5 

£**•”*■1  *  5*  J55***  u  •  ***$*«>  > 

RXYTL-1,5-(YTLRX*(.5*YSIXE) ) 
RYYTL-(YLBH-(8m*YSIXBl)/2  *-1 

RYTTt ( 4 1 -YLEK+i . 2 

S5  ? J  *52’  *i .S**SX*8 (3) 

00  1040  r-1,4 

00  1040  JPI*1,1RT 

PJ-PLOATtJPI-l) 

TPA-i^BJUSSBO^PJ/S. 

AX».008*COS(WA) 

AY-.OM*SXJ»(trA> 

IPUM.CQ.l)  AX-0,0 
IPURT.BQ.l)  AY-0.0 
-31 

’ 0)  °°  w  10,0 

Zu&&r„%j°  “« 

oo  ioso  jpi«i*mt 

PJ-PtOATtJPt-l) 

WftNf . 283 18S3C0*PJ/S. 
AX-.008*COS(T?A) 

AY-t008*SXN(TfA) 

AX-0.0 

ipaxr.BQ.u  AY-o.o 
CAU  PUJT<AX,AY,-»J) 

low 

‘l&.SP®®#'*’  “  ” 1M1 

00  1080  JM-l,m? 
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FJ-FLOAT( JPI-1) 

TFA-6. 283 1853E0*FJ/5. 

AX* . 0  08 *COS  <TFA ) 

AY*,008*SIN(TFA) 

IFtIRT.EQ.l)  AX-0.0 
XF(XRT.EQ.l)  AY-0.0 
CALL  PLOTtAXiAY , -3) 

IF(NXTL.EQ.O)  GO  TO  1060 

CALL  8<2iTTL(XTL,NXTL,XSI2E,RXXTL,RYXTL,l.,0.) 

1060  COOT  IN  UE 

1061  RXTTL3-XLEN+3 . 

CALL  NEWPEN ( IPXTL) 

DO  1070  JPI-1, IRT 
FJ-PLOATtJPI-1) 

TFA-6 .  283 18  53  E0  *FJ/  5 . 

AX* . 0  08  *COS l TFA ) 

AY* . 008*SIN (TFA) 

IF( XRT.EQ.l)  AX-0.0 
IFtIRT.EQ.l)  AY-0.0 
CALL  Pl>QTtAX,AY»-3) 

CALL  BG8TTt.tTTL3 ,15,  .14 , RXT7L3 ,4, 0,1. ,0.) 

CALL  BGNTTLtTTL4, 20, .14  , RXTTL3 , 3 . 7 , 1 .  .0.) 

CALL  SYMBOL tRXTTL3 ,4.  07  , .14 ,1 , 0.0, -1) 

CALL  SYMBQLtRXTTL3»3.77,,14»24»0.0»-l) 

1070  COOTIHUE 
CALL  XEVPCNtlPCRVtl)) 

OO  1080  JPI-1* IRT 
FJ-FLOATtJPt-1) 

TFA-6. 283 1I53E0*FJ/S. 

AX-.008*C05tTPA) 

AY-. 008*318 tTFA) 

IFtIRT.EQ.l)  AX-0.0 
IFtIRT.EQ.l)  AY-0.0 
CALL  PLOT) AX, AY, -3) 

CALL  &PLOTtI$CL,MOEU>  ,0,XYl,mt80EU) *1)  ,1.0 
1  l.SYKB, .14) 

CALL  80VPLT  t ISCL) 

1080  COOT  XX  Uf 

CALL  MWPEXU) 

CALL  8EOTC8tIFCStV{2) ) 

00  1010  JPX-l,XRY 
FJ*PLOATtJPt«l) 

TFA-6 . 2S3 1853£0*FJ/S . 

AX-.008*COStTFA> 

AY-. 008*818 tTFA) 

IFtIRT.EQ.l)  AX-0.0 
IFtIRT.EQ.l)  AY-0.0 
CALL  PLOTtAX*AY*-3) 

CALL  lPU>TttSCL*80S(») ,l.XY2,XY2lliOtt2)  *1) .1*0 
1  0.0..0.) 

CALL  BOVPLTtlSCL) 

1040  COOT  XROE 

CALL  SEMPER  tl) 

CALL  MOPLTtO) 

trtIRH.CQ.l.OR.IRM.EQ.2)  CALL  CLOSE  (1) 
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CALL  CLOSE  (2) 

CALL  CLOSE  (4) 

WRITE(5,S55)  ESC, HOME, ESC, CLR 
WRITE (5, 35) 

1100  READ(5,105)  LEN,ANS 

IF(LEN.EQ.O)  ANS-'N’ 

IF (ANS.EQ.YES.OR.ANS.EQ.NO)  GO  TO  1110 
WRITE(5,155) 

GO  TO  1100 
1110  COOT  IN  UE 

IF!ANS.EQ„YES)  GO  TO  11 
CALL  NEWPEN(l) 

CALL  BNDPLT(l) 

WRITE (5 ,555)  ESC, HOME, ESC, CLR 
WRITE (5, 25) 

9999  CALL  EXIT 

5  FOSMAT(Q,53X, F14.ll) 

15  FORMAT! ‘  You  must  enter  60  characters  or  less.') 

25  FOIMAT! ///////////» 29X, 'Successful  completion. * *//////////> 

35  FORMAT !///////////,'  Create  another  graph  (Y/Nl? 

45  FORMAT!/, '  Input  file  for  Graph  ’,I<NOIG>,'  IfilespecJ?  *,i) 

55  F0RMAT!Q,13A1) 

65  FORMAT!///, '  Output  file  tfileapec)?  >,S) 

75  FORMAT!//,'  Your  data  file  contains  more  than  1000  X/Y* 

1  '  coordinate*. ',//, •  Program  stop.  No  output  produced.') 

•5  FOIMAT!  14) 

95  FORMAT! 'at  x  and  Y  axes  specifications! ') 

105  FORMATtQ.lAl) 

115  FORMAT!'  That  file  does  not  exist.') 

125  FORMAT!//,'  Plot  sptci flcationsi ' ) 

135  FORMAT!//,*  "Actual  data"  curve.') 

145  FORMAT!/,1  "Estimated  values"  curve.') 

155  FORMAT!'  What)?  *,l) 

165  FORMAT!/,'  Rat  race  Option  (Y/Nl?  ',») 

175  FORMAT!'  (1)  SLACK  (2)  ROD  (3)  GREEN',/, 

1  '  Color?  '»>) 

Its  FORMAT!/,'  Scale  option  (Y/Nl?  ',*) 

195  FORMAT!/,'  X-Axis  length  in  inchta?  ',*) 

20S  FORMAT!/ • '  Distance  be  tween  blips  in  inches?  ',*) 

as  FORMAT!/,'  Y-A*ie  length  in  inches?  ',#> 

225  FORMAT!/,'  (11  BUCK  (2)  RED  (31  GREEN'./, 

1  '  Color?  ',*) 

235  FORMAT! '  Character  Site  n  Inches  (.07-.51?  ',») 

245  FORMAT!/,'  » t 

^  1  •s«»»M»sm#e*ewee«»«e«e»Mse«ietM!ees«»e«»«»e * 

255  FORMAT! '»t  X  axis  specifications!',/) 

265  FORMAT! '♦*»1A1) 

275  FOJMA.i  '♦  * ,F4S . 7) 

265  FOIMAT! '♦* ,ill) 

295  FOIMAT! '  '  ,60A1) 

305  FORMAT! 'a*,40Al) 

315  FORMAT!'  X-Axis  Titls  (1-40  chars.)?  ',») 

325  FORNAT(Q,60A1) 

335  FORMAT!/,’  Y~Axi*  Title  tl-40  chert. »?  *,il 

345  FOWATC  Graph  title  (line  Ml,'!  (1-60  chars. 1 7  ') 


■  sawKis  was*****- 
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355  FOWATC  Graph  ntsaber  ',I<NDIG>,$> 

365  FORMAT!/,'  (11  linear  (21  Seai-log  13]  Log-log' 

1  Type  of  scale?  ',*) 

375  FORMAT! Y  axis  speci f ications: ' ,/> 

385  FORMAT* '+ ' ,E15.7) 

395  FOIMAT!/,'  Number  of  cycles  {>01?  ',$) 

405  FORMAT!  XI) 

415  FORMAT < '  tl)  Read  responses  fraa  GRAEB1.COM 

1  '  and  do  not  updar.e  GRAEB1 .COM. ' ,//, 

2  '  123  Read  responses  f  "a  keyboard  and',/, 

3  '  update  GRAiBl.CUM.  * ,//, 

4  '  331  Read  response  from  keyboard  and',/, 

5  '  ao  not  upda.  GRAPB1.COM. ',//, '  Option?  ',$) 

425  FORMAT (F15 .7) 

435  FORMAT(E15.7> 

445  FORMAT ( '+  Title  s:  ci f ications: ' ,/) 

455  F0RMAT(Q,I1) 

4*;.  FORMAT!/,'  X-Axis  mini  bub/ maximum  overide  option  tY/Nl?  ',$> 

475  FORMAT!/, '  Y-Axis  ainiaum/maxiaum  overide  option  IY/N] ?  ',$i 

485  FORMAT  ( '  Minin:*  value?  ',?> 

495  ?0*tIAT(  '  Maximum  vtJue?  ',$) 

555  FORMAT,"  *,4A1) 

£65  FORMAT!/,'  Character  size  in  inches  I.07-.51?  ',$) 

575  FORMAT!Q,F5.4) 

535  FORMAT ( 6 0A1) 

595  FORMAT (F5. 3) 

END 

SUBROUTINE  BOPTMM( XMIN, XMAX , AXLEN , QIBB) 


C  BOPTMM  will  determine  a  scale  for  plotting  values  between  XMIN  and  XMAX 
C  such  that  each  l/20th  inch  tic  mark  has  an  easily  READ  decimal  value. 

C  If  QIBB  specifies  tic  marks  on  the  axis,  '0.'  must  fall  on  a  blip  even 
C  if  off  the  plot  page.  AXLEN  is  truncated  to  the  nearest  l/20th  inch  and 
C  XMIN  and  XMAX  are  changed  to  the  endpoints  of  the  axis  to  be  drawn. 

C 

C  Arguments: 

C 

C  XMIN  »  Minimum  value  in  array  to  be  plotted. 

C 

C  XMAX  ■  Maximum  value  in  array  to  be  plotted. 

C 

C  AXLEN  ■  Axis  length  in  inches. 

C 

C  QIBB  -  Number  of  quarter  inches  between  tic  marks  on  axis. 

C 

C  If  QIBB  is  positive  t>0]  -  fit  scale  to  XMIN  and  XMAX 

C  using  tic  mark  interval. 

Q 

C  If  QIBB  is  negative  £<01  -  permit  a  2.5%  overscale  on 

c  either  end  using  tic  marks. 

C  _ 

- - 


REAL *8  XL0W,XEI, POWER 
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IXHI  -XXHI  /10 
FXMAX-FXHAX/10. 

IEXP-ISXP+i 
GO  TC  120 

130  IF  (IZ.EQ.O)  GO  TO  140 

IX-IABS ( IXuOW-lO* ( IXLOW/1 0) ) 

I?  (IX.EQ.O)  GO  TO  140 
IF  (IXLOW.LT. 0)  IX-10-IX 
IL-IXHI-IX 

IF  (FXMAX.GT.IL) GO  TO  140 

IXHI-IL  { 

IXLOW-IXLOW-IX 

GO  TO  120  !■ 

140  POWER-10. D0**IEXP  ! 

XLOW-IXLOW*POWER  l' 

XHI  -IXHI  *  POWER 

IF  t RANGE. LT.O.)  GO  TO  ISO  i 

XMIN-XLOW  b 

XMAX-XHI 
RETURN 

ISO  XKIN-XHI 

XFAX-XLOW  p 

RETURN 

5  FOHMATf  Invalid  BOPTMM  arguments: ' ,//, 

1  '  Minimum  value  *  ',F14.7,/, 

2  '  Maximum  value  «  1 ,F14.7,/,  ■» 

3  *  Axis  length  ■  '  ,FI4.7,/,  '■ 

4  •  Blip  intervals  «  ’,P14.7,//,'  Program  stop.') 

END 
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12*0 

RANGE* XMAX-XKIN 
IF  (RANGE. NE.O.)  GO  TO  10 

IF(IRW.£Q.2)  WRITE(5,5)  XMIN, XMAXrAXLEN,QIBB 
CALL  EXIT 

10  R-ABSI RANGE) 

IF  (RANGE. GT.O.)  GO  TO  20 

X-XMIN 

XMIN-XMAX 

XMAX*X 

20  AXL20-INT(AXLEN*20.) 

AXLEN-AXL20/20. 

UNP20-R/AXL20 

IF  (CIBB.GT.O.)  GO  TO  30 

X*.Q25*R 

XNAX-XMAX-X 

XHIN-XMIN+X 

R-XMAX-XMIN 

UNP20-R/AXL20 

3C  TWPBLP«INT(ABS(QIBB) )*5 

IF  (TWPBLP.EQ.O.)  GO  TO  40 
IF  (AXL20/TWPBLP.GT. 1.0001)  GO  TO  50 
IF  (R.LE.ABS(XMAX+XKIN))  GO  TO  50 
40  TWPBLP-1. 

IZ-1 

50  IEXP— 1 

IF  (ONP20.GE.1.)  IEXP*9 
60  X-0. 

Y1-10.**IEXP 
70  X-X+l . 

FONP20«X*Y1 

IF  (UNP20.GT.FONP20)  GO  TO  70 
IF  (X.NE.l.)  GO  TO  80 
IF  (ONP20.EQ.PUNP20)  GO  TO  80 
IEXP-IEXP-1 
GO  TO  60 
80  IX*X 

IAXL20-AXL20 
tTWPBL-TWPBLP 
FXMAX-XMAX/Y1 

90  IUNPBL*ITWPBL*IX 

Y-XMIN/(IUNPBL*Y1) 

IY-Y 

IF  (Y.E (•  IY)  GO  TO  100 
IF  (Y.LT.O.)  IY-IY-1 
100  IXLOW«IY*IUNPBL 

IXHI-IXLO;/ .’.20* IX 
IF  (FXMAX.LT  IXHI)  GO  TO  110 
IX-IX+1 
GO  TO  90 

110  IL-IONPBL*INT(  (FL0AT( IXHI)  -FXMAX)  /<2*IONPBL) ) 

IXLOW-IXLOW-IL 
IXHI  -IXHI  -IL 

120  IF  <1 0*(IXLOW/10) .NE.IXLOW.OR.10*(IXHI/10) .NE.IXHI)  GO  TO  130 

IXLOW-IXLOW/1 0 
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CALCOMP  Generated  Plot  (Semi logarithmic)  of  the  GRPHX.XXX  Data 
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KEbPIMM  UUANI 1IAT1VE  Ell  lEbllllb 

1I1S1 KUL I i UNb  R)K  IK lU(i  THE  CUMPUTEK  lEWWllALS  111  USAFSATVVN  TO  PKOlESS  lift.  DATA  COLLECTED  01  THE  SALT  FOG 
lIlblKUTlLNl 

KEY:  -  COMPUTER  GENERATED  INFORMATION  (CRT  SCREEN) 

-  Piogiam  LUei  Ge.ne.ita.tzd  Intimation  \Enteied  \>ia  the  keyboard;  intimation  diiptaijed  on  the  CRT 
6cieenl 

-  Sequential  Step  Numbers;  not  displayed  on  CRT  screen 


—  Comments  to  help  the  program  user 

CRT  SI  KEEN  DISPLAY 

1)  >  _ 

2)  >  HELLO  _ 

3)  ACCOUNT  OR  NAME;  _ 

A)  ACCOUNT  OR  NAME:  XXXXXX_ 

5)  PASSWORD:  _ 

6)  PASSWORD:  XXXX  _ 

.  (SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 

ENTER  TERMINAL  TRACKING  NUMBER  (WHITE  (AG  ON 
RIGHT  FRONT)? 


CONTENT 

CRT  display  status  normally  found  on  an  Idle 
terminal. 

Type  in  'HKl.LO';  depress  "Return"  key  on  keyboard. 
Computer  response. 

Type  in  your  lust  name;  depress  “Return"  key  on 
keyboard. 

Computer  response. 

Type  in  yonr  password;  password  characters  are  not 
displayed  on  CRT  screen;  depress  "Return"  key  on 
keyboard. 

No  response  required  on  your  part. 

Last  message  of  the  group. 
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LKl  blKttN  DlbPLAY 


7)  ENTER  TERMINAL  NUMBER  (WHITE  TAG  ON  RIGHT 
FRONT)?  XX  _ 


(SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 


>_ 

8)  >  Run  I)R3 :  [  305 , 4JNACLKQFT 


(PROGRAM  STATEMENTS) 


Qrtltlll 

Enter  two  dibits;  depress  "Return"  key  on  keyboard. 
No  response  required  on  your  part. 

Computer  Is  ready. 

Type  In  "Run  l)K3:  i  30b, 4  jNACLIUjFT;  depress  the 
"Return"  key  on  the  keyboard. 

No  response  required  on  your  part 


9)  ENTER  THE  FOLLOWING:  UUl  FOR  THE  FIRST  DATA  La»t  °f  program  statements. 

SETj  UU?  FOR  THE  SECOND  DATA  SET;  UU3  FOR  THE 

THIRD  DATA  SET,  ETC* 

10)  eniky  =  ^X  Enter  a  three  di^lt  mir.tber  to  name  this  set  of 

.  ~  data.  Depress  the  "Return"  key  on  keyboard. 


U) 


(PROGRAM  STATEMENTS) 


No  response  required  on  your  part. 


12)  ENTER  the  FOLLOWING!  UUl  FOR  THE  FIRST  RESIDUAL  Usc  uf  Pr°l5r«‘l;'  Statements. 
set;  UU?  FOR  THE  SECOND  RESIDUAL  SET;  UUJ  FOR 
THE  THIRD,  ETC. 

12)  ENTRY  =  XXX  _ 


Enter  a  three  dip. It  number  to  name  this  set  of 
residuals;  depress  the  "Return"  key  on  the 
keyboard. 
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LKl  bLKttN  lllbPLAY  UlULHf 

14)  . 

.  (PROGRAM  STATEMENTS)  No  response  required  on  your  part. 


15)  ENTER  THE  FOLLOWING  001  FOR  THE  FIRST  GRAPH  Last  of  program  statements. 

SET;  UU2  FOR  THE  SECOND  GRAPH  SET;  003  FOR  THE 

THIRD,  ETC- 

16)  £NTRY  3  XXX  _  Enter  a  three  digit  number  to  name  this  graph  set; 

depress  the  "Return"  key  on  the  keyboard. 

17)  ENTER  THE  NUMBER  OF  SODIUM  CHLORIDE  CALIBRATION  There  are  seven  sodium  chloride  calibration  con- 

CONCEHTRATION  STANDARDS-  centration  standards.  Enter  the  number  '7'; 

depress  the  “Return"  key  on  the  keyboard. 

ENTRY  =  ^  _ 

l8'  DEPRESS  RETURN  KEY  AFTER  A  VOLTAGE  MEASUREMENT  No  response  required  on  your  part. 

19)  lihER  the  DATA  POINTS 


SAMPLE  NUMBER 

1 

2 

5 

‘I 

b 

U 

/ 


VOLTAGE 
MEASUREMENT 
(X  DATA) 

_XXX_ 

XXX" 

XXX 

-XXX~ 

XXX 

~XXX" 

xxx~ 


SODIUM  CHLORIDE 
CONCENTRATION 
(Y  DATA) 

1.0 

0.1 

"0.0T 

~O.OOT 

~o.ooor 

~O.OOOOT 

'o.ooooor 


Enter  a  voltage  measurement,  then  depress  the 
“Return"  key  on  the  keyboard.  Type  the  correspond¬ 
ing  sodium  chloride  calibration  concentration  then 
depress  the  “Return"  key  on  the  keyboard.  Repeat 
this  procedure  until  all  data  points  have  been 
entered. 
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Ik  I  LI.KLIJJ  III  jflAY 

20)  SUBJECT  NAME:  XXXXX_ 

21)  TYPE  CF  MASK:  _XXXXX_ 

22)  DATE  TESTED:  _XXXXX_ 

23)  TIME  TESTED:  XXXXX_ 

24)  THE  USER  IS  FREE  TO  SELECT  ONE  OF  TWO  GROUPS 
OF  EXERCISE  PROTOCOLS- 

-  (PROGRAM  STATEMENTS) 

25)  TO  SPECIFY.  THE  EXERCISE  PROTOCOL  llROUP  OF 

Interest,  type  either:  Group  l  or  Group  2 

ENTRY  =  CiYOup  x 

26)  - 


LOfTILM  I 


Enter  subject's  name  (up  to  30  characters  long); 
depress  the  "Return''  key  on  tlx;  keyboard. 

Enter  mask  nomenclature  (up  to  30  characters  long), 
depress  "Return"  key  on  keyboard. 

Enter  date  subject  was  tested;  depress  the  "Return” 
key  on  tlx;  keyboard. 

Enter  tlx;  time  ot  djy  subject  was  tested  (for 
F.xainple:  1430  hours);  depress  "Return"  key  on 

keyboard. 

First  of  a  series  of  program  statements  explaining 
the  two  groups  of  exercises.  No  response  required 
on  your  part. 


Last  of  program  statements. 

Enter  either  ’(.'roup  1’  or  'Croup  2'j  depress  the 
"Return"  key  on  the  keyboard. 


The  program  statements  explain  how  to  enter 
Integrator  data.  No  response  on  your  part. 


(PROGRAM  STATEMENTS) 
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IK  I  .>t,KLLU  UlbltAY 


ILI'VlLIJl 


21)  EXERCISE  COUNT  DATA: 

INTEGRATOR  TIME  PEKlOU 
EXERCISE  COURT  (IN  SECOIJDS) 

NORMAL  BREATHING 

STRAIGHT  AHEAD  XXXXXX  _XX_ 


Enter  Integrator  cxurclsc  count  data  as  previously 
specified  In  Hie  program  statements;  depress  the 
''Return"  key  on  the  keyboard  after  each  Integrator 
count  value.  Kilter  the  time  period  for  the  exer¬ 
cise  using  no  decimal  points;  depress  the  "Return" 
key  on  the  keyboard.  Repeat  until  data  for  each 
exercise  lias  been  input. 


( PROGRAM  STATEMENTS) 


The  program  statements  explain  the  functional 
definitions  of  the  variable  (X);  no  response  Is 
required  on  your  part. 


2D)  TO  SPECIFY  A  FUNCTIONAL  DEFINITION  OF  (X) 
SELECT  THE  CORRESPONDING  NUMBER  INSIDE  THE 
BRACKETS 

ENTRY  =  3 


30)  ENTER  THE  ORDER  OF  THE  DESIRED  POLYNOMIAL 
(MAXIMUM  =  THE  NUMBER  OF  SODIUM  CHLORIDE 
CALIBRATION  CONCENTRATION  STANDARDS  ”1) 

ENTRY  =  6 


31) 


(PROGRAM  STATEMENTS) 


Last  of  program  statements. 


Enter  the  number  depress  the  "Return"  key  on 

the  keyboard. 

Enter  the  number  '6';  depress  Che  "Return"  key  on 
tile  keyboard. 


A  list  of  coefficient  numbers  Is  printed  on  the  CRT 
screen;  no  response  required  on  your  part. 


32)  Do  YOU  WANT  A  LIST  OF  RESIDUALS? 

(answer  Yes  or  no)  Ye* 


Enter  'Yes’,  unless  otherwise  Instructed;  depress 
the  "Return"  key  on  the  keyboard. 
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luntui 

A  list  of  residuals  Is  displayed  on  the  CKT  screen. 
No  response  required  on  your  part.  If  you  want  to 
stop  Che  flow  of  information  on  the  CKT  screen, 
depress  simultaneously  the  'CTRL’  and  'S'  keys  on 
the  keyboard.  To  resume  the  flow  of  information  on 
the  CRT  screen,  depress  simultaneously  the  'CRTL' 
and  ’ Q *  keys  on  the  keyboard. 

The  descriptive  and  protection  factor  calculations 
are  displayed  on  CRT  screen;  no  response  required 
on  your  part. 


35)  UO  YOU  WISH  TO  USE  THE  SAME  SODIUM  CHLORIDE  in  or  ""  depending  upon  your  desires; 

STANDARD  CALIBRATION  VOLTAGE  MEASUREMENTS  AND  depress  the  Return  key  on  the  keyboard. 

EXERCISE  INTEGRATOR  COUNT  DATA,  BUT  CALCULATE  A  -  lf  You  *71“  1,1  IJu  tlle  result  will  be:  36) 

DIFFERENT  DEGREE  OR  FORM  OF  THE  POLYNOMIAL  CURVE  ~  lf  You  W0  1,1  Yes  the  result  wiu  1165  29> 

FITTING  FUNCTION?  (ANSWER  YES  OR  Uo)  _XXX_ 

36)  1)0  YOU  WISH  TO  CALCULATE  PROTECTION  FACTORS  FOR 
ANY  OR  ALL  OF  THE  FOLLOWING  CONDITIONS: 
i)  A  DIFFERENT  SUBJECT  2)  A  DIFFERENT  SET  OF 
SODIUM  CHLORIDE  CALIBRATION  STANDARD  VOLTAGE 
MEASUREMENTS  5)  A  DI FFERENT ' SET  OF  EXERCISE 
INTEGRATOR  COUNT  DATA?  (ANSWER  YES  OK  llo)  _XX*_ 

37)  >  PRINT  DATA. XXX  If  you  want  a  computer  print  out  on  paper  of  the 

data  that  you  input,  type  ’PRINT  DATA.XXX',  where 
XXX  is  the  number  you  entered  from  the  keyboard  in 
step  9);  depress  the  "Return"  key  on  keyboard. 


Type  in:  'Yes'  or  'No'  depending  upon  your 
desires;  depress  the  "Return"  key  on  the  keyboard. 

—  If  you  type  in  'No'  the  result  is  37) 

—  If  you  type  in  'Yes'  the  result  is  9) 
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IKl  bLKhLN  UlbPLAY 


lUfv'iLIJ  i 


38)  >  PRINT  CXLCX.XXX 


39)  >  PR/VT  ORPHX.XXX 


40)  >  BYE 


If  you  wont  .i  cm iput er  print  out  on  |vtper  of  the 
results,  type  '  PRINT  CAI.CX.XXX',  where  t  lie  first  X 
of  CAI.fIX.XXX  is  assigned  the  number  I  hy  t he  com¬ 
puter  and  incremented  by  I  each  time  <i  form  of  the 
polynomial  curve  fitting  function  Is  calculated, 
and  XXX  is  tlie  numher  you  entered  from  the  keyboard 
In  step  13);  depress  the  "Return”  key  ou  the  key¬ 
board. 

If  you  want  a  computer  print  out  on  paper  of  the 
graphing  data,  type  'PRINT  CRPIIX.  XXX '  ,  where  the 
first  X  of  (IK  I'll  X.  XXX  Is  assigned  the  number  I  by 
the  computer  and  Incremented  by  I  each  time  a  form 
of  tlie  polynomial  curve  fitting  function  is  calcu¬ 
lated,  and  XXX  Is  the  numher  you  entered  from  the 
keyboard  in  step  lb);  depress  tl>e  "Return"  key  on 
the  keyboard. 

If  you  are  finished  with  the  computer,  type  In: 
’BYE'  from  Che  keyboard;  depress  tlie  "Return"  key 
on  keyboard  and  you  will  he  automatically  logged 
off  the  computer  terminal.  Computer  printed 
results  can  he  picked  up  In  computer  center,  Bldg 
150,  approximately  2  hours  after  you  log  off  the 
computer  terminal. 


0s 

c? 


^r 

So 
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Z -J1  HLU1  f  Hlh  KOUTINL 

ll&lKtU I IDilb  t-(«h  W>IN<»  THt  CfMPIIfLK  iLlfllNALS  111  INlFS/Tl/VU  Tl)  HJ)I  A  (.KAMI  IN  IWU  DIMtHSlUNS- 

KtY :  -  COMPUTER  GENERATED  INFORMATION  (CRT  SCREEN) 

-  Pioctfam  ll-ie  1  Generated  Intoimatior.  {Entered  via-  the  henhoaid;  information  di-iptaued  on  the.  APT 


icreeni 

-  Sequential  Step  Numbers;  not  displayed  on 

_  Comments  to  help  the  program  user 

LKI  SlkLLN  DISPLAY  . 

1)  >  _ 

Z)  >  HEUO  _ 

3)  ACCOUNT  OR  NAME:  _ 

4)  ACCOUNT  OR  NAME:  XXXXXX_ 

5)  PASSWORD: 

h)  PASSWORD:  XXXX  _ 

.  (SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 


CRT  screen 


caiiLin 

CRT  display  status  normally  found  on  an  Idle 
terminal. 

Type  In  'HELLO';  depress  "Return"  key  on  keyboard. 
Computer  response. 

Type  In  your  last  name;  depress  "Return"  key  on 
keyboard. 

Computer  response. 

Type  In  your  password;  password  characters  are  not 
displayed  on  CRT  screen;  depress  "Return"  key  on 
keyboard. 

No  response  required  on  your  part. 


ENTER  TERMINAL  TRACKING  NUMBER  (WHITE  TAG  Oil  RIGHT  FRONT )?_  l-ust  nos an go  of  the  group. 
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Lkl  SCKllifJ  DISPLAY  ClMMHJl 


7)  BULK  TERMINAL  TRACKING  NUMBER  (WHITE 
TAG  ON  RIGHT  FRONT)?  XX_ 


(SERIES  OF  COMPUTER  SYSTEM  messages) 


8)  Run  VR3:[30S,4]HACLGRAPH 


9)  '111  READ  RESPONSES  FROM  GRAPHl -COM' AND  DO 
NOT  UPDATE  GRAPHl.COM. 

121  READ  RESPONSES  FROM  KEYBOARD  AND  UPDATE 
GRAPHl.COM. 

I  5 1  HEAD  RESPONSES  FROM  KEYBOARD  AND  DO  NOT 
UPDATE  GRAPHl.COM. 

UP11UN?  *_ 


10)  OUTPUT  FILE  IfilespecIV  xxxxxx.xxx 


U)  GRAPH  NUMBER  (I 

X-AXIS  HI Ml  MUM/MAX I  HUM  OVERRIDE  OPTION 

IY/HJ?  V 


Enter  two  digits;  depress  "Return"  key  on  keyboard. 


No  response  required  on  your  part. 


Computer  is  ready. 

Type  in  "Run  Dll'):  1 305,4  )NACLGRAPH ;  depress  tbe 
“Return"  key  on  keyboard 

Computer  response  explaining  the  three  graphing 
options  of  the  program.  No  response  required  on 
your  part. 


Enter  the  number  * l *  ’2’  or'V;  depress  the 
"Return"  key  on  the  keyboard.  If  you  enter 
the  number  '2*  or  'V  the  result  is  10).  It  you 
enter  the  number  ’1‘  the  result  la  S9). 

Enter  the  name  u(  the  output  file  in  the  form 
'XXXXXX.XXX';  the  file  name  may  have  up  to  nine 
characters  before  the  Period  ana  must  have  three 
characters  after  the  period.  For  example, 
' PLOT. DU l ' l  depress  the  "Return"  key  on  the 
keyboard. 

Enter  the  teller  'Y';  depress  the  “Return"  key  on 
the  keyboard. 
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12)  MINIMUM  VALUE?  XXX_ 

13)  MAXIMUM  VALUE?  XXX_ 

14)  Y-AXIS  MINIMUM/MAXIMUM  OVERRIDE  OPTION 

t  Y/M I? 

15)  MiMIMUM  VALUE?  XXX_ 

16)  MAXIMUM  VALUE?  XXX_ 

1?)  PLOT  SPECIFICATIONS: 

111  LINEAR  |?1  SEMI-LOG  lil  LOG-LOG 
TYPE  OF  SCALE?  ?_ 

IS)  NUMBER  OF  CYCLES  l >01?  7_ 

19)  RETRACE  OPTION  tY/NlV 

20)  graph  number  II:  x-axis  specifications: 

X-AXIS  LENGTH  IN  INCHES?  *X_ 

21)  DISTANCE  BETWEEN  BLIPS  IN  INCHES?  XX_ 


UWItlJl 


Enter  the  smallest  value  on  the  X-Axis;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  largest  value  on  the  X-axis.  Depress  the 
“Return"  key  on  t he  keyboard. 

Enter  the  letter  'Y';  depress  the  "Return"  key  on 

the  keyboard. 

Enter  the  smallest  value  on  the  Y-axis;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  largest  value  on  the  Y-axis;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  number  '2'  unless  otherwise  instructed; 

depress  the  "Return"  key  on  the  keyboard. 

Enter  the  number  unless  otherwise  Instructed; 

depress  the  "Return"  key  on  the  keyboard. 

Enter  the  letter  * Y * ;  depress  the  “Return"  key  on 

the  keyboard. 

Enter  the  length  of  the  X-axis  in  Inches 
number  not  a  wurd;  depress  the  “Re- 
keyboard. 

Enter  the  distance  In  in-- ' 

X-axis.  Enter  a  no- 
“Return"  key  un  tlie 
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22)  SCALE  UPriON  IV/iJl?  w_ 


23)  GRAPH  IIIIMBCR  It:  Y-AXIS  SPECIFICATIONS: 
Y-AXIS  LENGTH  IN  INCHES?  **  _ 


24 )  DISTANCE  HETWEEN  BLIPS  III  INCHES?  **_ 


2b)  SCALE  OPTION  lY/IJl?  W 


26)  liRAPH  NUMBER  IJ :  X  AND  Y  AXIS  SPECIFICATIONS: 
111  BLACK  I 2 I  RED  I  5  I  GREEN 

COLOR?  L 

27)  CHARACTER  SI2E  IN  INCHES  I  - U/  *  -b  I? 

2ti)  GRAPH  NUMBER  IJ  TITLE  SPECIFICATIONS: 

X-AXIS  TITLE  ll-'IIJl?  XXXXX_ 

2lJ)  CHARACTER  SIZE  IN  INCHES  |.U/-.‘j|?  XX_ 

30)  !  J I  I  Slack  121  Ked  1 5 1  (ireen 

COLOR?  L 


urni.jji 


Enter  the  letter  *N'  unless  otherwise  instructed; 
depress  the  "Return''  key  on  the  keyboard. 

Enter  the  length  <>l  the  Y-axis  in  Inches.  Enter  a 
number  not  a  word;  depress  the  "Return"  key  on 
keyboard. 

Enter  the  distance  in  inches  between  points  on  the 
Y-axis.  Enter  .1  number  not  a  word;  depress  the 
"Return"  key  on  keyboard. 

Enter  the  lottoi  unless  otherwise  instructed; 

Repress  the  "Return  key  on  the  keyboard. 

Enter  the  number  'I'.  Depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  size"  nt  the  number  to  be  displayed  along 
the  X  and  Y  axis.  Enter  a  number  between  0.07  and 
0.5  in  Inches.  Depress  the  "Return"  key  on  the 
keyboard. 

Type  in  the  title  lor  the  x-axis.  Limit  the  number 
of  characters  to  AO;  depress  the  "Return”  key  on 
the  keyboard. 

Enter  the  size  of  the  characters  for  the  X-axis 
title  in  inches.  Type  a  number  between  0.07  and 
0,5;  depress  the  "Return"  key  on  the  keyboard. 

Enter  the  number  depress  the  "Return”  key  on 

the  keyboard. 
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31)  Y-AXIS  TITLE  Il'AU  CHARS- IV  XXXXX_ 

3?)  CHARACTER  SIZE  IN  IIICHES  l  -0/~ -t> IV  XX_ 

33)  in  Black  IVI  Keo.  131  Breen 

lolor?  f_ 

34)  GRAPH  TITLE  (LINE  1)  11‘bU  CHAkCS- |V  XXXXX 

3b)  CHARACTER  SIZE  III  INCHES,  j  ._U/" -5 IV  XX_ 

36)  111  Black  I VI  Ked  131  Breeii 

LolorV  J_ 

3H  GRAPH  title  (line  V)  ll-UU  chars.  IV  xxxxx_ 


3«)  CHARACTER  SIZE  IN  INCHES  (  -U/~ -b  JV  XX 


UkVILHl 


Type  "in  the  1 1  to  for  the  Y-axis.  Limit  the  number 
of  characters  to  AO;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  size  ol  the  characters  for  the  Y-axis 
title  in  inches.  Inter  a  number  between  0.07  and 
0.5;  depress  the  “Return"  key  on  the  keyboard. 

Enter  the  number  '1';  depress  the  "Return"  key  on 
the  keyboard. 

Type  in  the  first  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  size  "of  the  characters  In  the  first  line 
of  the  graph  title  in  Inches.  Enter  a  number 

between  0.07  and  0.5;  depress  the  "Return”  key  on 
the  keyboard. 

Enter  the  number  'I*;  depress  the  "Return"  key  on 
the  keyboard. 

If  a  line  of  the  graph  title  is  to  be  centered 

under  a  preceding  line.  Include  an  appropriate 
number  of  blank  spaces  before  typing  the  title. 
Type  in  'the  second  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 
"Return”  key  on  the  keyboard. 

Enter  the  size  of  the  characters  In  the  second  line 

of  the  graph  title  in  inches.  Enter  a  number 

between  0.07  and  0.5;  depress  the  “Return"  key  on 
the  keyboard. 
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39)  II]  hack  121  Ked  131  Green 
Color?  L 

40)  GRAPH  title  (line  i )  []*bU  chars. I?  xxxxx 


41)  CHARACTER  SIZE  IN  INCHES  |.U/-.b|?  XX_ 


42)  111  Black  12]  Bed  13]  Creek 
Color?  L 

^3 )  graph  title  (line  9)  11-faO  charsI?  xxxxx_ 


44)  CHARACTER  SIZE  IN  INCHES  1 -0/- -b I V  XX_ 


4b)  111  Black  121  Bed  131  Green 

COLOR?  L 

46)  GRAPH  NUMBER  II 

INPUT  FILE  FOR  GRAPH  N  I FILE:'.!'!  <:|? 
XXXXXX.XXX 


Enter  the  lumber  '1';  depress  the  "Return"  key  on 
the  keyboard. 

Type  in  the  third  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 

"Return"  key  on  the  keyboard. 

Enter  the  size  of  the  characters  In  the  third  line 
of  the  graph  title  In  Inches.  Enter  a  number 
between  0.07  and  0.5;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  number  '1';  depress  the  "Return"  key  on 
the  keyboard. 

Type  in  the  fourth  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 

"Return"  koy  on  the  keyboard. 

Enter  the  size  of  the  characters  in  the  fourth 
line  of  the  graph  title  in  Inches.  Enter  a  number 
between  0.07  and  0,5;  depress  the  “Return"  key  on 
the  keyboard. 

Enter  the  number  '1';  depress  the  "Return"  koy  on 
the  keyboard. 

Enter  the  name  of  the  input  file  in  the  form 

'XXXXXX.XXX';  the  file  name  may  have  up  to  nine 
characters  before  the  period  and  must  have  three 
characters  after  the  period.  For  example, 
•CRI'lll.OOl ' ;  depress  the  "Return"  key  on  the 
keyboard. 
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47)  'actual  data'  curve. 

Ill  Black  121  Bed  151  Breen 
(.ot  or?  >_ 

48)  'estimated  values'  curve. 

Ill  Black  121  Bed  15)  Breen 
Color?  i_ 

49)  create  another  graph  IY/H1?  x_ 


M)  successful  comflet ion 


91)  >_ 


s?)  uwpusi.au) 


S3)  >*msK  tiRx  (si:  x_ 


54)  >*tlic  XXX.XX  (si:  XXX,XX_ 


55)  >*FILENAME.EXT  Isll  XXXXXX.XXX_ 


.MU  -miaillM 


awLtn 


Enter  the  number  *  l  * ;  depress  the  "Return'*  key  on 
the  keyboard. 


Enter  the  number  *  l';  depress  the  "Return"  key  on 
the  keyboard. 


Enter  either  the  le- ter  'V1  or  'N'i  depress  the 
"Return"  key  on  the  keyboerd. 

Computer  response.  H  no  other  graph  ip  to  he 
crested.  No  response  required  on  your  part.  If 
another  graph  Is  to  be  created,  the  computer  util 
return  to  step  T)  of  the  Erections. 

Computer  is  ready  to  plot  the  graph  on  the  Cslcoop 
plotter. 

Enter  *&>PMSC.Cm>* ;  depress  the  “Return"  key  on  the 
keyboard. 

Enter  th"  disk  -.umber  of  the  account  that  you 
a i giiod  on  the  terminal  with;  depresa  the  "Return* 
key  on  the  keyboard. 

E.  tr  the  UlC  number  of  the  account  that  you  signed 
on  the  terminal  with;  use  the  fora  XXX.XX:  depresa 
the  "Return"  key  on  the  keyboard. 

Enter  the  name  of  the  output  file  for  the  CRAEH. 
Use  the  for*  \mXXX.XXX\  This  will  be  th*  sane 
name  os  In  step  Hi)  ot  the  directions. 
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IKI  hlKt.t.ll  iUSPLAY 

%)  *  PlAlfi/UNED  1st:  Plain 

•  (SERIES  OF  COMPUTER  STATEMENls) 

s?)  >a  <tiiF  ^ 
sa)  >8«e_ 

"»*)  OUTPUT  Ft^e  i  FILE  SPEC  IV 

00!  f.KAPll  (IIWI0.R  (I 

INPUT  Fill  FOP  n#AP4«  N  I  HU  spec  I  V  XXX 


urTtt.ui 


Enter  the  word  'I't.iin'  unless  otherwise  Instructed; 
depress  the  "Return"  key  on  the  keyboard. 

No  response  required  on  your  |St  r  t . 


hast  ot  computer  cl ntements;  GRAPH  has  been  created 
and  stored. 

11  you  are  tlnlshed  us l hr  the  computer,  type  In 
*Hye‘  t  rum  the  keyboard;  depress  the  “Return"  key 
on  the  keyboard  and  you  will  he  automatically 
lopped  oil  I  In-  .  ouputer  terminal. 

The  GRAPH  can  he  picked  up  In  the  cuoputer  center, 
Uldi;  ISO,  in  appro* (mutely  1  hours. 

l>l rect  Inns  tor  option  I 

Enter  the  name  ut  the  output  file  In  the  fora 

‘  XXJtXXX.  XXJt  * ;  t  to*  tile  name  may  have  up  to  nine 

characters  he  (ore  the  period  aad  must  have  three 

eharaetets  dtit  the  period.  For  esaople, 
‘GRAPH  l,  util ' ;  depress  the  “return"  key  on  the 
keyboard. 

Enter  the  name  of  the  Input  file  In  the  form 

‘XXXXXX.XXN’;  the  file  name  may  have  np  to  nine 

characters  before  the  period  and  must  lutve  three 

eharaeters  after  the  period;  depress  the  “Return" 
key  on  live  Wyhuafd. 
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lumtm 


6U  CREATE  ANOTHER  GRAPH  lY/HlV  *_ 


6t)  SUCCESSFUL  COMPLETION 


Enter  the  letter  *Y'  or  * N * ;  depress  the  "Return" 
key  on  the  keyboard.  If  you  enter  the  letter  'Y1 
the  computer  goes  to  direction  9);  If  you  entor  the 
letter  *K *  the  computer  goes  to  62). 

Computer  response;  no  response  required  on  your 
part.  Computer  nuv  returns  to  direction  $1). 


ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


ASCC 

C 

Ca 

COE 

cm 

Ctt 

OB 

dc 

OTL 

Ein 

Eos 

Eout 

e* 

FG 

t 

*s 

tc 

i.d. 

i  rr 

K  or  k 
kfa 

UO 

vh 

vf 

vm 

vs 

vV 

mA 

min 


Air  Standardization  Coordinating  Committee 
capacitor 

ambient  challenge  atmosphere  concentration 
Chemical  Defense  Establishment  (United  Kingdom) 
centimeter 

sampled  leakage  concentration 

chemical  warfare • 

deep  breathing 

direct  current 

diode-transistor  logic 

input  voltage 

offset  voltage 

output  voltage 

exponential  of  l 

facial  grimacing 

current;  or,  the  ith  exercise 

Input  bias  current 

integrator  count 

inside  diameter 

current  leakage  through  the  integrator  capacitor 

international  Telephone  aftd  Telegraph 

kilo  (1000) 

kilo-Pascals 

light  emitting  diode 

logarithm  to  the  bat*  ten 

microampere 

microfarad 

micrometer 

microsecond 

microvolt 

mittlampere 

minute 

(Coat'd,  on  facing  page) 
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